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This  work  is  comprised  of  two  separate  investigations, 

both  related  to  spin-lattice  relaxation,  or  Tl,  information 

in  nuclear  magnetic  resonance  CNFIRD  imaging.  One  study 

explored  the  ability  of  commercially  available  NtlR  imagers 

to  accurately  and  precisely  determine  the  Tl  value  of  imaged 

objects.  The  specific  goal  was  to  evaluate  empirically  the 

advice  found  in  the  NflR  spectroscopy  literature  on  Tl 

determination,  and  how  well  this  advice  applied  to  NMR 

imaging  with  its  unique  set  of  experimental  constraints.  The 

primary  conclusion  was  that  if  one  wished  to  obtain  a  direct 

estimate  of  the  actual  Tl  value  in  an  object  which  might  be 

ix 


within  a  spatially  Inhomogeneous  radio  fraquancy  field,  the 
most  accurate,  precise,  and  time-effective  technique  to  use 
luas  three  fast  inversion  recovery  images,  with  suitably 
chosen  values  of  the  inverting  time,  whose  signal 
intensities  were  fitted  by  a  three-parameter  exponential 
function. 

The  other  study  concerned  itself  with  the  detailed 
theory,  practical  considerations,  and  possible  applications 
of  the  stimulated  echo  CSTED  in  NflR  imaging.  Whereas 
conventional  NMR  imaging  techniques  rely  upon  the  spin  echo, 
which  has  solely  TS  relaxation  weighting  unless  the  NMR 
signal  is  saturated,  STE  NMR  imaging  is  unique  in  that  the 
STE  has  intrinsic  Tl  weighting.  Possible  applications 
abound.  In  addition  to  generating  Tl  contrast  images,  it  is 
possible  to  calculate  quantitative  Tl  information  from  a 
series  of  STE  images.  Additionally,  STE  images  effectively 
enhance  objects  with  long  Tl  values  over  those  with  shorter 
Tl  values,  whereas  spin  echo  images  do  not.  Also,  it  was 
demonstrated  that  the  STE  easily  integrates  into  chemical 
shift  imaging  schemes.  Of  prime  interest  are  two  STE  imaging 
sequences  which  permit  the  acquistion  of  a  series  of  STE 
images  within  one  imaging  sequence,  where  each  image  has 
progressively  increased  Tl  weighting.  Finally,  a  method  of 
in  vivo  determination  of  diffusion  coefficients  is 
proposed,  which  utilizes  STE  imaging  to  lessen  the  effect  of 
T8  weighting. 


CHAPTER  I 
INTRODUCTION 


From  modest  but  promising  beginnings  in  the  1940s  and 
1950s,  nuclear  magnetic  resonance  CNMR^  spectroscopy  has 
developed  into  an  important  research  tool .  The  first 
applications  of  NMR  yielded  insights  into  the  properties  of 
the  atomic  nucleus.  Early  in  its  history,  NMR  spectroscopy 
uias  adapted  from  the  sole  use  by  physicists,  to  the  realm  of 
chemists,  tuho  saw  the  potential  of  the  chemical  shift 
phenomenon  as  a  structural  probe.  This  first  useful 
parameter  has  now  been  supplemented  by  many  other 
experimentally  accessible  quantities,  making  NMR  a  powerful 
and  versatile  technique  that  yields  information  related  to 
molecular  structure,  interactions,  and  dynamics.  Most  new 
applications  of  NMR  have  been  derived  from  parallel 
improvements  in  instrumentation  and  methods,  resulting  in 
applications  in  physics,  chemistry,  biology,  geology,  and 
medicine.  Particular  applications  include  two-dimensional 
fourier  transform  NMR  CAu76],  high  resolution  NflR  in 
solids  CHa7BD,  chemically  induced  dynamic  nuclear 
polarization  CKa78Il,  multiple  quantum  NflR  CUe773,  and 
NtIR  imaging  CLa733  .  The  conclusions  of  the  present 
investigation  are  particularly  applicable  to  NMR  imaging. 
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The  tBchniquB  of  NflR  imaaing  is  one  that  interests  not 
only  scientists,  but  also  nonscientists,  For  it  promises  to 
provide  a  safe  and  noninvasive  method  for  diagnosing 
dysfunction  or  disease  in  human  tissue.  Most  uhole  body  NMR 
imagers  in  use  today  detect  the  NMR  signals  from  the  protons 
in  the  object.  In  in  vivo  applications  the  major  source  of 
protons  is  water.  Many  applications  of  NflR  imaging  are  based 
on  the  fact  that  different  regions  of  the  object  Ci.o. 
different  tissues  or  organs)  have  different  uiater  contents 
and  different  values  of  the  characterizing  times,  Tl  and  12. 
For  soft  tissue  differentiation,  uihere  relative  uiater 
content  is  essentially  a  constant,  it  is  the  uida 
variability  in  the  values  of  Tl  and  12  uihich  makes  possible 
images  luith  anatomical  detail,  and  more  importantly,  images 
uhich  contain  pathological  detail.  These  images,  uthich 
contain  pathological  detail,  are  of  particular  interest 
because  of  the  proposal  CDa71D  that  Tl  values  of  water 
protons  in  cancerous  or  damaged  cells  are  longer  than  those 
of  protons  in  normal  calls  CLeBlD .  The  present 
investigation  explored  the  role  of  Tl,  the  spin-lattice 
relaxation  time,  in  NMR  imaging. 

Chapter  II  offers  a  brief  account  of  the  theory  of 
spin-lattice  relaxation  and  discusses  the  impetus  for  the 
present  investigation:  the  dual  role  of  Tl  in  NflR  imaging. 
One  role  of  Tl  in  NMR  is  qualitative  in  nature,  Tl  as  a 
source  of  image  contrast.  The  other  role  is  quantitative, 
for  Tl  values  have  the  potential  of  being  diagnostic. 
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Chapter  III  explores  the  ability  of  commercially 
available  NdR  imagers  to  accurately  and  precisely  determine 
the  Tl  value  of  imaged  objects.  The  specific  goal  uias  to 
evaluate  empirically  the  advice  on  Tl  determination  found  in 
the  NMR  spectroscopy  literature  and  to  see  hoiu  uiell  it 
applied  to  NMR  imaging  uiith  its  unique  set  of  experimental 
constraints.  The  primary  result  presented  in  Chapter  III  is 
concerned  with  the  direct  estimate  of  the  actual  Tl  value  in 
an  object  uihich  may  be  uithin  a  spatially  inhomogeneous  rf 
field.  The  most  accurate,  precise,  and  time  effective 
technique  to  use  is  three  fast  inversion  recovery  images, 
uiith  suitably  chosen  values  of  the  inverting  time,  uihose 
signal  intensities  are  fitted  by  a  three-parameter 
exponential  function. 

Chapter  lU  presents  the  detailed  theory,  practical 
considerations,  and,  by  way  of  examples,  the  possible 
applications  of  the  stimulated  echo  CSTE)  CHaSOD  in  NMR 
imaging.  Uhereas  conventional  NflR  imaging  techniques  rely 
upon  the  spin  echo,  STE  NtlR  imaging  incorporates  the  unique 
Tl  dependence  of  the  STE.  The  results  reported  in  Chapter  lU 
take  the  form  of  specific  applications  of  the  STE  to  NflR 
imaging.  First,  it  is  shown  that  in  addition  to  generating 
Tl  contrast  images,  it  is  possible  to  calculate  quantitative 
Tl  information  from  a  series  of  images.  Second,  a  novel 
application  of  the  Tl  weighted  STE  image  is  demonstrated, 
either  the  enhancement  or  the  suppression  of  elements  in  the 
object  with  different  Tl  values.  Third,  it  is  shown  that  the 
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5TE  is  easily  intsgratBd  into  chemical  shift  imaging 
schemes.  Fourth,  two  stimulated  echo  imaging  methods  are 
presented  which  permit  the  acquisition  of  a  series  of  5TE 
images  within  one  imaging  sequence,  where  each  image  has 
progressively  increased  11  weighting.  Finally,  a  method  of 
in  vivo  determination  of  molecular  translational 
self-diffusion  coefficients,  which  utilizes  STE  imaging  to 
lessen  the  effect  of  TS  weighting,  is  proposed. 

In  the  final  chapter.  Chapter  U,  a  summary  of  results 
is  presented,  as  well  as  the  final  conclusions,  which 
identify  the  thread  which  runs  through  this  entire  work,  a 
thread  that,  at  once,  connects  and  binds. 


CHAPTER  II 
UTILIZATION  OF  THE  SPIN-LATTICE  RELAXATION  TIME 

IN  NHR  IMABINB 


Introduction 

The  two  NflR  relaxation  times,  Tl  and  T2,  play  a 
pivotal  role  in  the  understanding  of  the  organization  of 
biological  systems  at  the  molecular  level  in  general,  and  in 
NriR  imaging  in  particular.  Differences  in  proton  NMR 
relaxation  times  of  normal  and  pathological  tissue  are  the 
key  to  NflR  image  contrast  and  the  discrimination  of  disease, 
a  fact  responsible  for  their  widespread  use  as  diagnostic 
parameters  in  clinical  NflR  imaging.  Tl  and  TB  directly 
affect  the  selection  of  imaging  pulse  sequence  timing 
parameters,  and  consequently,  the  total  imaging  times  and 
patient  throughput.  Also,  Tl  may  influence  the  choice  of  the 
optimum  operational  magnetic  field  strength  for  NflR  imaging, 
due  to  Tl's  significant  variation  with  frequency  CBo84D . 
Theory  of  Spin-Lattice  Relaxation 

A  detailed  description  of  the  complete  theory  of 
spin-lattice  relaxation  is  beyond  the  scope  of  this  work, 
rather,  those  elements  of  the  theory  crucial  to  the 
understanding  of  this  investigation  are  presented.  Consider 
an  ensemble  of  identical,  interacting  atoms,  whose  nuclei 
contain  either  an  odd  number  of  protons,  an  odd  number  of 
neutrons,  or  both.  At  room  temperature,  the  majority  of 
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nuclei  mill  reside  in  the  loiuest  energu  level,  the  so  called 

ground  state.  In  the  absence  of  a  magnetic  field  the  nuclear 

spin  states  are  degenerate  all  with  the  same  energy,  though 

the  application  of  a  magnetic  field  removes  this  degeneracy. 

Additionally,  the  nuclei  uill  precess  about  the  direction  of 

the  applied  magnetic  field,  a  concept  uhich  is  more  fully 

explained  in  Chapter  lU.  Radiofrequency  electromagnetic 

radiation  stimulation  uiill  cause  the  nuclei  to  absorb 

energy,  raising  them  to  an  excited  state.  The  nuclei  in  an 

excited  state  can  return  to  the  ground  state  only  by 

dissipating  the  excess  energy  to  their  surroundings.  Return 

to  the  ground  state,  or  relaxation,  also  requires  a 

stimulating  rf  field.  The  fields  causing  this  spin-lattice 

relaxation  are  provided  by  the  surrounding  nuclear 

environment,  the  so-called  lattice.  The  term,  lattice,  was 

first  introduced  to  describe  the  positions  of  molecules  in 

crystalline  solids,  but  has  since  been  extended  to  other 

phases  in  addition  to  the  solid  phase,  and  nou  simply 

indicates  the  magnetic  environment  of  the  nuclei . 

The  rf  fields  provided  by  the  lattice  for  relaxation 

result  from  the  presence  of  other  magnetic  nuclei, 

paramagnetic  ions  and  moleculBS,  and  molecular  magnetism 

uihich  is  the  result  of  the  fast  rotation  of  electronic 

charges.  The  most  common  source  of  lattice  fields  is  the 

dipole  field  produced  by  neighboring  magnetic  nuclei .  For 

example,  in  the  uiater  molecule,  one  of  the  hydrogen  nuclei 

produces  a  magnetic  field,  thus  affecting  the  adjacent 


7 
proton.  The  lattice  Fiald  must  fluctuate  to  transfer  energy 
effectively  from  the  excited  proton  to  the  lattice.  Thus, 
these  fluctuations  must  occur  at  a  rate  uihich  matches  the 
transitional  frequency  of  the  excited  protons. 

In  liquids,  the  fluctuations  in  the  lattice  field  are 
the  result  of  molecules  undergoing  Brouinian  motion,  uihich 
may  be  either  translational  or  rotational  in  nature.  Both 
intramolecular  and  intermolecular  relaxation  processes 
occur.  Uith  intramolecular  relaxation,  energy  is  transferred 
betiueen  nuclei  uiithin  the  same  molecule,  uihereas 
intermolecular  relaxation  involves  nuclei  of  different 
molecules.  The  protons  in  water  and  lipid,  the  primary 
sources  of  protons  in  the  human  body,  relax  predominantly  by 
the  intramolecular  dipole-dipole  mechanism. 

Typically,  the  average  rate  at  which  the  molecules 
reorient  themselves  is  related  to  the  size  of  the  molecule. 
Small  molecules,  such  as  water,  reorient  more  quickly  than 
larger  molecules,  such  as  lipids,  with  correlation  times  on 

-11      -B 

the  order  of  10    and  10   seconds  respactively .  Indeed,  the 

large  macromolecules,  such  as  DNf^  and  proteins,  tumble 
rather  slowly,  with  correlation  times  three  or  four  orders 
of  magnitude  slower  than  lipid.  The  frequency  of  rotation 
for  the  medium-sized  molecules,  such  as  lipids,  most  closely 
corresponds  to  the  transitional  frequency  of  the  excited 
protons,  at  typical  nuclear  magnetic  resonance  magnetic 
field  strengths.  Hence,  lipid-based  protons  will  relax 
faster  than  watei — based  protons,  which  rotate  at  a  frequency 
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that  is  typicallu  Qreatsr  than  the  transitional  frequency  of 
the  protons.  Similarly,  the  macromaleculos  are  inefficient 
in  causing  relaxation,  for  they  rotate  at  frequencies  uihich 
are  much  less  than  the  transitional  frequency. 

Efficient  relaxation  correlates  to  short  Tl  values, 
uihereas  inefficient  relaxation  results  in  long  Tl  values. 
For  example,  in  fat,  uihich  has  a  high  lipid  content,  Tl  is 
typically  of  the  order  of  a  feui  hundred  milliseconds,  yet 
the  Tl  of  pure  uiater  is  about  three  seconds.  Although  free 
uiater  relaxes  slouily,  the  uater  in  biological  tissue  tends 
to  relax  much  faster,  luith  typical  Tl  values  of  only  several 
hundred  milliseconds.  In  an  attempt  to  explain  this 
phenomenon,  it  uias  postulated  that  a  fraction  of  uiater  in 
tissues  is  bound  to  the  surface  of  proteins  CZi57D .  Hence, 
the  motion  of  bound  uater  is  reduced,  thus  more  closely 
matching  the  transitional  frequency  of  the  protons.  This 
enhanced  relaxation  results  in  shorter  Tl  values.  In 
practice,  an  equilibrium  exists  betuieen  bound  and  free 
water.  It  is  thought  that  this  equilibrium  is  perturbed  in 
certain  pathological  conditions,  resulting  in  the  clinically 
observed  elevation  in  Tl  values  of  certain  tumors  CDa71D . 
Clinical  Use  of  Spin-Lattice  Relaxation 

Although  spin  density  images  are  useful  in  a  number  of 
clinical  situations,  there  is  nearly  universal  agreement 
that  images  which  depend  significantly  on  the  relaxation 
parameters  Tl  and  TS  show  considerably  greater  soft  tissue 
contrast.  This  increased  contrast  allows  improved 
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differentiation  and  recognition  of  anatomical  detail  and 
helps  to  demonstrate  and  assess  mass  effects  CStSSD . 

The  explanation  for  this  increased  contrast  lies  in  the 
fact  that  relaxation  times  of  tissue  cover  a  uider  range  of 
values  than  the  range  of  proton  densities  in  the  same 
tissues.  For  example,  at  low  magnetic  field  strengths,  the 
difference  in  Tl  between  liver  and  kidney  is  over  505i, 
whereas  the  difference  in  proton  density  is  less  than  105s. 
Perhaps  more  important  is  the  large  alteration  in  relaxation 
time  that  occurs  in  various  disease  states,  even  when  the 
proton  density  itself  is  not  altered  significantly.  For 
example,  tumors  often  have  Tl  values  which  are  increased  by 
EOQyi   or  300/S  compared  to  the  surrounding  normal  tissue. 

The  determination  of  absolute  numbers  for  Tl  promises 
to  aid  in  tissue  specification,  although  there  is  still 
considerable  variability  of  values  being  reported  from 
investigator  to  investigator.  Recently,  Bottomley  at  al . 
have  collated  a  vast  body  of  relaxation  time  information  for 
the  purpose  of  establishing  the  range  of  normal  values 
CBoBHD .  Deficiencies  in  measurement  techniques  were 
identified  as  a  major  source  of  data  irreproducibility . 
Additionally,  there  is  significant  overlap  between  normal 
and  abnormal  tissue  in  some  types  of  pathology,  for  example, 
diffuse  liver  disease  CDoBBD .  Thus,  to  date,  measured  Tl 
values  are  still  too  variable,  and  have  yet  to  improve 
significantly  the  tissue  specificity  of  NMR  imaging. 


CHAPTER  III 
INUESTIGATION  INTO  Tl  DETERniNATION  ON  A 
WHOLE  BODY  NHR  IMAGER 


Introduction 

Proton  nuclear  magnetic  resonance  CNMR)  imaging  may 
yield  both  qualitative  images  uihich  are  evaluated  visually 
and  quantitative  information  uihich  is  evaluated  numerically. 
The  quantitative  information  may  be  used  to  generate  values 
of  localized  in  vivo  NUR  parameters  such  as  the  spin 
density,  the  spin-spin  relaxation  time,  TE,  and  the 
spin-lattice  relaxation  time,  Tl .  The  NflR  spectroscopy 
literature  is  full  of  suggestions  on  hotu  best  to  determine 
these  parameters,  particularly  Tl,  in  conventional  samples 
studied  by  physicists  and  chemists  CGrB3,  HaBlD . 
Similarly,  the  NflR  imaging  literature  offers  a  multitude  of 
suggestions  for  in  vivo  Tl  determination,  based  on  both 
empirical  and  theoretical  arguments  CPy83D .  The  goal  of 
this  work  was  to  evaluate  empirically  the  advice  found  in 
the  NflR  spectroscopy  literature  as  to  how  well  it  applied  to 
NMR  imaging,  with  its  unique  experimental  constraints. 

The  pulse  programming  capabilities  of  commercial  NflR 
imagers  are  typically  limited,  yet  there  are  no  fewer  than 
four  distinct  pulse  sequences  available  which  may  be  used  to 
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provide  estimates  of  the  actual  in  vivo  Tl:  the  spin  echo 
CSE)  sequence  CHaSOD,  the  progressive  saturation  CPS) 
sequence  CFr71D,  the  inversion  recovery  CIR)  sequence 
CUoEBD,  and  the  fast  inversion  recovery  CFIR)  sequence 
CKo773.  Each  of  these  sequences  has  at  least  tuo  timing 
parameters,  values  of  uhich  need  to  be  optimized  so  as  to 
reduce  the  overall  error  in  the  determined  Tl  value.  In 
addition,  a  rational  choice  of  the  number  of  images  to  be 
generated,  the  spacing  of  the  variable  timing  parameters, 
and  the  form  of  the  fitting  function  must  be  made  with 
consideration  toutards  the  resultant  precision  and  accuracy 
in  measurement,  and  the  net  imaging  time  required.  Precision 
is  used  in  this  work  to  mean  the  degree  of  exactness  with 
which  a  quantity  is  stated,  while  accuracy  is  used  to 
indicate  the  conformity  of  an  indicated  value  to  an  accepted 
standard  value  CUa753 .  A  precise  Tl  estimate  is  not 
necessarily  accurate,  whereas  it  is  hoped  that  all  accurate 
Tl  estimates  are  precise.  Furthermore,  other  constraints  of 
the  system  must  be  incorporated  into  the  decision  process, 
namely  the  spatially  inhomogeneous  rf  field,  which  is  the 
result  of  rf  coil  design  and  the  rf  attenuating  properties 
of  the  object  CBo78D. 

In  the  present  study,  experiments  were  conducted  on  a 
phantom  consisting  of  an  array  of  vials  containing 
paramagnetically  doped  water,  whose  Tl  values  spanned  the 
range  of  clinical  concern.  Each  previously  mentioned 
parameter  was  systematically  varied,  and  single  slice  images 
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of  the  phantom  uere  generated  uihlle  the  phantom  uas  First 
immersed  in  air,  and  second,  in  a  saline  solution  of 
physiological  concentration.  The  effect  of  the  saline 
solution  was  to  dissipate  the  rf  energy  in  a  similar  fashion 
to  the  human  body,  resulting  in  a  spatially  inhomogeneous  rf 
field.  The  raui  data  mere  computer  processed,  resulting  in  a 
large  population  of  reduced  data,  upon  uihich  correlative 
studies  mere  performed . 

This  investigation  produced  two  principal  conclusions. 
The  primary  conclusion  was,  accepting  the  assumptions  given 
in  the  methods  and  materials  section  concerning  the  phantom 
material,  that  to  obtain  a  direct  estimate  of  the  actual  Tl 
value  in  an  object  which  may  be  within  a  spatially 
inhomageneaus  rf  field,  the  most  accurate,  precise,  and  time 
effective  technique  to  use  was  three  FIR  images,  with 
suitably  chosen  values  of  the  inverting  time,  TI,  fitted  by 
a  three-parameter  exponential  function.  The  use  of  IR  images 
was  equally  accurate  and  precise,  but  not  as  temporally 
efficient,  i^  second  conclusion  was  that  the  use  of  three  SE 
images,  with  suitably  chosen  values  of  the  pulse  sequence 
repetition  time,  TR,  and  fitted  by  a  two-parameter 
exponential  function  can  be  more  precise  and  time  effective 
than  the  FIR  technique  in  estimating  Tl,  but  was  always  much 
less  accurate .  These  conclusions  may  be  used  a  priori  to 
design  the  series  of  images  which  will  yield  the  most 
accurate,  precise,  and  time  efficient  estimate  of  an  in 
vivo  Tl  value,  or  they  may  be  used  a  posteriori  to  help 
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evaluate  the  accuracy  of  a  calculated  Tl  value  From  a  given 
series  of  images.  Examples  of  each  type  of  application  are 
given  in  the  results  section. 

Methods  and  flatsrials 
NfIR  Imaginn  System 

All  experiments  in  this  investigation  were  conducted  on 
a  Teslacon  whole  body  NfIR  imager,  manufactured  by  Technicare 
Corporation.  The  main  magnetic  field  was  produced  by  a  six- 
coil,  air-core,  water-cooled  resistive  electromagnet.  The 
magnet  was  nominally  operated  at  0.15  Tesla,  but  could  be 
adjusted  by  altering  the  current  in  the  electromagnet.  The 
clear  bore  diameter  was  one  meter,  allowing  access  of  an 
entire  human  body,  with  the  long  axis  of  the  body  aligned 
with  the  z  direction  of  the  main  magnetic  field. 

The  imager  used  separate  rf  receiver  and  transmitter 
coils,  allowing  separate  performance  optimization.  The 
transmitter  coil  was  for  all  intents  and  purposes  a  part  of 
the  imager,  although  it  was  passible  to  utilize  a  myriad  of 
rf  receiver  coils.  The  manufacturer  supplied  body  rf  coil 
was  used  for  all  measurements.  This  permitted  the 
investigation  to  be  conducted  over  a  large  field  of  view, 
0.75  meters,  which  facilitated  the  examination  of  spatially 
dependent  variables. 

The  rf  coil  used  in  NMR  is  part  of  a  resonant  circuit, 
and  it  was  possible  to  measure  the  quality  factor,  Q,  of  a 
circuit  containing  the  body  rf  coil,  on  the  bench.  A  sweep 
wave  generator,  manufactured  by  Uavetec,  was  coupled  to  the 
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rf  coil.  The  frequency  dependent  response  was  monitored  with 
an  oscilloscope,  and  the  Q  was  determined  by 

Q  -  f/6F  C3-1D 

where  6f  was  the  width  of  the  resonance  curve  at  70.7^  of 
the  amplitude  of  the  response  at  the  resonance  frequency,  f 
CKrBlD.  The  □  was  determined  with  the  coil  both  physically 
unloaded,  and  loaded  with  a  dielectric  material  (O.S/i  NaCl 
aqueous  solution) .  The  ratio  of  the  unloaded  Q  to  loaded  Q 
was  approximately  0.7. 

The  Teslacon  imager  made  use  of  a  DEC  PDP  11/24 
computer  for  pulse  sequence  programming,  data  acquisition, 
and  data  processing.  Additionally,  a  FPS  floating  point 
array  processor  aided  in  data  processing.  All  experiments 
were  conducted  with  the  standard  software  supplied  by  the 
manufacturer  Csof tware  release  C) ,  except  for  modifications 
implemented  which  allowed  for  continuous  data  acquisition 
with  software  controlled   incremental  timing  parameters. 
Phantom 

The  phantom  used  for  all  measurements  was  designed  to 
model  certain  properties  of  the  human  body .  Due  to  the 
spatial  extent  of  the  body  it  was  of  interest  to  examine  the 
possible  effects  of  spatially  dependent  parameters,  such  as 
main  magnetic  field  and  rf  magnetic  field  inhomogeneities , 
upon  in  vivo  spin-lattice  relaxation  measurements.  Hence, 
the  phantom's  geometric  design,  as  seen  in  figure  3-1, 
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Figure  3-1 


Phantom  used  for  all  Tl  measurements.  Inner 
structure  holds  eight  vials  coaxial  with  the 
main  magnetic  field  luhile  outer  structure  could 
be  filled  with  saline  solution. 
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permitted  comparative  measurements  to  be  made  on  identical 
samples  uihich  were  spatially  distributed  within  the  imager's 
Field  of  view.  The  external  container  could  be  filled  with  a 
saline  solution,  immersing  the  phantom. 

The  range  of  Tl  values  within  the  body  typically  spans 
from  100  milliseconds  to  1000  milliseconds,  with  a  few 
singular  exceptions  Ce.g.  cerebral  spinal  fluid).  As  given 
in  table  3-1,  all  experiments  were  conducted  on  aqueous 
solutions  of  copper  sulfate  with  varying  concentrations, 
spanning  a  Tl  range  of  70  to  1100  milliseconds. 

Three  assumptions  were  made  concerning  this  simple 
phantom  material.  First,  the  spin-lattice  relaxation  decay 
process  was  a  monotonically  decreasing  function  of  time. 
Second,  the  exponential  decay  had  a  single  time  constant. 
Third,  the  noise  spectrum  was  white,  Gaussian,  and  had  a 
zero  mean.  The  third  assumption  has  been  verified  for  the 
NMR  imager  over  a  wide  range  of  experimental  conditions 
CSaB4D . 

The  actual,  or  "gold  standard"  value  of  Tl  for  each  of 
the  phantom  materials  was  determined  in  the  following 
manner:  a  single  cylindrical  glass  vial  Cwith  a  diameter  of 
1.0  centimeters  and  a  length  of  10.0  centimeters) 
containing  the  doped  water  solution  under  analysis,  was 
placed  at  the  position  of  maximum  sensitivity  of  the  rf 
transmission  coil,  and  coaxial  to  the  main  magnetic  field. 
Owing  to  the  small  filling  factor  resulting  from  the 
relatively  small  sample  volume,  the  receiver  gain  was 
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Table  3-1 

Concentration  of  Copper  Sulfate  Doped  Uater  and 
Resultant  Spin-Lattice  Relaxation  Time  for  Phantom  Material 


Material  Concentration  Cmfl)            Tl  Cmsec) 

BCKey  1)                   5.0  73.E±14.3 

CCKey  E)                   3.5  154.0±13.a 

DCKey  3)                   B.O  E77.B±E3.1 

ECKey  4)                   1.0  540.4±57.1 

FCKay  5)                   0.5  110E.E1109.3 
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adjusted  to  make  optimum  use  of  the  dynamic  range  of  the 
analog-to-digital  converters.  The  imager  was  set  up  to 
acquire  an  IR  image,  luith  the  pulse  sequence  repetition 
time,  TR,  set  to  10  times  the  expected  Tl  value  so  as  to 
avoid  saturating  the  signal.  All  magnetic  field  gradients, 
required  for  spatial  encoding  of  the  NMR  signal,  mere  turned 
off.  This  resulted  in  the  imager  being  used  as  a 
conventional  spectrometer.  The  inverting  time,  TI ,  ujas 
incrementally  varied  for  over  twenty  values  uihich  spanned 
the  estimated  Tl  value.  For  each  value  of  TI,  the  Fourier 
transform  of  the  time  domain  signal  uias  observed,  and  ten 
values  of  maximum  amplitude  and  mean  noise  uere  recorded. 
The  arithmetic  means  of  these  values  were  fitted  with  a 
three-parameter  exponential  function  using  the  data 
reduction  routine  of  the  program  "NUR",  resident  on  a 
Nicolet  1180E  computer  CNmBSD .  This  procedure  was  repeated 
for  each  different  doped  solution  on  two  occasions, 
separated  by  over  a  year,  with  the  results  given  in  table 
3-1. 

In  addition  to  effecting  the  □  of  the  rf  coil,  a 
dielectric  substance  will  also  dissipate  incident  rf  energy. 
An  otherwise  spatially  homogeneous  rf  field  will  become 
spatially  inhomogeneous  if  a  dielectric  fills  the  space. 
Hence,  the  mere  presence  of  the  body  in  the  NflR  imager  may 
alter  the  homogeneity  of  the  transmitted  rf  field.  The 
macroscopic  implication  of  this  effect  was  a  variation  in 
pulse  tip  angle  from  point  to  point  within  the  body.  For 
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example,  for  a  given  transmitted  rf  pulse,  a  region  close  to 
the  surface  of  the  body  might  experience  a  IBO  degree  tip 
angle,  while  an  interior  region,  owing  to  rf  attenuation  by 
the  tissue,  might  experience  a  170  degree  tip  angle.  Since 
Tl  measurement  accuracy  is  often  sensitive  to  misset  tip 
angles,  the  problem  of  a  spatially  inhomogeneous  rf  field 
was  studied  as  part  of  this  investigation.  Figure  3-1  shows 
the  phantom  holder  positioned  within  a  large  cylindrical 
bottle.  All  experiments  were  conducted  on  the  phantom  within 
the  empty  bottle,  and  also  within  the  bottle  while  full  of 
saline  solution  of  physiological  concentration  C0.9/i  NaCl). 
The  dimensions  of  the  battle  approximated  those  of  a  human 
abdomen . 
Spatial  flapping  of  Intrinsic  Inhomogeneities 

Tl  is  a  calculated  value,  determined  from  intensity 
measurements.  ThereforB  any  source  of  intensity  variation 
might  introduce  errors  into  Tl  calculations.  Two  common 
sources  were  main  magnetic  field  inhomogeneity  and  rf  field 
inhomogeneity .  Before  attempting  any  experiments,  the 
intensity  of  the  main  magnetic  field  and  rf  transmitted 
field  were  spatially  mapped  utilizing  the  specially 
constructed  field  mapper  depicted  in  figure  3-2.  The  results 
are  presented  in  figures  3-3  and  3-4  respectively . 

Since  all  images  were  single  slice,  it  was  sufficient 
to  examine  the  z  equal  zero  plane  for  main  magnetic  field 
and/or  rf  transmitted  field  inhomogeneities.  A  56  centimeter 
diameter  circular  piece  of  lucite  with  thirty-seven  small 
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Figure  3-2.  Field  mapping  apparatus,  flultiple  point  source 

phantom  with  specially  designed  rf  receiver  coil 
shoun  at  center.  Designed  for  measurements  in 
the  transverse  plane  over  a  56  centimeter 
diameter  circular  region. 
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Figure   3-3 


Spatial  hamogenBity  mapping  of  the  main  magnetic 
field.  Each  isastrength  curve  is  S.5  ppm  of  the 
main  magnetic  field,  fleasurements  were  made  in 
the  transverse  plane  over  a  56  centimeter 
diameter  circular  region. 
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Figure  3-4 


Spatial  homogeneity  mapping  of  the  transmitted 
rf  field.  Each  isostrength  curve  is 
approximately  0  .ES^i   of  the  field  at  the  central 
region.  Heasurements  were  made  in  the  transverse 
plane  over  a  56  centimeter  diameter  circular 
region. 
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vials  of  doped  utater  embedded  ulthln  It  uas  placed  in  the  z 
equal  zero  plane.  A  one  centimeter  diameter  rf  receiver  coil 
uas  constructed  ,  as  seen  in  figure  3-E,  which  could  be 
placed  in  turn  over  each  small  vial.  The  imager  tuas  operated 
as  a  spectrometer,  uith  the  only  detectable  source  of  signal 
being  the  vial  over  uhich  the  special  rf  coil  uas  situated, 
resulting  in  essentially  point  source  measurements. 

To  evaluate  the  main  field  inhomogeneity ,  the  coil  was 
placed  on  the  center  vial,  and  the  main  magnetic  field  luas 
adjusted  to  insure  resonance.  The  coil  uias  then  moved 
systematically  from  vial  to  vial  while  the  registered 
deviation  from  the  resonance  frequency  was  recorded.  Since 
no  imaging  gradients  were  in  use,  there  was  a  one-to-one 
correlation  between  resonance  offset  and  main  field 
inhomogeneity.  Over  the  useful  field  of  view,  a  coaxial 
circle  of  thirty  centimeters  in  diameter,  the  main  magnetic 
field  was  homogeneous  to  approximately  B5  parts  per  million 
of  the  main  field.  This  value  was  within  the  manufacturer's 
specifications  and  was  deemed  experimentally  acceptable. 

Although  this  investigation  considered  the  effects  of  a 
spatially  inhomogeneous  rf  field,  it  was  concerned  primarily 
with  spatial  variations  resulting  from  rf  attenuation  by  the 
object,  and  not  intrinsic,  and  therefore  constant 
inhomogeneities  of  the  transmitted  rf  field.  This  intrinsic 
variation  was  determined,  for  documentation  purposes,  in  the 
following  manner. 
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The  special  rf  coil  was  sustematically  moved  From  vial 
to  vial,  and  while  on  each  vial  the  main  magnetic  field 
would  be  slightly  altered  so  as  to  insure  resonance.  This 
was  accomplished  simply  by  altering  the  current  in  the 
electromagnet.  With  the  resonance  condition  holding, 
numerous  intensity  measurements  were  recorded  and  later 
averaged.  This  procedure  was  repeated  For  each  vial.  After 
all  intensity  values  had  been  recorded,  the  averaged 
intensities  were  compared  to  determine  the  transmitted  rF 
homogeneity.  IF  each  identical  point  source  sample 
experienced  the  same  tip  angle,  then  all  intensity  values 
should  have  been  equal  on  resonance.  Any  variation  in 
intensity  was  attributed  to  variation  in  tip  angle.  Dver  the 
useFul  Field  oF  view,  a  coaxial  circle  about  thirty 
centimeters  in  diameter,  the  transmitted  rF  Field  varied  by 
approximately  3.0^. 
Data  Acquisition  flethodologu 

The  investigation  was  aimed  at  determining  the  most 
eFFicient  and  accurate  method  of  obtaining  actual  Tl  values 
from  NflR  images,  with  consideration  towards  the  constraints 
a  spatially  inhomogeneous  rf  field  imposses .  Nonuniform  rf 
irradiation  was  a  concern  for  it  results  in  a  spatially 
dependent  systematic  error  in  rf  pulse  tip  angles,  and  hence 
in  measured  Tl  CFr71D .  Evaluation  of  four  pulse  sequences 
was  conducted,  the  spin  echo  CSE),  the  inversion  recovery 
CIR),  the  progressive  saturation  CPSD,  and  the  fast 
inversion  recovery  CFIR) .  The  degree  of  freedom  of  the 
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fitting  function  dictated  the  Iciuer  bound  on  the  number  of 
images  required  for  a  unique  determination  of  the  Tl  value, 
bjhile  the  upper  bound  uias  investigated  as  to  its  dependence 
upon  pulse  sequence  utilized,  desired  accuracy  in 
measurement,  and  net  imaging  time.  Similarly,  the  spacing  of 
the  variable  timing  parameter  ujas  evaluated  as  a  function  of 
the  measurement  accuracy  and  net  imaging  time  for  a  given 
pulse  sequence. 

For  all  experiments,  a  nonlinear  least  squares  fitting 
algorithm  was  used,  of  the  form 

SCTP)  -  K  -  CCexpC-TP/TDD  C3-2D 

luhere  K  and  C  uiere  constants  and  SCTP)  was  the  signal 
intensity  as  a  function  of  a  timing  parameter,  TP . 

In  PS  and  SE  images,  Tl  weighting  was  introduced  by 
saturating  the  signal  with  rapid  pulsing.  Therefore,  TP  was 
the  pulse  sequence  repetition  time,  TR .  The  spectroscopy 
literature  CFr71D  suggests  using  a  fitting  function  with 
two  degrees  of  freedom,  given  by  setting  C  in  equation  3-S 
equal  to  K,  to  determine  Tl  from  either  PS  or  SE  images. 

Tl  information  was  incorporated  into  IR  images  by 
inverting  the  equilibrium  magnetization,  allowing  some  time 
to  pass  during  which  spin-lattice  relaxation  occurred,  and 
sampling  the  remaining  magnetization  by  bringing  it  into  the 
transverse  plane  where  detection  took  place.  For  IR  images, 
TP  was  equal  to  the  inverting  time,  Tl  .  The  NflR  spectroscopy 
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literature  CUoBBD  suggests  using  a  Fitting  function  with 
two  degrees  of  freedom,  given  by  setting  C  in  equation  3-E 
equal  to  2K,  for  IR  images  where  the  inverting  pulse  was 
exactly  IBO  degrees.  When  the  inverting  pulse  was  misset, 
possibly  due  to  a  spatially  inhomogeneous  rf  field,  the 
literature  CKo773  suggests  the  use  of  a  fitting  function 
with  three  degrees  of  freedom,  such  as  equation  3-S . 

The  FIR  images  were  identical  to  the  IR  images,  except 
that  a  rapid  IR  was  used,  with  one  result  being  a  reduction 
of  the  total  image  time.  Typically,  TR  for  IR  images  was 
five  times  Tl  or  longer,  while  in  FIR  images  TR  was  commonly 
two  to  three  times  Tl .  As  suggested  in  the  literature 
CCa75D,  this  rapid  TR  requires  the  use  of  a  fitting 
function  with  three  degrees  of  freedom,  such  as  equation 
3-E,  whether  the  inverting  pulse  was  misset  or  not. 

Details  of  the  exact  experimental  procedure  follow.  It 
should  be  noted  that  the  specified  method  was  performed  for 
each  of  the  five  phantom  samples  as  given  in  table  3-1,  both 
surrounded  by  air  and  surrounded  by  the  dielectric  saline 
solution.  This  permitted  the  analysis  to  span  the  entire 
clinically  useful  Tl  range,  and  also  to  evaluate  the  effects 
of  misset  rf  pulse  tip  angles  resulting  From  an 
inhomogeneous  rf  field.  All  acquired  images  were  single 
slice,  at  the  z  equal  zero  plane,  with  slice  thickness  of 
one  centimeter.  The  Teslacon  imager  typically  gathered  256 
data  points  in  the  readout  direction  and  128  data  points  in 
the  phase  encoding  direction.  The  displayed  image  was  512 
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pixels  bu  512  pixels  which  uias  derived  from  the  stored  image 
data  uihich  uias  dimensioned  256  by  256.  Hence  the  data  in  the 
phase  encoding  direction  were  interpolated.  To  avoid  the  use 
of  interpolated  data,  the  data  acquisition  routine  was 
modified  to  permit  the  use  of  256  phase  encoding  gradient 
steps.  In  all  experiments  the  echo  time,  TE,  was  maintained 
at  its  shortest  value,  30  milliseconds,  to  minimize 
contributions  from  spin-spin  relaxation  and  molecular 
self-diffusion.  Two  averages  were  taken  for  each  image. 

The  PS  and  SE  experiments  differed  in  only  two 
respects.  First,  the  PS  images  relied  upon  the  bulk 
magnetization's  reaching  a  steady-state  value  in  the 
presence  of  the  rapid  pulsing  with  Tt/E   rf  pulses.  It  is 
suggested  that  the  system  reaches  this  steady-state  within 
four  pulses  CFr713,  hence  all  PS  images  were  preceded  by 
four  pulses  prior  to  image  acquisition.  Second,  PS 
experiments  yield  consistent  values  of  Tl  for  TR  in  the 
range  of  0.5T1  to  2.0T1,  whereas  SE  experiments  yield 
consistent  values  of  Tl  for  TR  in  the  range  of  0.5T1  to 
3.0T1.  Hence,  images  were  collected  with  TR  ranging  from  its 
minimum  value,  dictated  by  the  manufacturer  to  be  50ms,  to  a 
maximum  of  three  times  the  actual  Tl .  Additionally,  a  TR 
value  of  5.0T1  was  used  to  sample  the  unsaturated  initial 
magnetization.  In  total,  no  less  than  twenty  images  were 
acquired  with  varying  TR  values.  In  data  processing  images 
with  TR  less  than  or  equal  to  1.5T1  were  considered  PS 
images,  and  those  with  TR  greater  than  1.5T1  were  SE  images. 
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The  IR  and  FIR  experiments,  as  acquired  with  the 
imager,  differed  in  only  two  respects.  First,  the  IR  images 
bjere  acquired  assuming  that  the  bulk  magnetization  uias  at 
its  initial  equilibrium  value  prior  to  the  commencement  of 
the  pulse  sequence.  Therefore,  TR  for  all  IR  images  uas  set 
to  at  least  five  times  the  actual  Tl  value.  This  condition 
bjas  not  required  for  FIR  images,  hence  the  origin  of  "fast" 
in  fast  inversion  recovery,  with  TR  nominally  set  to  twice 
the  actual  Tl  value  in  the  FIR  images.  Second,  since  the 
maximum  value  of  Tl  must  be  less  than  TR,  the  range  of  Tl 
for  IR  images  was  from  its  minimum  value,  dictated  by  the 
manufacturer  to  be  25  milliseconds,  to  Just  less  than  5.0T1. 
In  the  FIR  images,  Tl  ranged  from  25  milliseconds  to  Just 
less  than  2.0T1.  In  total,  no  less  than  ten  IR  images  or  ten 
FIR  images  mere  acquired  with  varying  Tl  values. 
Data  Processing  flethodologu 

All  images  acquired  with  the  Teslacon  system  were 
displayed  on  a  high  resolution  monochromatic  CRT  monitor, 
with  a  maximum  of  1024  gray  levels.  It  was  also  possible  to 
display  the  images  on  a  high  resolution  color  monitor,  as 
shown  in  figure  3-5.  Each  image  was  reconstructed,  and  the 
spatially  dependent  signal  intensities  were  displayed  on  the 
CRT  monitor  for  analysis.  Reconstruction  of  the  PS  and  SE 
images  differed  from  that  of  the  IR  and  FIR  images. 

The  Imager  acquired  all  raw  data  using  quadrature 
detection,  effectively  resulting  in  each  data  point's  being 
defined  by  a  complex  number .  PS  and  SE  experiments  sampled 
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Figure  3-5, 


Color  NflR  images  of  phantom.  Ca)  With  eight 
identical  vials,  Cb)  with  all  vials  surrounded 
by  saline  solution. 
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the  magnetization,  uhose  value  lay  betuieen  zero  and  its 
equilibrium  value.  Since  this  is  a  one  sided  range,  all 
positive  in  sign,  signal  intensities  in  PS  and  SE  images 
were  represented  by  the  magnitude  oF  the  corresponding 
complex  number.   IR  and  FIR  experiments  sampled  the 
magnetization  uihile  it  uias  in  the  range  oF  plus  or  minus  its 
equilibrium  value.  To  retain  the  signed  InFormation,  signal 
intensities  in  IR  and  FIR  images  uiera  determined  and 
represented  taking  into  account  the  phase  oF  the 
corresponding  complex  number. 

tilhile  in  the  display  mode,  quantitative  InFormation  uias 
obtained  From  the  image.  Signal  intensity  InFormation  uas 
gathered  by  positioning  a  soFtware  controlled  region  oF 
interest,  RQI,  about  the  spatial  area  oF  concern.  The  system 
made  available  the  number  oF  pixels  enclosed  in  the  ROI,  the 
mean  intensity  value  oF  those  pixels,  and  the  corresponding 
standard  deviation.  For  all  images  taken  with  the  phantom 
designed  For  this  investigation,  no  Fewer  than  one  hundred 
pixels  were  used  to  deFined  each  vial .  The  intensity 
InFormation  From  all  eight  vials  in  the  Field  oF  view,  and 
also  the  intensity  InFormation  From  a  representative  region 
oF  background  noise,  were  determined  and  recorded.  These 
measurements  Formed  the  raw  data  base  on  which  all 
subsequent  analysis  was  conducted. 

A  Fortran  program,  which  perFormed  a  nonlinear  least 
squares  Fit  without  the  need  For  initial  guesses  oF  the 
Fitted  values,  was  written  and  implemented  on  a  IBM  470 
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computer  system.  The  general  purpose  of  this  program  uias  to 
determine  the  optimal  number  of  images,  and  the  values  of 
their  associated  timing  parameters,  required  to  accurately 
and  precisely  determine  Tl  for  a  given  pulse  sequence.  The 
details  of  hou  this  uas  accomplished  follou. 

The  program  accepted  as  input  the  type  of  pulse 
sequence  and  the  signal  intensities  as  a  function  of  the 
corresponding  timing  parameter.  This  information  uas  given 
for  each  vial,  for  each  different  concentration  of  solution, 
and  for  both  immersion  in  air  and  immersion  in  saline 
solution.  In  NfIR  spectroscopy,  where  the  sample  under 
investigation  may  be  examined  for  any  given  period  of  time, 
typically  ten  to  thirty  data  pairs  are  used  in  fitting  the 
empirical  data  to  the  theoretically  expected  function, 
resulting  in  an  estimate  of  Tl .  In  NfIR  imaging,  the  time 
constraint  is  more  restrictive.  Ill  patients  cannot  remain 
in  the  imager  for  an  extended  period  of  time,  and  the 
physician  is  usually  not  willing  to  spend  an  inordinate 
amount  of  time  on  a  single  procedure. 

This  study  aimed  to  determine  whether  the  acquisition 
of  two,  three,  or  four  images  resulted  in  the  most  accurate 
and  precise  Tl  estimate.  If  each  data  set  contained,  say, 
ten  data  pairs,  the  analysis  program  considered  each 
possible  combination  of  two  data  pairs  at  a  time,  three  at  a 
time,  and  four  at  a  time.  These  two,  three,  or  four  data 
pairs  were  then  fitted  to  the  appropriate  functional  form  of 
equation  3-B  for  the  given  pulse  sequence. 
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If  a  function  with  three  degrees  of  frBedom  luas 
required,  then  only  the  combinations  involvino  three  and 
four  data  pairs  mere  used.  For  example,  an  IR  data  set, 
consisting  of  ten  data  pairs  and  fitted  to  both  a  function 
with  two  degrees  of  freedom  and  a  function  with  three 
degrees  of  freedom,  generated  seven  hundred  and  five  fitted 
values  of  Tl .  Multiplying  this  by  the  eight  vials,  the  five 
different  concentration  samples,  and  the  passible  immersion 
in  air  or  saline  resulted  in  over  five  thousand  estimated  Tl 
values.  These  estimated  values  were  compared  to  the  actual 
Tl  value,  and  also  to  the  estimated  Tl  value  obtained  if  all 
of  the  data  pairs,  for  example  ten,  were  fitted  to  equation 
3-2.  Comparison  with  the  actual  Tl  tested  for  precision  and 
accuracy,  while  comparison  with  the  estimated  Tl  tested 
simply  for  precision.  It  required  over  five  minutes  of  CPU 
time  to  process  all  the  raw  data  and  transfer  in  excess  of 
twelve  megabytes  of  reduced  data. 

Table  3-2  presents  some  typical  raw  data  which  were 
acquired  for  Tl  determination.  The  experiment  performed  was 
an  inversion  recovery  experiment,  conducted  on  phantom 
material  C.  The  TR  and  TE  values  were  kept  constant  at 
1250msec  and  SOmsec  respectively.  Data  are  presented  for 
both  the  phantom  material  immersed  and  not  immersed  in  the 
saline  solution.  Column  one  depicts  the  values  of  Tl  used 
for  Tl  determination,  chosen  to  properly  span  the  actual  Tl 
value  of  phantom  material  C  of  154msBC.  The 
region-of-interest  derived  signal  intensities  for  vials  1 
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Table  3-2 
Typical  Raw  Data  Acquired  for  Tl  DBtermination 


TI CmsBcD 


SICvial  1) 


SICvial  5D 


Background 


PHANTOn  NOT  inHERSED  IN  SALINE  SOLUTION 


50 
75 
1S5 
175 
2B5 
275 
325 
400 
1000 


-43.7 

-7.1 

44.1 

B7.4 

113.5 

14B.7 

153.6 

186.2 

223.3 


-90.4 
-26.4 
BO.l 
173.2 
229.5 
273.9 
326.0 
361.4 
430.0 


0.38 

0.12 

0.06 

-0.26 

-0.27 

-1.32 

1.15 

-0.07 

0.05 


PHANTOn  innERSED  IN  SALINE  SOLUTION 


50 
75 
125 
175 
225 
275 
325 
400 
1000 


16.3 
21.8 
27.4 
33.7 
38.4 
33.8 
40.2 
47.8 
51.2 


24.0 
36.0 
48.4 
60.0 
68.9 
78.1 
82.7 
86.5 
95.1 


0,06 
-0.33 


1 
1 
1 
1 

-0 
-2 
-0 


.40 
,90 
.10 
,40 
.74 
,40 
,44 


Inversion  recovery  experiments  mere  conducted  on  phantom 
material  C,  with  TR-1250msBC  and  TE-30msec,  resulting  in 
the  stated  signal  intensities  CSID  and  background  values. 
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and  5  are  prasentBd  in  columns  two  and  three  respectively, 
in  arbitrary  units.  Uial  1  was  in  a  region  of  minimum 
receiver  coil  sensitivity,  while  vial  5  was  in  a  region  of 
higher  sensitivity .  Column  four  contains  the  background 
noise  intensity  values.  Two  points  of  interest  are  noted. 
One,  the  magnitude  of  all  signal  intensities  of  vial  5  are 
greater  than  the  corresponding  values  for  vial  1 .  This  is  a 
direct  result  of  the  receiver  coil  sensitivity.  Second,  the 
signal  intensities  for  the  phantom  immersed  in  the  saline 
solution  vary  dramatically  from  those  signal  intensities  for 
the  phantom  not  immersed  in  the  saline  solution.  A  possible 
explanation  is  that  the  saline  solution  absorbs  a  portion  of 
the  transmitted  rf  power,  thus  resulting  in  misset  tip 
angles  at  the  vials'  positions. 

Figures  3-5  and  3-7  present  examples  of  fitted  curves 
to  typical  data  which  were  acquired  for  Tl  determination. 
Figure  3-Ba  illustrates  the  fitting  of  ten  inversion 
recovery  data  sets  by  a  three-parameter  exponential 
function.  The  estimated  Tl  value  was  132±BmsBC,  which 
happens  to  be  within  155i  of  the  actual  Tl  value  of  ISHmsec. 
In  the  subsequent  analysis  of  this  data,  only  groups  of 
three  data  sets  and  groups  of  four  data  sets  will  be  used  to 
determine  Tl .  Each  Tl  value  determined  in  this  manner  will 
be  compared  to  both  the  natural  Tl  value  of  IS^msec,  and 
also  the  estimated  Tl  value  of  13EmsBC. 

Figure  3-Bb  illustrates  the  fitting  of  only  three  out 
of  the  possible  ten  data  sets  to  a  three-parameter 
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T|=  132±8ms 


Tj  =   131ms 


Figure  3-6.  Examples  of  fitted  data  sets  yielding  goad  Tl 
estimates.  Curve  CaD  uas  generated  by  fitting 
all  points  while  Cb)  was  generated  by  fitting 
only  the  points  denoted  by  "D". 
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T  "  61ms 


T|  -  70ms 


X        X 


T^  -   101ms 


Figure  3-7.  Examples  of  fitted  data  sets  which  yielded  poor 
Tl  estimates.  Curves  Ca),  Cb),  and  Cc)  were 
fitted  only  to  the  data  points  denoted  as  "□". 
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exponential  function.  The  three  data  points  used  for  fitting 
are  denoted  as  "D",  while  the  unused  data  points  are 
presented  as  "X".  For  this  particular  choice  of  three  data 
points,  the  estimated  Tl  value  is  131msec.  Thus,  this 
particular  value  is  both  within  15/J  of  the  all  points 
estimated  Tl  value  of  13Smssc,  and  also  within  155;  of  the 
natural  Tl  value  of  154msec.  Figure  3-7  presents  three  other 
possible  choices  of  three  data  sets,  fitted  by  a 
three-parameter  exponential  function.  In  each  case,  the 
estimated  Tl  value  is  neither  within  155i  of  the  all  points 
estimated  Tl  value  nor  the  natural  Tl  value.  These 
particular  choices  of  three  data  sets  would  be  deemed  poor. 
The  different  grouping  of  the  data  points  in  figure  3-7a,  b, 
and  c,  is  an  artifact  of  the  curve  fitting  program,  done  so 
that  the  full  fitted  curve  might  be  displayed. 

The  good  choice  of  three  data  sets  of  figure  3-6b,  and 
the  bad  choices  of  figure  3-7  illustrate  certain  basic 
characteristics.  The  good  choice  of  figure  3-Bb  depicts 
three  Tl  values  which  span  the  entire  relaxation  curve.  One 
point  defines  the  signal  intensity  at,  essentially,  time 
equal  zero.  Another  Tl  value  results  in  a  signal  intensity 
corresponding  to  almost  the  completly  relaxed  state.  Finally 
the  third  Tl  value  comes  Just  about  at  the  time  the 
relaxation  curve  changes  the  most.  Of  course,  these  are  not 
the  only  choices  of  Tl  values  which  resulted  in  good  fits, 
indeed,  there  were  many.  One  of  the  purposes  of  this 
investigation  was  to  identify  and  quantify  how  much 


variation  could  be  toloratod  from  the  near  ideal 
distribution  of  TI  values  as  depicted  in  figure  3-Bb. 
Obviously,  the  variations  illustrated  in  figure  3-7  could 
not  be  tolerated.  In  each  case  the  three  TI  values  mere 
grouped  very  near  each  other,  and  therefore  were  unable  to 
characterize  the  entire  relaxation  curve,  rather,  they  only 
characterized  the  small  region  of  the  curve  where  they  were 
located . 
Data  Analysis  flethodology 

To  facilitate  analysis  of  the  numerical  data,  a  program 
was  developed  on  an  APPLE  S-plus  microcomputer,  which 
down-loaded  the  processed  data  from  the  IBM  3330  disk  pack, 
and  permitted  graphical  display  of  correlated  parameters.  An 
example  of  one  such  output  is  given  in  figure  3-B.  The 
graphical  presentation  of  correlated  data  allowed  for  rapid 
qualitative  data  analysis.  By  this  method,  conclusions 
concerning  positional  dependence,  dependence  on  the  actual 
TI  value,  optimal  number  of  images  required,  variation  in 
accuracy  and/or  precision  in  estimated  TI  value  as  a 
function  of  phantom  immersion  in  either  air  or  saline 
solution,  choice  of  fitting  function,  and  relative  merit  of 
each  pulse  sequence  were  efficiently  and  accurately 
determined . 

The  criterion  used  in  all  evaluations  was  as  follows: 
for  a  given  set  of  parameters  Ce.g.  FIR  experiment,  actual 
TI  of  S77.8  milliseconds,  vial  number  three,  phantom 
immersed  in  saline  solution,  using  three  images,  and  fitted 
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Figure    3-B .    Sample    computer    output    of    data    correlation 
program . 
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by  a  function  with  three  degrees  of  freedom),  the  criterion 
kias  that  the  estimated  Tl  value  fitted  with  less  than  a  155i 
relative  error  to  either  the  actual  Tl  value  or  the  Tl 
estimate  obtained  by  fitting  all  the  data  pairs.  In 
comparing  different  sets  of  parameters,  what  was  actually 
compared  was  the  percentage  of  all  possible  permuted  data 
pairs  which  met  this  criterion.  This  percentage  was  denoted 
as  "range  in  Ji"  in  figure  3-B. 

The  other  features  of  figure  3-8  are  as  follow.  The 
parameter  "RUN"  refers  to  which  phantom  material  was  used  in 
the  particular  run  illustrated.  The  phrase  "RUN-KEY"  implies 
that  the  plot  is  an  overlay  plot  of  many  runs,  and  that  one 
is  refered  to  the  "Key"  which  indicates  which  platting 
character  corresponds  to  which  phantom  material.  The  "Key" 
indicates  that  the  results  of  five  different  phantom 
materials  are  presented  on  this  single  plot.  "RUN"  1  through 
5  corresponds  to  phantom  material  B  through  F. 

The  parameter  "UIAL"  indicates  from  which  of  the  eight 
passible  vials  the  data  came  from.  Similarly,  "PSEQ" 
indicates  which  pulse  sequence  was  used  for  the  Tl  estimate 
Ci.e.  PS,  SE,  IR,  FIR).  "#PAR"  refers  to  the  number  of 
parameters  in  the  fitting  function  used  for  Tl 
determination,  either  two  or  three.  "ttPTS"  refers  to  the 
number  of  data  paints  fitted  to  determine  Tl .  "ttPTS"  was 
either  two,  three,  or  four  for  "#PAR"  equal  to  two,  or 
"»PTS"  was  three  or  four  for  "#PAR"  equal  to  three.  The 
parameter  "SALT"  indicates  if  the  phantom  was  immersed  in 
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the  saline  solution  C "SALT-YES") ,  or  if  the  phantom  was  not 
immBrsed  in  the  saline  solution  ("SALT-NO").  The  parameter 
"TCOfl"  indicates  what  the  estimated  Tl  was  compared  to, 
either  the  all  points  estimated  Tl  C "TCDM-ALL") ,  or  the 
natural  Tl  C "TCOn-GLD") .  Finally,  "LEU<X"  indicates  that  the 
acceptance  levels  are  given  on  the  x  axis. 

Figure  3-9  is  an  aid  which  illustrates  how  to  interpret 
the  graphical  displays  of  correlated  parameters.  Figure  3-Sa 
illustrates  a  favorable  situation.  This  figure  illustrates 
that  approximately  50^  of  all  experiments  resulted  in  an 
estimated  Tl  value  within  5?i  of  the  standard  value.  Indeed, 
over  955i  of  all  experiments  were  within  155i  of  the  standard. 
If  this  hypothetical  curve  corresponded  to  an  inversion 
recovery  experiment,  then  it  could  be  interpreted  as 
follows.  Although  many  different  sets  of  Tl  values  were 
considered,  the  Tl  value  which  was  estimated  appeared  to 
remain  relatively  constant.  Thus,  it  would  not  be  very 
crucial  in  practice  to  optimize  the  choice  of  the  Tl  values 
used  to  gather  data  for  Tl  determination,  for  so  many 
different  combinations  were  equally  able  to  generate  an 
accurate  and  precise  estimate. 

Figure  3-9b  illustrates  a  poor  situation.  In  this 
figure,  less  than  S'<   of  all  experiments  would  result  in  an 
estimated  Tl  value  which  was  within  554  of  the  standard. 
Indeed,  it  appears  that  less  than  305;  of  all  experiments 
would  result  in  estimated  Tl  values  within  25/S  of  the 
standard.  If  this  hypothetical  curve  corresponded  to  real 
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Figure  3-9.  Hypothetical  output  of  data  correlation  program. 
Output  CaD  depicts  a  nearly  ideal  output,  while 
Cb)  depicts  a  poor  output. 
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data  it  iiiould  indicate  that  Just  a  vary  few  choices  of  the 
timing  parameters  utould  result  in  an  acceptable  Tl  estimate. 
Typically,  most  of  the  actual  curves  fell  betiueen  the  tiuo 
hypothetical  curves  of  figure  3-9. 

Typically,  natural  biological  variation  will  far  exceed 
any  machine  error,  hence  total  measurement  errors  within  105i 
for  in  vitro  experiments  are  not  uncommon  CBbB4D . 
Indeed,  total  measurement  errors  could  exceed  105s  for  in 
vivo  experiments,  where  the  investigator  has  less  control 
over  certain  biological  variables.  Although  certain 
combinations  of  parameters  resulted  in  estimated  Tl  values 
of  high  precision  Cthe  estimated  value  had  much  less  than  a 
155i  relative  error),  acceptance  at  the  IB^i  relative  error 
level  was  chosen  to  coincide  with  typical  in  vivo 
biological  variability. 

After  the  qualitative  data  analysis  based  upon  the 
graphical  display  of  correlated  parameters  was  completed,  a 
more  detailed  quantitative  analysis  was  conducted  in  order 
to  determine  those  values  of  the  variable  timing  parameter 
which  permitted  the  most  accurate  and  precise  estimate  of  Tl 
to  be  made.  A  description  of  the  analysis  method  follows. 

First,  for  a  given  set  of  experimental  parameters  Ci,e. 
pulse  sequence,  phantom  solution  concentration,  immersion  in 
air  or  saline,  particular  vial,  number  of  images  used  for 
fit,  and  the  form  of  the  fitting  function),  the  values  of 
the  variable  timing  parameter  which  resulted  in  a  fit  with 
less  than  15^  relative  error  were  noted.  By  way  of 
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illustration,  if  one  of  the  fixed  parameters  was  the  use  of 
three  images,  then  the  analysis  culminated  in  a  set  of,  for 
example,  one  hundred  groups  of  three  numbers.  Each  group  of 
three  numbers  tuas  actually  three  values  of  the  variable 
timing  parameter  CTR  for  PS  and  SE  images,  TI  for  IR  and  FIR 
images^  which  resulted  in  a  good  fit.  Each  value  was  then 
scaled  to  the  Tl  value  which  the  fit  was  being  compared 
to.  Next,  a  linear-multiple-regression  analysis  CSpBlD  was 
performed  on  this  set  of  groups,  yielding  a  regression 
equation  which  related  the  three  values  of  the  variable 
timing  parameters  to  each  other.  This  is  the 
multidimensional  analog  to  the  least-squares  fit  line  used 
for  two  dimensional  data  sets.  Hence,  for  three  images,  a 
least  squares  fit  plane  was  determined  which  incoporated  the 
empirical  data  into  an  analytical  expression.  This 
analytical  expression  could  be  used  a  priori  in  selecting 
values  of  the  variable  timing  parameter,  or  a  posteriori 
in  evaluating  the  group  of  timing  parameters  used  in  a 
series  of  images. 

Results 
There  were  five  basic  results  of  this  investigation, 
two  of  which  are  of  primary  importance.  Some  of  the  results 
were  of  a  general  nature,  such  as  the  dependence  of 
measurement  upon  the  position  within  the  field  of  view,  and 
the  variation  in  the  accuracy  and  precision  in  the  estimated 
Tl  value  as  a  function  of  the  actual  Tl  value.  Other  results 
were  more  specific,  such  as  the  determination  of  how  many 
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images  acquired  uith  which  pulse  sequence,  and  fitted  by 
which  function,  resulted  in  accurate  and  precise  estimated 
Tl  values,  independent  of  the  presence  of  the  rf  attenuating 
saline  solution.  Finally,  particular  results  uiere 
quantitative,  for  example,  the  linear  multiple  regression 
analysis  of  the  values  of  the  variable  timing  parameter 
which  resulted  in  accurate  and  precise  11  estimates. 

The  performance  of  the  PS  experiment  was  so  poor  in 
contrast  to  the  SE,  IR,  and  FIR  experiments,  for  the  reasons 
offered  in  the  discussion  section  of  this  chapter,  that  it 
was  discounted  as  a  viable  method  of  Tl  determination.  The 
fallowing  results  apply  only  to  SE,  IR,  and  FIR  experiments. 
Representative  data  are  presented  in  support  of  all  results. 
The  reproducibility  of  signal  intensities  for  a  given  set  of 
parameters  was  at  all  times  greater  than  357i,  and  typically 
greater  than  97.55s.  The  reproducibility  was  determined  from 
a  series  of  measurements  which  were  all  repeated  ten  times. 
Positional  Dependence 

All  estimated  Tl  values  were  constant  within  155i  to  the 
position  within  the  field  of  view  from  which  the  individual 
signal  intensities  were  recorded.  That  is,  although  spatial 
variations  in  signal  intensity  occurred,  calculated  values 
of  Tl  did  not  exhibit  these  variations,  within  acceptable 
experimental  limits.  The  signal  intensity  variations 
resulted  from  the  intrinsic  magnetic  field  inhomogeneity , the 
intrinsic  transmitted  rf  field  inhomogeneity,  and  mostly, 
the  spatially  inhomogeneous  rf  receiver  coil  response. 
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The  representativB  data  of  figure  3-10  support  this 
result.  Uial  three  was  positioned  at  the  site  of  maximum 
receiver  sensitivity,  uhile  vial  five  utas  located  at  a 
position  of  poor  sensitivity .  A  comparable  percentage  of 
experiments  met  the  acceptance  criterion  (fitted  Tl  value 
had  less  than  a  155;  relative  error 5  for  both  vials.  Although 
similar  estimated  Tl  values  luere  calculated  for  both  vials, 
typically  the  standard  deviations  in  the  fits  for  vial  five 
were  larger  than  those  for  vial  three,  owing  to  the  smaller 
signal  to  noise  ratio  of  the  intensities  at  that  position. 
Dependence  Upon  Actual  Tl  Ualue 

The  percentage  of  experiments  which  met  the  acceptance 
criterion  was  constant  within  155i  to  the  actual  Tl.  That  is, 
a  given  set  of  parameters  Cwith  the  values  of  the  variable 
timing  parameter  suitably  chosen  for  the  actual  Tl) 
generated  a  similar  number  of  good  fitting  Tl  estimates, 
independent  of  the  actual  Tl .  Also,  for  a  given  set  of 
parameters  characterizing  experiments  which  met  the 
acceptance  criterion,  the  values  of  the  variable  timing 
parameter  normalized  to  the  actual  Tl  were  constant  within 
155i  to  the  actual  Tl .  Thus,  it  was  possible  to  characterize 
the  optimal  values  of  the  timing  parameter  independent  of 
Tl .  The  representative  data  of  figure  3-11  support  this 
result.  Each  different  "RUN"  represents  a  different  actual 
Tl  value  of  the  phantom  material .  The  close  grouping  of  the 
data  along  the  "RANGE  IN  '<"    axis,  as  a  function  of 
"ACCEPTANCE  LEUEL",  indicates  the  insensitivity  to  Tl  value. 
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Figure  3-10.  Representative  data  indicating  independence  to 
position  in  f ield-of-view .  Ca)  Uial  5  output, 
Cb)  identical  output  for  vial  3. 
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Figure  3-11.  RepresBntativB  data  indicating  independancB  to 
actual  Tl  value. 
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DBpendence  Upon  Number  oF  Imacjes 

For  a  given  set  of  parameters,  the  percentage  of 
experiments  which  met  the  acceptance  criterion  utilizing  a 
two-parameter  fitting  function  was  constant  within  15/i  to 
the  acquisition  and  use  of  two,  three,  or  four  images  for  Tl 
estimation.  Similarly,  the  percentage  of  experiments  which 
met  the  acceptance  criterion  utilizing  a  three-parameter 
fitting  function  was  constant  within  155j  to  the  acquisition 
and  use  of  three  or  four  images  for  Tl  estimation.  That  is, 
a  given  set  of  parameters,  with  the  degree  of  the  fitting 
function  constant,  generated  a  similar  number  of  good 
fitting  Tl  estimates,  independent  of  the  number  of  processed 
images . 

There  was  a  corollary  result.  As  explained  previously, 
IR  experiments  were  fitted  twice,  once  by  a  function  with 
two  degrees  of  freedom  C given  by  setting  C  in  equation  3-E 
equal  to  2K),  and  once  by  a  function  with  three  degrees  of 
freedom  Cgiven  by  equation  3-2) .  For  each  form  of  the 
function,  all  IR  experiments  were  constant  within  15/i  to  the 
number  of  images  used,  although  fitting  with  the 
three-parameter  function  resulted  in  a  substantially  higher 
percentage  of  experiments  which  met  the  acceptance 
criterion . 

The  representative  data  of  figure  3-12  support  these 
results.  Additionally,  the  results  of  fitting  FIR 
experiments  with  a  function  with  two  degrees  of  freedom  are 
also  presented  in  figure  3-13. 
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Representative  data  indicating  independence  to 
number  of  images.  CaD  Using  three  images  for  Tl 
determination,  Cb)  using  four  images  for  Tl 
determination . 


51 


lOOi 


100 


OjOiCOCvJOO  cmc>jco<nco 

^    E    E    E    EE  '^    "^    E    E    E 

CO  U-      u.  w  u.      u. 

Salt  Salt 

Estimated  I,  Natural  T. 


CMCMCOCMCO  OJOJCOCMCO 

lU     cc     (T     i     cc  tU     a:     of     i     of 

No  Salt  No  Salt 

Estimated  T,  Natural  T, 


Figure   3-13.    Summary    of    results, 
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Optimal  Sets  of  Parameters 

To  obtain  a  direct  estimate  of  the  actual  Tl  value  in 
an  object  tuhich  may  be  uiithin  a  spatially  Inhomogeneous  rf 
field,  the  most  accurate,  precise,  and  time  effective 
technique  to  use  was  three  FIR  images,  with  suitably  chosen 
values  of  TI,  fitted  by  a  three-parameter  exponential 
function.  The  use  of  IR  images  luas  equally  accurate  and 
precise,  but  not  as  temporally  efficient. 

Three  SE  images,  with  suitably  chosen  values  of  TR,  and 
fitted  by  a  two-parameter  exponential  function  can  be  more 
precise  and  time  effective  than  the  FIR  technique  in 
estimating  Tl,  but  were  always  much  less  acurate,  and  prone 
to  error  if  an  inhomogeneous  rf  field  was  present.  The 
representative  data  of  figure  3-13  support  these  results. 
Optimum  Uariable  Timing  Parameter  Ualues 

For  each  of  the  three  above  mentioned  optimum  sets  of 
parameters,  a  linear  multiple  regression  analysis  was 
conducted  on  the  empirical  values  of  the  variable  timing 
parameters  which  resulted  in  estimated  Tl  values  of  155i 
relative  error  or  less.  For  each  regression  analysis,  a 
coefficient  of  linear  multiple  correlation  was  determined. 
The  coefficient  may  lie  between  0  and  1.  The  closer  it  was 
to  1,  the  better  was  the  linear  relationship  between  the 
variables,  with  a  value  of  1  indicating  a  perfect 
correlation.  The  closer  it  was  to  0,  the  worse  was  the 
linear  relationship.  For  all  cases,  the  coefficient  of 
linear  multiple  correlation  was  0.5  or  greater. 
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The  use  of  three  FIR  images,  fitted  by  a  three 
parameter  exponential  function,  was  the  optimal  method  of 
estimating  the  actual  Tl  value,  in  the  presence  of  a 
spatiallg  inhomogeneous  rf  field.  The  relationship  between 
the  three  values  of  TI  in  these  FIR  experiments  was 

H  -  1.1  +  O.Bn  -  O.IL  C3-3] 

where  the  three  TI  values,  H,  n,  and  L  were  the  highest, 
middle,  and  lowest  values  all  scaled  to  Tl  by  dividing  the 
specific  timing  parameter  by  the  actual  Tl  value.  The 
standard  error  of  estimate  in  H  was  0.3,  for  fl  it  was  0.3, 
and  for  L  it  was  0.2. 

The  use  of  three  IR  images,  fitted  by  a  three-parameter 
exponential  function,  was  another  method  of  estimating  the 
actual  Tl  value,  in  the  presence  of  a  spatially 
inhomogeneous  rf  field.  The  relationship  between  the  three 
values  of  TI  in  the  IR  experiments  was 

H  -  l.B  +  l.in  +  O.BL  C3-4D 

where  the  three  TI  values,  H,  fl,  and  L  were  the  highest, 
middle,  and  lowest  values  normalized  to  the  actual  Tl  value. 
The  standard  error  of  estimate  in  H  was  1.2,  for  fl  it  was 
0.4,  and  for  L  it  was  0.3. 

A  large  percentage  of  SE  experiments  did  not  meet  the 
acceptance  criterion  when  the  estimated  Tl  value  was 
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comparBd  to  the  actual  Tl  value.  Houiever,  there  uias  a  large 
percentage  of  SE  experiments  mhich  did  meet  the  acceptance 
criterion  when  the  calculated  Tl  value  uias  compared  to  the 
value  of  Tl  estimated  by  the  fitting  of  all  data  pairs.  This 
indicated  that  uhile  estimated  Tl  values  from  SE  images  were 
not  accurate,  they  were  precise.   The  relationship  betueen 
the  three  values  of  TR  for  this  set  of  parameters  in  the  SE 
experiment  uias 

H  -  1.7  +  l.in  -  0.3L  C3-5D 

where  the  three  values,  H,  H,  and  L  tuere  the  highest, 
middle,  and  loiuest  values  normalized,  in  this  case,  to  the 
Tl  value  estimated  by  fitting  all  data  pairs.  The  standard 
error  of  estimate  in  H  was  1.1,  for  M  it  was  0.7,  and  for  L 
it  was  0.3. 

Discussion 

The  in  vivo  determination  of  Tl  values  is  complicated 
by   many  factors,  some  of  which  were  addressed  by  this 
investigation. The  specific  aim  of  this  investigation  was  to 
evaluate  empirically  the  suggestions  offered  by  the  NfIR 
spectroscopy  literature  on  Tl  determination  and  as  to  how 
well  they  applied  to  NfIR  imaging,  with  its  unique  set  of 
experimental  constraints. 

The  NfIR  spectroscopy  literature  offered  suggestions  on 
which  pulse  sequence  to  use,  the  number  and  values  of  the 
variable  timing  parameter  to  use,  and  the  form  of  the 
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Fitting  function  required  to  calculate  Tl  from  the 
individual  signal  intensity  measurements.  The  added 
constraints  imposed  by  NMR  imaging  were  the  follouiing: 
measurements  could  be  obtained  from  any  position  uiithin  a 
large  field  of  view,  typical  Tl  values  in  the  human  body 
span  a  range  of  over  an  order  of  magnitude,  total  imaging 
times  had  to  be  minimized,  and  the  human  body  has  dielectric 
properties  uihich  resulted  in  spatially  dependent  rf 
attenuation . 
Positional  Dependence 

It  was  determined  that  all  estimated  Tl  values  were 
constant  within  15Ji  to  the  position  within  the  field  of  view 
from  which  the  individual  signal  intensities  were  recorded. 
If  there  were  no  spatial  variation  in  signal  intensity  Cfor 
fixed  experimental  parameters),  then  this  would  be  an 
expected  result.  As  demonstrated  by  the  signal  intensity 
variations  in  figure  3-5,  this  was  not  the  case. 

The  predominant  cause  of  spatial  variations  in  signal 
intensity  measurements  was  the  spatially  inhomogeneous  rf 
receiver  coil  response.  The  half  saddle  shape  of  the  rf 
receiver  coil  produced  an  axially  asymmetric  spatial 
response,  as  seen  in  figure  3-5.  Thus,  for  a  hypothetical 
sample  that  filled  the  field  of  view  and  produced  a 
constant,  homogeneous  signal  from  each  position,  the 
recorded  signal  intensity  would  be  the  constant  intensity 
convoluted  with  the  rf  coil's  sensitivity  response  at  that 
location.  This  would  influence  only  K  and  C  of  equation  3-2, 
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and  not  the  exponential  time  constant.  Thus,  this 
investigation  verified  experimentally  the  theoretical 
prediction . 

This  result  is  significant  for  the  following  reason. 
Although  it  is  tempting  to  make  a  differential  diagnosis 
based  upon  differences  in  signal  intensities  Ce.g.  if  the 
liver  is  more  intense  than  the  spleen  then  diagnosis  A,  if 
vice  versa  then  diagnosis  B),  these  differences  may  not  be 
entirely  organic  in  nature.  Uariations  in  calculated  Tl 
values  Cfor  a  constant  set  of  imaging  parameters)  are  a  more 
reliable  indicator  of  true  clinical  variation. 
Dependence  Upon  Actual  Tl  Ualue 

It  uas  determined  that  the  percentage  of  experiments 
uhich  met  the  acceptance  criterion  uas  constant  within  155i 
to  the  actual  Tl  value,  for  properly  chosen  values  of  the 
variable  timing  parameter.  Since  experiments  were  conducted 
over  an  actual  Tl  range  of  approximately  100  milliseconds  to 
1000  milliseconds,  this  result  is  valid  only  over  this 
range.  To  justify  this  empirical  result,  it  is  necessary  to 
consider  two  factors,  the  physical  model  of  spin-lattice 
relaxation,  and  any  instrumentational  dependence  upon  Tl 
determination . 

Spin-lattice  relaxation  theory  predicts,  for  such  a 
simple  phantom  material  Cparamagnetically  doped  water),  a 
monoexponential  relationship  between  the  signal  intensity 
and  the  variable  timing  parameter.  The  spin-lattice 
relaxation  rate  is  directly  obtainable  from  the  exponential 
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time  constant.  This  relatively  simple  relationship  holds  no 
Tl  dependent  bias.  That  is,  all  other  factors  being  equal, 
Tl  is  simply  a  scaling  Factor  in  the  exponential  argument, 
and  does  not  influence  the  general  form  of  the  function. 
Thus,  the  empirical  invar iance  to  Tl  was  to  be  expected,  on 
the  basis  of  the  physical  model . 

Although  instrumental  effects  on  Tl  determination  can 
have  many  causes,  ranging  from  rf  field  production  to 
computer  roundoff  errors,  there  could  not  be  any  direct 
instrumentational  dependence  upon  Tl ,  for  obviously  the 
instrument  could  have  no  knoiuledge  of  the  phantom's  Tl 
value.  Thus,  any  Tl  dependent  bias  uould  have  to  have  been 
indirectly  related.  Since  Tl  calculations  are  based  upon 
signal  intensities  and  the  values  of  a  variable  timing 
parameter,  any  direct  instrumental  bias  toiuards  these 
parameters  mill  indirectly  affect  Tl  determinations. 

The  accuracy  in  defining  a  timing  parameter  bias 
governed  by  the  computer's  CPU  clock.  If  me  consider  a  clock 
frequency  of  one  megahertz  Ca  value  much  less  than  even 
modern  personal  computers,  let  alone  the  PDP  11/24),  then 
timing  events  could  be  controlled  to  within  a  microsecond. 
Since  typical  values  of  timing  parameters  were  of  the  order 
of  tens  or  hundreds  of  milliseconds,  it  seems  unlikely  that 
incorrectly  set  timing  parameters  were  a  large  source  of 
error  in  the  determination  of  Tl  values. 

On  the  other  hand,  the  signal  intensity  was  influenced 
by  a  myriad  of  variables,  most  of  which  were  instrumental  in 
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nature.  For  example,  all  experiments  in  this  investigation 
utilized  a  single  slice  mode  of  acquisition.  The  use  of 
multislice  acquisition  introduced  other  variables  which 
influenced  signal  intensities,  if  all  other  factors  remained 
constant . 

Perhaps  the  greatest  instrumental  influence  upon  signal 
intensity  uas  the  required  use  of  spin  echo  formation  and 
imaging  magnetic  gradients  for  image  signal  acquisition. 
Conventional  Fourier  imaging  CEdBOD  relies  upon  the 
detection  of  spin  echoes.  This  did  not  introduce  further 
complications  in  the  SE  images,  but  PS,  IR,  and  FIR 
experiments  conventionally  generate  a  free  induction  decay 
CFID)  signal,  and  not  an  echo.  Indeed,  this  complication 
could  explain  the  poor  performance  of  the  PS  sequence.  The 
PS  experiment  relied  upon  the  net  magnetization's  reaching  a 
steady-state  value  while  being  subjected  to  repetitive  w/E 
rf  pulses,  but  by  its  very  nature  a  spin  echo  had  a 
dynamically  varying  net  magnetization.  Although  the  IR  and 
FIR  experiments  do  not  force  the  net  magnetization  to  a 
steady-state  value,  altering  their  pulse  sequence  to  include 
echo  formation  could  have  been  a  source  of  signal  intensity 
error.  The  formation  of  spin  echoes  required  the 
magnetization  to  remain  in  the  transverse  plane  Cin  the 
rotating  coordinate  system)  for  an  appreciable  period  of 
time.  This  introduced  TE  damping  of  the  signal  intensity, 
and  in  the  presence  of  magnetic  field  gradients,  damping  due 
to  molecular  self-diffusion  CStBSD . 
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The  empirically  determined  invariance  to  the  actual  Tl 
value  is  significance  for  two  reasons.  First,  it  permitted 
the  development  of  gfeneral  recommendations  concerning  the 
optimal  values  of  the  variable  timing  parameter.  Second,  it 
indicated  that  concerns  raised  about  instrumental  effects  on 
Tl  determination  did  not  introduce  any  systematic  errors  in 
measurement,  possibly  only  statistical  errors.  Thus,  the 
results  of  this  investigation  could  be  applied  to  in  vivo 
Tl  determination  without  loss  of  accuracy  due  to  any  Tl 
dependence,  to  the  extent  that  the  phantom  used  in  this 
investigation  modeled  the  in  vivo  object. 
Dependence  Upon  Number  of  Images 

It  was  determined  that  for  a  given  set  of  parameters, 
and  with  a  fitting  function  with  two  degrees  of  freedom,  the 
percentage  of  experiments  which  met  the  acceptance  criterion 
was  constant  within  155i  for  the  acquisition  and  use  of  two, 
three,  or  four  images  for  Tl  estimation.  Similarly,  for  a 
fitting  function  with  three  degrees  of  freedom,  the 
percentage  of  experiments  which  met  the  acceptance  criterion 
was  constant  within  lS>i   to  the  acquisition  and  use  of  three 
or  four  images  for  Tl  estimation.  This  result  agrees  with 
the  NtlR  spectroscopy  literature,  which  states  that  to  a 
first  approximation,  the  error  in  estimated  Tl  values  does 
not  depend  upon  the  number  of  data  pairs  used  in  the 
calculation,  provided  the  rms  error  of  the  experimental  data 
was  less  than  one  tenth  of  the  signal  intensity  at  time 
equal  to  infinity. 
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In  this  investigation,  an  sstimated  Tl  value  utas  deemed 
a  goad  approximation  to  the  actual  Tl  value  if  the  relative 
error  uias  less  than  15/S.  For  this  acceptance  level,  only 
First  order  eFFects  ujere  oF  suFFicient  magnitude  to 
influence  the  error  in  the  estimated  Tl  value.  Therefore  the 
observed  invariance  to  the  number  of  data  pairs  used  in 
calculating  the  spin-lattice  relaxation  time  corresponds  to 
the  NflR  spectroscopy  literature's  suggestion  concerning  this 
point  CBeBOD . 

This  result  uas  of  particular  significant  for  the 
following  reason.  In  NtIR  imaging,  time  is  at  a  premium.  This 
investigation  sought  to  identify  the  method  of  spin-lattice 
relaxation  time  determination  which  maximized  accuracy  and 
precisian,  and  minimized  the  total  imaging  time.  One  way  to 
minimize  the  total  imaging  time  was  to  take  the  least  number 
of  images  required  to  generate  mathematically  unique  and 
statistically  significant  results.  Thus,  for  those 
techniques  which  were  fitted  by  a  function  with  two  degrees 
of  Freedom,  two  images  sufficed,  while  three  were  required 
for  fitting  functions  with  three  degrees  of  freedom. 

Since  the  earliest  beginning  of  NMR  imaging,  some  have 
made  the  suggestion  that  the  best  method  of  in  vivo  Tl 
determination  was  the  fitting  of  two  SE  images  with  a 
function  with  two  degrees  of  freedom  LHrB3,  na79D .  This 
result  indicates  that,  if  the  signal  to  noise  ratio  was 
adequate  in  each  image,  two  images  would  indeed  be  the 
resonable  number  to  acquire  and  process.  This  does  not  imply 
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that  this  methad  was  optimal,  but  simply  that  it  conFormBd 
to  this  "least  number  of  images"  result.  Indeed,  results  of 
this  work  indicated  that  although  the  two  image  SE  method 
was  optimal  in  reducing  the  total  imaging  time,  it  was  far 
from  optimal  in  terms  of  its  accuracy  in  estimating  Tl . 

The  summary  of  results  presented  in  figure  3-13 
contains  much  information  concerning  the  appropriate  choice 
of  pulse  sequence/fitting  function  for  a  given  situation. 
The  precision  of  each  sequence  was  indicated  by  the 
corresponding  percent  of  experiments  meeting  the  acceptance 
criteria  when  the  calculated  Tl  was  compared  to  the 
"Estimated  Tl"  value.  The  least  demanding  situation  was  when 
the  phantom  was  not  immersed  in  the  saline  solution.  For 
this  situation  all  combinations  of  pulse  sequence/fitting 
function  were  essential  equal  in  precision.  This  was  not 
true  when  the  phantom  was  immersed  in  the  saline  solution. 
In  this  situation  the  SE:B  method  was  slightly  more  precise 
than  any  other  methad  of  Tl  determination,  all  of  which 
shared  an  essentially  equal  precision. 

Those  experiments  whose  calculated  Tl  values  met  the 
acceptance  criteria  when  compared  to  the  "Natural  Tl" 
represented  methods  of  determining  precise  and  accurate  Tl 
estimates.  For  this  situation  both  the  IR:2  and  FIR:2 
methods  were  inadequate,  whether  the  phantom  was  immersed  in 
the  saline  solution  or  not.  Additionally,  both  the  IR:3  and 
FIR:3  methods  of  Tl  determination  were  best,  while  the  SE:2 
method  performed  marginally. 
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It  was  of  importance  to  document  to  what  degree  a 
particular  pulse  sequence/fitting  function  varied  from  the 
"Salt"  to  "No  Salt"  case.  To  facilitate  that  analysis  the 
percent  relative  difference  was  calculated  from  figure  3-13 
for  both  the  "Estimated  Tl"  and  "Natural  Tl"  results.  The 
percent  relative  difference  was  determined  by  subtracting 
the  "^  fleeting  Criteria"  of  the  "Salt"  case  from  the  "Ss 
deeting  Criteria"  of  the  "No  Salt"  case,  and  then  the  result 
was  divided  by  the  "\   Heeting  Criteria"  of  the  "No  Salt" 
case.  This  was  dons  for  each  pulse  sequence/fitting 
function,  for  both  the  "Estimated  Tl"  and  the  "Natural  Tl" 
results,  with  the  outcomes  given  in  figure  3-14. 

Ideally,  a  given  pulse  sequence/fitting  function  would 
perform  as  well  whether  the  phantom  was  immersed  in  the 
saline  solution  or  was  not  immersed  in  the  saline  solution. 
Since  this  is  the  ideal  situation,  the  ideal  percent 
relative  difference  would  be  zero.  The  SE:E  method  of  Tl 
determination  suffered  the  smallest  loss  of  precision,  as 
depicted  in  the  "Estimated  Tl"  result  in  figure  3-14.  The 
"Natural  Tl"  results  indicated  that  the  IR:3  was  least  prone 
to  lasses  of  accuracy  and  precision,  while  the  SE:B  and 
FIR: 3  methods  were  nearly  as  lossless.  These  results  must 
not  be  considered  out-of-context .  For  example,  the  result 
that  the  SE:2  method  had  a  low  percent  relative  difference 
for  the  "Natural  Tl"  case  viewed  in  conjunction  with  the 
results  of  figure  3-13  indicated  that,  in  this  particular 
situation,  the  SE:E  method  went  from  being  a  poor  method  of 
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Figure  3-14.  Effect  of  rf  attenuating  material  upon 
mBasurement  precisian  and  accuracy . 
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dBtermining  accurate  and  precise  Tl  values  when  the  phantom 

uas  not  immersed  in  a  saline  solution,  to  being  a  slightly 

poorer  method  when  the  phantom  uias  immersed  in  the  saline 

solution! 

ft  Priori  Recommendationa 

All  of  the  results  of  this  investigation  dealt  with  the 
determination  of  the  optimal  method  of  in  vivo  Tl  values. 
Of  primary  inportance  were  the  results  concerned  luith  the 
optimal  pulse  sequence  and  optimal  values  of  the  variable 
timing  parameter.  These  particular  results  may  be  used  a 
priori  to  design  an  imaging  scheme  which  results  in 
reliable  Tl  values,  or  a  posteriori  to  evaluate  the 
reliability  of  a  calculated  Tl  value.  The  a  priori 
recommendations  folloiu. 
Optimal  sets  of  parameters 

The  first  decision  to  be  made  in  the  design  of  a  Tl 
imaging  series  is  the  choice  of  pulse  sequence.  This 
decision  is  based  upon  three  factors,  total  imaging  time, 
and  the  accuracy  and  precision  in  Tl  determination.  The 
results  indicate  that  there  mere  three  choices,  three  FIR 
images  fitted  by  a  three-parameter  exponential  function, 
three  IR  images  fitted  by  a  three-parameter  exponential 
function,  or  three  SE  images  fitted  by  a  tujo-parameter 
exponential  function. 

The  FIR  series  of  Tl  images  had  a  number  of  advantages 
over  the  IR  or  SE  series,  which  made  it  the  recommended 
method.  First,  the  method  required,  at  most,  half  the  total 
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time  of  the  corresponding  IR  series.  Second,  it  yielded 
consistent  results,  independent  of  whether  the  saline 
solution  was  present  or  not.  Finally,  it  generated  estimates 
of  Tl  which  were  both  accurate  and  precise.  This  meant  that 
in  vivo  Tl  values  could  be  obtained  directly,  without  the 
need  of  a  calibrated  set  of  measurements. 

There  were  also  some  disadvantages  to  the  FIR  series  of 
Tl  images.  First,  being  an  inversion  type  of  experiment, 
phase  reconstruction  of  the  images  was  required  to  make  use 
of  the  full  dynamic  range  offered  by  the  technique.  In 
practice,  phase  reconstruction  is  not  always  easily 
accomplished.  Second,  it  did  not  always  offer  the  minimum 
total  imaging  time.  Often,  the  SE  series  of  images  would 
result  in  less  total  imaging  time.  Finally,  the  SE  series 
often  produced  more  precise  measurements  of  Tl  than  the  FIR 
series  did.  Although  the  FIR  series  had  some  flaws,  overall 
it  was  the  most  rugged  method,  yielding  accurate  and  precise 
estimates  of  Tl  under  the  most  adverse  of  situations. 

The  IR  series  shared  many  of  the  attributes  of  the 
recommended  FIR  series,  but  had  an  intrinsic  tragic  flaw.  An 
IR  image  necessitated  the  use  of  a  long  total  image  time. 
This  long  TR  would  often  clash  with  the  clinical  necessity 
of  speed.  An  advantage  the  IR  series  had  over  the  FIR  was 
that,  due  to  the  extended  TR  time,  less  saturation  of  signal 
occurred.  Thus,  the  IR  images  had  slightly  more  contrast,  an 
aid  for  visual  evaluation,  but  no  real  advantage  in  the 
quantitative  analysis  required  for  Tl  determination. 
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The  enigma  of  tha  investigation  was  the  SE  method  of 
spin-lattice  relaxation  time  determination.  The  SE  series 
had  many  advantages.  One,  it  often  resulted  in  the  shortest 
total  imaging  time.  Tuo,  only  simple  magnitude 
reconstruction  uas  required.  Finally,  the  SE  series  oftened 
attained  estimated  spin-lattice  relaxation  time  values  with 
very  high  precision.  Additionally,  although  not  a  result  of 
this  investigation,  the  SE  technique  is  often  the  clinical 
technique  of  preference,  displacing  the  inversion  techniques 
consistently.  For  these  reasons,  the  SE  method  of 
spin-lattice  relaxation  time  determination  appears  to  enjoy 
the  most  widespread  use.  Additional  results  of  this 
investigation  indicated  that  the  SE  method  had  some  serious 
shortcomings,  and  yielded  quality  results  only  in 
restrictive  situations. 

The  most  important  disadvantage  of  the  SE  method  was 
its  lack  of  accuracy  in  estimating  Tl  values.  The  SE  method 
always  fared  best  when  its  estimated  Tl  value  was  compared 
to  the  Tl  value  determined  by  fitting  all  points,  rather 
than  when  it  was  compared  to  the  actual  Tl  value.  Thus,  in 
vivo  Tl  values  could  not  reliably  be  determined  directly. 
Instead,  some  type  of  calibration  curve  would  need  to  be 
generated  to  permit  actual  Tl  values  to  be  gleaned  from  SE 
estimated  Tl  values.  This  conclusion  agrees  with  other 
reported  investigations  into  the  use  of  SE  images  for  Tl 
determination  CPyB3D . 
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Optimum  variable  timing  parameter  values 

Dnce  the  decision  as  to  uihich  pulse  sequence  to  be 
utilized  is  made,  either  the  SE,  IR,  or  FIR  imaging 
sequence,  it  is  necessary  to  specify  the  values  of  the 
variable  timing  parameter.  For  FIR  and  IR  images  this  is  TI, 
while  for  SE  images  it  is  TR .  To  insure  optimal  accuracy  and 
precision,  the  results  of  the  linear  multiple  regression 
analysis  should  be  used.  The  derived  linear  relationships 
among  the  three  values  of  the  variable  timing  parameter  for 
the  FIR,  IR,  and  SE  methods  have  been  given  in  aquations 
3-3,  3-4,  and  3-5  respectively. 

There  are  some  general  guidelines  to  be  offered  as  to 
houi  best  to  use  the  suggested  results  a  priori .  First,  the 
greater  the  uncertainty  in  the  Tl  value  under  investigation, 
the  greater  should  be  the  spread  in  the  three  values.  The 
lowest  value,  L,  could  be  picked  to  correspond  to  the 
shortest  permissible  value  of  the  particular  timing 
parameter.  For  the  instrumentation  used  in  this 
investigation,  that  luas  E5  milliseconds  for  TI  and  100 
milliseconds  for  TR .  The  highest  value  of  the  timing 
parameter,  H,  mould  vary  depending  upon  the  pulse  sequence. 
For  the  FIR  and  IR  methods,  H  could  correspond  to  some  value 
close  to  TR,  say  0.3TR.  The  solution  is  not  as  straight 
forward  for  the  SE  method,  where  there  is  no  fixed  timing 
parameter  of  relavance  on  which  to  base  the  choice.  In  this 
case,  H  should  be  chosen  with  consideration  towards  the 
largest  value  in  the  expected  Tl  range. 
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Often  it  is  desired  to  merge  the  Tl  determining  images 
with  the  images  desired  for  visual  evaluation.  In  this  case, 
one  or  tuio  of  the  values  for  H,  f1,  or  L  may  be  taken  to 
correspond  to  the  values  associated  with  the  images  required 
for  visual  evaluation,  with  the  remaining  timing  parameter 
determined  by  the  linear  multiple  regression  equation.  For 
example,  tuio  SE  images  uiith  the  TR  values  of  500 
milliseconds  and  1000  milliseconds  are  desired  for  visual 
evaluation.  Additionally,  an  in  vivo  Tl  determination  is 
desired  of  a  tissue  uihose  Tl  is  thought  to  be  about  650 
milliseconds.  By  setting  L  and  fl  to  0.77  and  1.54, 
respectively,  and  utilizing  the  linear  multiple  regression 
equation  for  the  SE  method,  H  is  determined  to  be  3.16, 
corresponding  to  a  TR  value  of  S055  milliseconds.  Note  that 
nonsense  answers  would  result,  in  this  particular  case,  if 
the  two  existing  TR  values  were  assigned  to  L  and  H,  or  n 
and  H.  Additionally,  if  the  Tl  of  concern  was  thought  to  be 
350  milliseconds,  then  assigning  n  to  be  0.77,  and  H  to  be 
1.54  yields  the  third  value  of  TR  as  145  milliseconds.  By  a 
similar  process,  the  optimal  values  of  the  timing  parameters 
in  FIR  and  IR  series  are  determined. 

In  the  manner  Just  described,  an  optimal  method  of  in 
vivo  Tl  determination  may  be  developed  a  priori,  based 
upon  the  available  imaging  time,  desired  accuracy  and 
precision,  and  range  of  Tl  under  investigation. 
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ft   Posteriori  Recommendations 

It  frequently  happens  that  the  desire  for  a 
quantitative  determination  of  Tl  is  not  expressed  until 
after  the  imaging  series  is  completed.  If  the  images  mere 
acquired  along  the  lines  of  either  the  FIR,  IR,  or  SE 
methods  outlined  earlier  Ci.e.  three  images  uere  acquired, 
all  with  the  same  number  of  signal  averages),  then  it  would 
be  possible  to  generate  an  estimated  Tl  value,  and  with  the 
aid  of  the  linear  multiple  regression  equations,  determine 
the  reliability  of  the  estimate.  The  procedure  is  explained 
in  the  following  scenario. 

Three  inversion  type  images  were  acquired,  with  TE 
equal  to  30  milliseconds,  TR  equal  to  1500  milliseconds,  and 
three  values  of  Tl  equal  to  50,  550,  and  B50  milliseconds. 
All  images  were  acquired  with  four  signal  averages.  A 
three-parameter  exponential  is  fitted  to  the  data,  and  an 
estimated  Tl  value  of  700  milliseconds  is  calculated.  The 
question  is:  is  this  a  reasonable  Tl  value  to  expect  this 
series  of  images  to  properly  characterize,  or  is  it  Just  the 
result  of  some  mathematical  fitting  routine  which  is 
insensitive  to  certain  physical  realities? 

Since  the  value  of  TR  is  approximately  twice  the  Tl 
estimate,  this  series  is  comparable  to  the  FIR  method  of  Tl 
determination.  The  linear  multiple  regression  equation  for 
the  FIR  series  indicates  that  the  optimal  Tl  value,  for  this 
set  of  Tl  values,  is  477  milliseconds.  This  is  obviously  not 
equal  to  the  estimated  value  of  700  milliseconds,  but  the 
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linear  multiple  regression  equation  Is  not  without  error, 
presented  in  the  form  of  standard  error  of  estimates  in  L, 
n,  and  H. 

By  a  simple  propagation  of  error  analysis  of  equation 
3-3,  the  standard  error  of  estimate  in  Tl  for  the  FIR  method 
of  Tl  determination,  cCTlD,  may  be  determined  to  be  152 
milliseconds.  Thus,  the  optimal  Tl  value  to  be  determined  by 
the  set  of  Tl  values  used  in  this  example  is  4771152 
milliseconds.  The  estimated  value,  700  milliseconds,  is 
nearly  1.5  standard  errors  greater  than  the  optimal  mean 
value,  and  it  is  therefore  deemed  an  unreliable  estimate. 
This  analysis  could  be  conducted  as  easily  for  any  IR  or  SE 
series  of  images. 


CHAPTER  lU 
EXPLOITING  THE  STIMULATED  ECHO  IN  NMR  IMAGINB 

Introduction 

There  is  an  adage  that  applies  equally  to  all 

multipulsB  NriR  experiments  uhich  states  it  is  easier  to 

induce  spin  echoes  than  not.  Rather  than  ignore  or  suppress 

these  additional  echoes  in  NflR  imaging  experiments  CDuSHD, 

this  investigation  sought  to  glean  added  information  From 

them.  In  particular,  this  study  exploited  the  unique 

properties  of  the  stimulated  echo  CSTE) ,  as  first  identified 

by  Hahn  CHaSOJ,  and  further  quantified  by  Uloessner 

CWoBlD.  Although  neu  to  NMR  imaging  CFrBB,  Ha85,  SaBSa, 

SaBBbD ,  stimulated  echoes  have  been  successfully  applied  by 

Tanner  to  the  measurement  of  translational  self-diffusion 

coefficients  CTa70D,  by  Lausch  and  Spiess  to  study 

infrequent  Jumps  of  complex  molecules  CLaBO,  SpBOD ,  and 

more  recently  to  analyze  sloui  rotational  motions  of 

molecular  solids  by  Sullivan  et  al .  CSuBBU .  Furthermore, 

other  investigators  conducting  research  into  stimulated  echo 

NriR  imaging,  concurrently  luith  this  investigation,  have 

recently  reported  their  initial  findings  CFrBB,  HaBSD  . 

This  investigation  is  unique  in  that  it  is  the  first  to 

indicate  that  stimulated  echoes  may  be  applied  to  NflR 
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imaging,  to  specifically  outline  hoiu  the  stimulated  echo  may 
be  applied,  and  tc  present  actual  images  utilizing  the 
specified  methods. 

In  NflR  imaging,  image  contrast  from  area  to  area  in  the 
object  results  predominantly  from  the  differences  of  the 
spin  density,  the  spin-lattice  relaxation  time,  Tl ,  and  the 
spin-spin  relaxation  time,  TE.  With  current  instruments, 
contrast  due  to  relaxation  is  achieved  through  the  use  of 
either  the  spin  echo  CSE)  technique,  or  an  inverting 
technique,  such  as  the  inversion  recovery  CIRD  sequence.  In 
particular,  Tl  lueighting  is  introduced  in  the  SE  sequence  by 
the  rapid  repetition  of  the  entire  pulse  sequence,  resulting 
in  signal  saturation,  while  the  IR  technique  introduces  Tl 
weighting  by  inverting  the  equilibrium  magnetization, 
initially  aligned  along  the  positive  z  axis,  and  sampling 
its  recovery  with  a  rt/E   rf  pulse. 

STE  imaging  introduces  Tl  weighting  into  the  NMR  image 
in  the  following  manner.  Uiewed  from  the  rotating  frame  of 
reference,  an  initial  tt/B  rf  pulse,  at  time  equal  to  zero, 
rotates  the  equilibrium  magnetization  into  the  transverse 
plane.  While  in  the  transverse  plane,  the  net  magnetization 
is  reduced  due  to  TE  relaxation,  molecular  diffusion  and 
precession  within  an  inhomogeneous  magnetic  field.  A  second 
w/S  rf  pulse,  at  time  equal  to  t1 ,  will  split  the  net 
magnetization  equally  into  two  orthogonal  components,  one  of 
which  lies  in  the  transverse  plane  and  the  other  which  lies 
in  the  longitudinal  plane. 
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The  individual  isachromats  which  comprise  the  net 
transverse  magnetization  luill  constructively  interfere  to 
form  the  primary  echo  CPE)  at  a  time  equal  to  twice  t1 .  The 
net  longitudinal  magnetization  will  be  reduced  due  to  Tl 
relaxation  and  molecular  self-diffusion.  At  a  time  tE  after 
the  second  SO  degree  rf  pulse,  a  third  SO  degree  pulse  is 
applied  which  rotates  this  Tl  reduced  net  longitudinal 
magnetization  back  into  the  transverse  plans,  where  the 
individual  isochromats  constructively  interfere  to  form  the 
STE  at  a  time  t1  after  the  third  SO  degree  rf  pulse.  It  is 
precisely  this  ability  to  store  and  retrieve  magnetization 
along  the  longitudinal  direction,  where  Tl  relaxation 
occurs,  which  makes  the  application  of  the  STE  to  NMR 
imaging  unique. 

Conventional  Fourier  NMR  imaging  CEdBOD  relies  upon 
spin  echo  formation  for  data  acquisition.  This  investigation 
was  unique  in  that  it  introduced  the  use  of  the  STE  for  data 
acquisition  in  NMR  imaging.  STE  imaging,  with  its  unique  Tl 
dependence,  is  an  ideal  technique  for  Tl  contrast  imaging. 
As  indicated,  there  are  two  viable  methods  of  Tl  contrast 
imaging  currently  in  widespread  use,  the  SE  and  IR 
techniques.  Although  the  IR  sequence  has  produced  excellent 
results,  there  are  a  number  of  distinct  drawbacks  to  its 
implementation.  For  example,  one  must  insure  a  proper 
inverting  pulse  and  use  phase  sensitive  reconstruction  to 
fully  exploit  the  dynamic  range  afforded  by  the  technique. 
The  SE  sequence  is  intrinsically  a  IS   dependent  technique. 
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hence  imaQBS  acquired  luith  this  sequence  mill  frequently 
contain  a  high  degree  of  mixed  Tl  and  TS  contrast.  The 
results  of  this  investigation  indicate  that,  in  many  ways, 
STE  imaging  bridges  the  gap  betuieen  the  accuracy  of  the  IR 
technique,  and  the  efficiency  of  the  SE  imaging  technique. 

The  results  of  this  investigation  took  the  form  of 
specific  applications  of  the  STE  to  NtlR  imaging.  First,  it 
was  shown  that  in  addition  to  generating  Tl  contrast  images, 
it  was  possible  to  calculate  quantitative  Tl  information 
from  a  series  of  STE  images  in  which  the  storage  time  had 
been  systematically  varied.  Second,  a  novel  application  of 
the  Tl  weighted  STE  image  was  obtained:  the  enhancement  or 
suppression  of  elements  in  the  object  with  different  Tl 
values.  Third,  it  was  demonstrated  that  the  STE  was  easily 
integrated  into  chemical  shift  imaging  schemes.  Fourth,  two 
STE  imaging  methods  were  developed  which  permitted  the 
acquisition  of  a  series  of  STE  images  within  one  imaging 
sequence,  where  each  image  was  progressively  weighted  by 
increasing  Tl  relaxation  damping.  Finally,  a  method  of  in 
vivo  determination  of  molecular  translational 
self-diffusion  coefficients,  which  utilized  the  STE's  unique 
Tl  dependence,  was  proposed. 

Theory 
Introduction 

Echo  phenomena  have  long  held  a  prominent  role  in 
spectroscopy,  with  applications  in  various  fields  spanning 
from  magnetic  resonance  to  laser  spectroscopy.  Echoes  were 
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detBctBd  in  NHR  for  the  first  time  in  1950  by  Hahn  CHaSOD , 
and  spin  echoBS  havs  subsequsntly  bBBn  applisd  in  NMR  to 
various  snds,  including  the  mBasursmant  of  Tl  valuas 
CCa54D,  thB  invastigation  into  molacular  diffusion 
procBSSBS  CStBBD,  tha  dBtarmination  of  scalar  coupling 
constants  CFr75D,  the  indirect  dataction  of  magnatic 
rssonancB  CEmBOD,  cohBrsnce  transfar  Cna7B3,  and  NMR 
Imaging  CEdB03 .  Also,  the  sama  effects  have  been  exploited 
in  electron  spin  resonance  Cf1i723,  fnicrowave  spectroscopy 
CG17BD,  and  in  laser  optical  spectroscopy  CKuB4D . 

An  echo  is  usually  created  by  exciting  the  system  under 
investigation  at  least  tiuice,  where  the  excitation  is  often 
pulses  of  electromagnetic  radiation.  All  species  in  the 
system  experience  the  same  inital  pulse,  hence  a  coherence 
is  produced.  In  time,  inhomogeneous  interactions  within  the 
system  act  to  destroy  the  coherence.  This  is  accomplished  in 
NflR  by  an  inhomogeneous  magnetic  field,  by  an  inhomogeneous 
Stark  field  in  microwave  spectroscopy,  or  by  the  Doppler 
effect  in  optical  spectroscopy  Cna7BD .  A  second  pulse, 
applied  at  a  time  t,  inverts  the  accumulated  effects  of  the 
inhomogeneous  interaction.  Thus,  the  initial  coherence  is 
regained  and  an  echo  occurs  at  a  time  Bt .  Under  particular 
conditions,  portions  of  the  coherence  will  continue  to 
defocus,  even  after  the  application  of  the  second  pulse,  and 
hence  will  not  participate  in  echo  formation.  This  component 
is  dubbed  the  narcissus,  after  Narcissus  in  Breek  mythology, 
who  refused  tha  love  of  Echo  CHaBBD . 
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In  1946,  Bloch  introducBd  a  phenomenological  vector 
aquation  to  describe  NMR  CB146D .  It  accurately 
characterized  an  isolated  particle  of  spin  1/E  in  the 
presence  of  a  static  magnetic  field.  Feynman  et  al.  CFb57D 
have  demonstrated  that  this  description  is  complete,  that 
is,  a  geometric  representation  of  the  Schroedinger  equation 
is  possible.  Also,  Pegg  et  al .  CPeSlD  have  shown  this 
description  to  be  perfectly  rigorous,  because  the  vectors 
are  equivalent  to  Heisenberg  operators  in  the  Heisenberg 
representation  of  quantum  mechanics.  It  uias  advantageous  to 
utilize  this  graphical  method  of  analysis,  and  results 
derived  in  this  manner  are  correct  without  restriction. 

Bloch 's  model  of  NflR  assumed  that  the  magnetization  of 
bulk  material,  influenced  by  a  magnetic  field,  conformed  to 
the  laws  of  classical  electrodynamics.  Based  on  this 
premise,  a  vector  differential  equation  was  developed 
relating  the  bulk  magnetic  moment  vector,  fl,  to  the  applied 
magnetic  vector  field,  B,  such  that 

dfl/dt  -  rcn  X  B)  C4-1D 

where  r  is  a  proportionality  constant  called  the 
gyromagnetic  ratio.  The  geometric  interpretation  of  equation 
4-1  is  that  the  magnetic  moment  rotates  about  the  applied 
magnetic  field  with  the  frequency  n,  such  that 

|n|  -  |rB|  .  C4-S3 
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This  relationship  between  the  precession  Frequency  and  the 
applied  magnetic  Field  is  reFered  to  as  the  Larmor  equation. 
The  Frequency  n  is  the  Larmor  Frequency .  The  Larmor 
equation,  expressed  as  in  equation  4-2,  indicates  that  the 
precBssional  Frequency  n  is  proportional  to  the  magnetic 
Field  B,  where  T  is  the  constant  oF  proportionality. 

The  Larmor  equation  may  be  obtained  From  an  argument 
based  upon  classical  physics,  as  outlined  here,  or  derived 
in  identical  Form  From  a  quantum  mechanical  argument.  This 
unique  property  indicates  why  the  classical  Formulation 
oFFers  added  insight  into  the  NflR  phenomenon.  Additionally, 
the  absence  oF  the  Planck  constant  in  the  Larmor  equation 
given  in  equation  4-2  Further  justiFies  the  classical 
treatment  oF  the  resonance  phenomenon . 

Further  support  oF  the  classical  Formulation  oF  the  NMR 
phenomenon  is  oFFered  by  the  correspondence  principle  oF 
quantum  mechanics  CUJaSBD,  which  is  based  on  the  assumption 
that  quantum  theory,  or  at  least  its  Formalism,  contains 
classical  mechanics  as  a  limiting  case.  This  idea  was  First 
expressed  by  Planck  CPIOBD,  when  he  showed  that  in  the 
limit  the  Planck  constant  approaches  zero,  all  quantum 
theoretic  conclusions  converge  towards  classical  results. 
Formation  oF  the  Primary  Echo  Image 

Consider  a  spin  system  in  thermal  equilibrium  with  its 
surroundings,  subjected  to  the  rF  pulse  and  magnetic  Field 
gradient  experiment  displayed  in  Figure  4-1 .  In  the 
graphical  representation  oF  Figure  4-2a,  the  initial 
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Figure  4-1.  Basic  stimulated  echo  imaging  sequence, 
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Figure  4-2.  The  formation  of  a  primary  echo 
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equilibrium  magnetization  of  the  spin  system,  Mi ,  is 
depicted  as  being  initially  aligned  along  the  z  direction, 
coaxial  with  the  main  magnetic  Field.  The  net  magnetization 
is  rotated  into  the  plane  transverse  to  the  main  magnetic 
field  by  the  transmission  CTx)  of  a  90  degree,  or  w/E  rf 
pulse  applied  as  shown  in  figure  'i-Eb .  It  is  assumed  that 
the  rf  pulse  is  of  frequency  nCp)  and  width  tCp),  such  that 
tCp)  is  small  compared  to  Tl  and  IS,  and  excites  the  entire 
chemical  shift  frequency  bandwidth  equally.  If,  for  example, 
the  rf  pulse  has  a  phase  ft   equal  to  90  degrees  Ci.e.  along 
the  positive  x  direction  as  depicted  in  figure  4-Eb5,  then 
the  transverse  magnetization  will  initially  be  aligned  with 
the  negative  y  direction,  in  the  rotating  frame  of 
reference . 

The  rotating  frame  of  reference  refers  to  a  set  of  axes 
which  are  rotating  about  the  z  axis,  the  direction  of  the 
main  applied  magnetic  field.  The  z  axis  of  the  rotating 
frame  is  parallel  to  the  z  axis  of  the  laboratory  frame  of 
reference,  as  defined  by  the  main  applied  magnetic  field. 
The  two  rotating  axes  orthogonal  to  the  z  axis  rotate  with 
an  angular  speed  equal  to  the  effective  component  of  the  rf 
magnetic  field.  The  rotating  axis  in  the  direction  of  the  rf 
magnetic  field  is  refered  to  as  the  in-phase  component, 
while  the  other  axis  is  the  out-of-phase  component.  The 
rotating  frame  of  reference  is  a  useful  construction  for  two 
reasons.  One  reason  is  that  the  phenomenological  vector 
equation  which  describes  NHR  takes  on  a  simpler  form  when 
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expressed  In  terms  of  the  rotating  frame  of  reference.  Also, 
a  simpler  physical  picture  of  events  is  passible  uihen 
considered  in  the  rotating  frame  of  reference. 

For  a  given  nucleus,  theory  indicates  that  resonance 
occurs  at  a  single  frequency,  dictated  by  the  main  magnetic 
field  strength  as  in  equation  4-2.  In  practice,  resonance 
takes  place  over  a  range  of  frequencies  determined  by  the 
inhomogeneity  of  the  main  magnetic  field  throughout  the 
sample.  ThereforB,  the  object  may  be  considered  to  be 
comprised  of  an  ensemble  of  magnetic  moments,  uihoss 
resonance  frequencies  are  symmetrically  distributed  about 
the  Larmor  frequency,  o.  Figure  4-Ec  illustrates  the  free 
precession  Ci.e.  Larmor  type  precession  and  not  rf  pulse 
induced  rotation)  of  all  these  isochromats  during  the  time 
interval  t1 .  Since  figure  ^-Ec   depicts  the  dynamics  of  the 
magnetization  in  a  frame  of  reference  rotating  at  frequency 
o,  the  isochromatic  moment  pairs  maintain  a  symmetry  about 
the  y  direction,  but  rotate  in  opposite  directions.  This  is 
indicated  in  figure  H.-Ec    in  the  following  manner.  The  light 
gray  regions  represent  a  range  of  isochromats  which  deviate 
less  from  a   than  do  the  range  of  isochromats  represented  by 
the  dark  gray,  that  is,  the  light  gray  region  processes 
slower  than  the  dark  gray  region,  in  this  frame  of 
reference.  In  each  case,  solid  area  versus  hatched  area 
indicates  positive  versus  negative  deviation  from  q. 

Within  this  time  interval  t1,  pulsed  linear  magnetic 
field  gradients  are  applied  as  in  the  conventional  fourier 
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imaging  technique  CEdBO] .  The  preparatory  readout  gradient 
is  embodied  in  the  effective  x  gradient,  CGx),  whereas  Glj 
and  Gz  are  employed  for  phase  encoding.  The  effect  of  these 
pulsed  magnetic  field  gradients  is  the  spatial  encoding  of 
the  NHR  signal  in  a  precise  manner.  If  one,  two,  or  three 
orthogonal  pulse  field  gradients  are  used  for  image 
formation,  then  a  one,  two  or  three  dimensional  Fourier 
transform  of  the  time  domain  NflR  signal  will  result  in  an 
image  where  signal  intensity  is  a  function  of  one,  two,  or 
three  spatial  dimensions.  Obviously,  z  direction 
discrimination  could  also  be  achieved  with  a  selective  v/S 
rf  pulse  applied  in  the  presence  of  a  slice  selective  Gz 
CHo77D . 

After  the  time  interval  t1 ,  a  second  30  degree  rf  pulse 
is  applied,  as  indicated  in  figure  4-Sd .  Whereas  prior  to 
the  second  rf  pulse  all  magnetization  was  lying  in  the 
transverse  plane  (assuming  negligible  Tl  relaxation  during 
the  interval  t1),  after  the  second  rf  pulse  the  net 
magnetization  has  components  in  the  longitudinal  plane  as 
well  as  the  transverse  plane.  Figure  4-Ee  illustrates  the 
transverse  component  of  the  net  magnetization  immediately 
following  the  second  rf  pulse,  obtained  by  simply  projecting 
the  net  magnetization  of  figure  4-2d  onto  the  transverse 
plane.  Since  prior  to  the  second  rf  pulse  all  magnetization 
was  lying  in  the  xy  plane,  all  magnetization  after  the 
second  pulse  lies  in  the  xz  plane  Cif  the  rf  pulse  is 
applied  along  the  positive  x  direction  as  indicated),  honce 
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the  projection  of  the  net  magnetization  onto  the  transverse 
plane,  immediately  following  the  second  rf  pulse,  lies 
completely  on  the  x  axis. 

The  intrinsic  properties  of  the  isochromats  have  not 
been  altered  by  this  magnetization  gymnastics.  The  sense  and 
speed  of  free  precession  in  the  transverse  plane  for  the 
isochromats  following  the  second  rf  pulse  is  identical  to 
the  sense  and  speed  of  free  precession  prior  to  the  second 
rf  pulse  Ci.B.  as  indicated  in  figure  4-2c) ,  Hence,  from 
time  t1  on,  the  isochromats  will  freely  process  as  in  figure 
4-Sf .  From  time  t1  to  time  2t1  the  isochromat  vectors  will 
interfere  amongst  themselves,  with  maximum  constructive 
interference  occuring  at  time  StI .  This  constructive 
interference  constitutes  a  primary  echo  CPED,  with  maximum 
amplitude  at  time  2t1,  and  so  named  to  distinguish  it  from 
the  spin  echo  which  results  from  a  w/2-t— it  rf  pulse 
sequence.  The  maximum  amplitude  of  the  PE  at  time  EtI, 
nCPE),  is  given  by 

S 
nCPED  -  nisineisin  CeB/EJexpC-Srl/TEDf CG, D, t1)       C4-3: 

where  fli  is  the  equilibrium  magnetization,  9i  is  the  tip 

angle  of  the  i th  pulse,  and  fCG,D,Tl)  corresponds  to  the 

diffusional  damping  resulting  from  molecular  diffusion  in 

the  presence  of  magnetic  field  gradients.  For  a  constant 

steady  magnetic  field  gradient,  fCG,D,TlD  is  given  by 

E  E   3 
Bxp{!-CE/3)Dr  6  Tl  ;},  where  G  is  the  magnetic  field  gradient 
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and  D  is  the  translational  self  diffusion  coefficient.  It 
should  be  noted  that  pulsed  magnetic  field  gradients  have 
been  used  in  this  investigation,  hence  the  functional  form 
of  fCG.D.Tl)  will  be  different  CSt65D .  The  two  cases  are 
identical,  in  the  limit,  as  use   pass  from  pulsed  to 
continuous  application  of  the  gradient. 

As  illustrated  in  figure  4-1,  the  Gx  readout  gradient 
is  imposed,  centered  about  the  time  Sri,    to  frequency  encode 
the  PE  (iiith  x  direction  spatial  dependence.  Additionally, 
the  receiver  CRx)  is  gated  open  during  this  same  time,  to 
permit  acquisition  CA)  of  the  spatially  encoded  PE.  If  all 
pulses  are  ideal  w/S  rf  pulses,  then  the  PE  image  is 
identical  to  the  image  produced  by  conventional  spin  echo 
imaging,  except  for  a  factor  of  one  half  in  signal 
intensity.  This  reduction  in  the  signal  to  noise  ratio  would 
be  intolerable  unless  it  is  possible  to  recover  it,  or  reap 
some  compensating  benefit.  Fortunately  the  other  half  of  the 
magnetization  is  not  dissipated,  rather  it  has  been  stored 
as  longitudinal  magnetization  by  the  second  tt/S  rf  pulse. 
Formation  of  the  Stimulated  Echo  Image 

It  can  be  shown  that  the  solution  to  the  Bloch  equation 
in  the  rotating  frame  CB14BD  takes  on  the  form 

MxCt+tp)  -  nxCt)  C4-4J 

MyCt+tp)  -  nyCt)cose  -  HzCtDsinG  C4-53 

nzCt+tp)  -  HyCt^sine  -•■  HzCtDcose  C4-BD 
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in  response  to  a  rf  pulse  about  the  x  axis  commencing  at 
time  t  and  of  uidth  tp,  corresponding  to  a  tip  angle  of  9 
degrees.  In  this  representation  the  rotating  portion  of  the 
net  magnetization  is  decomposed  into  two  orthogonal 
components,  fix  is  taken  to  be  in  phase  with  respect  to  the  x 
axis  rotating  frame  of  reference,  while  fly  is  90  degrees  out 
of  phase.  Df  prime  interest  is  equation  4-B,  which 
characterizes  the  longitudinal  magnetization,  and  in 
particular  its  dependence  on  flyCt).  Since  the  spin  system  is 
initially  in  thermal  equilibrium,  f1yCt-0)  -  0.  Hence,  by 
equation  4-6,  HzCt-tp)  -  0  for  an  ideal  tt/E  rf  pulse.  As 
applied  to  the  experiment  of  figure  4-1,  the  implication  is 
that  during  the  interval  t1,  flz  simply  approaches  the 
equilibrium  magnetization,  if  the  affect  of  relaxation  is 
considered . 

If  t1  is  of  the  order  of  TE  or  less,  then  flyCt-Tl)  is 
surely  nonzero.  That  is,  at  time  t1  we  will  have  appreciable 
transverse  magnetization.  Therefore  the  second  ir/E  rf  pulse, 
in  addition  to  inducing  the  PE,  will  also  produce  net 
longitudinal  magnetization,  that  is,  flzCt-Tl+tpD  is 
nonzero.  This  becomes  quite  apparent  when  a  graphical 
analysis  is  conducted. 

Figure  4-3a  is  simply  the  graphical  representation  of 
the  net  magnetization  immediately  following  the  second  w/S 
rf  pulse  of  figure  4-1.  Indeed,  figure  4-3a  is  identical  to 
figure  4-Sd.  Whereas  we  considered  the  transverse  projection 
of  this  net  magnetization  in  order  to  describe  the  formation 


BB 


Figure  4-3.  The  formation  of  a  stimulated  echo 
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of  the  PE,  we  now  consider  the  longitudinal  component  of  the 
net  magnetizatian  immediately  following  the  second  rf  pulse, 
illustrated  in  figure  4-3b .  For  the  duration  of  tB  the 
longitudinal  magnetization  is  affected  solely  by 
spin-lattice  relaxation.  Furthermore,  even  the  readout 
pulsed  magnetic  field  gradient  for  the  PE  image  does  not 
influence  the  longitudinal  magnetization. 

The  third  90  degree  pulse,  which  comes  at  the  end  of 
the  tS  interval,  simply  rotates  the  stored  longitudinal 
magnetization  back  to  the  transverse  plane,  as  depicted  in 
figure  4-3c .  The  intrinsic  properties  of  the  isochromats 
have  not  been  altered  by  the  additional  magnetization 
gymnastics.  Indeed,  the  sense  and  speed  of  free  precession 
in  the  transverse  plane  for  the  isochromats  fallowing  the 
third  rf  pulse  is  identical  to  the  sense  and  speed  of  free 
precession  prior  to  the  second  rf  pulse,  which  is  indicated 
in  figure  4-Ec .  Hence,  after  the  application  of  the  third  rf 
pulse,  the  isochromats  will  freely  precess  as  in  figure 
4-3c .  Since  the  only  time  this  magnetization  was  influenced 
by  the  inhomogeneous  main  magnetic  field  was  during  the 
interval  t1 ,  the  isochromat  vectors  will  interfere  amongst 
themselves  during  the  interval  tE,  with  maximum  constructive 
interference  occur ing  at  a  time  t1  after  the  third  rf  pulse, 
for  tE>t1.  This  constructive  interference  constitutes  the 
stimulated  echo  CSTE),  with  maximum  amplitude  at  a  time  t1 
past  the  application  of  the  third  rf  pulse.  Thus,  the 
maximum  amplitude  of  the  STE  at  a  time  t1  after  the  third  rf 
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pulSB,  nCSTE),  is  given  by 

nCSTE)  -  l/Enisineisine2sine3BxpC-CSTl/TE  +  tE/TDD  • 

FCG,D,t1,t2)  C4-7D 

uihBrs  Oi  is  the  i  th  rf  pulse,  and  fCB.D.Tl.rS)  is  the 

diffusional  damping  term  for  the  STE .  For  a  constant  steady 

magnetic  field  gradient,  fCG,D,Tl,T2)  is  given  by 

S  E   3   E 
expC-CS/3)Dr  G  CtI  +t1  tE)>,  and  is  modified  for  pulsed 

magnetic  field  gradients  CTa703. 

The  relaxation  damping  term  tells  the  history  of  the 
magnetization  that  went  into  the  STE's  formation.  Since 
spin-lattice  relaxation  occurred  only  within  the  interval 
tS,  the  magnetization  must  have  been  stored  along  the 
longitudinal  direction  during  that  interval.  Likewise,  the 
magnetization  can  be  traced  to  the  transverse  direction  for 
both  the  t1  interval  between  the  first  and  second  rf  pulse, 
and  also  for  a  time  t1  subsequent  to  the  third  rf  pulse,  for 
a  total  time  of  TE  influence  amounting  to  EtI . 

As  illustrated  in  figure  4-1,  the  Gx  readout  gradient 
is  imposed,  centered  about  a  time  t1  after  the  third  rf 
pulse,  to  frequency  encode  the  STE  with  x  direction  spatial 
dependence.  Additionally,  the  receiver  CRx)  is  gated  open 
during  this  same  time,  to  permit  acquistion  CA)  of  the 
spatially  encoded  STE.  If  each  0i  was  an  ideal  rr/S  rf 
pulses,  and  ignoring  relaxation  and  diffusional  damping,  it 
is  noted  from  equation  4-7  that  (ICSTE)  is  proportional  to 
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Cl/E)f1i  .  This  is  the  other  factor  of  one  half  ue  noted 
earlier  after  the  formation  of  the  PE .  Whereas  the 
conventional  spin  echo  imaging  experiment  yields  a  single 
image  uihose  intensity  is  proportional  to  Mi  ,  the  stimulated 
echo  imaging  sequence  may  yield  ttuo  images,  each 
proportional  to  Cl/E)ni ,  and  each  in  spatial  registration 
with  the  other.  The  utility  of  these  images  lies  not  in  this 
proportionality,  but  rather  in  the  unique  Tl  dependence  of 
the  STE  image.  Applications  which  further  extend  and  exploit 
this  Tl  dependence  are  presented  later  in  this  chapter. 
Formation  of  Secondary  Echoes 

As  bias  previously  outlined,  the  rf  pulse  sequence  given 
in  figure  4-1  will  yield  both  the  primary  echo  and  the 
stimulated  echo.  Additionally,  the  application  of  the  three 
Tr/2  rf  pulses  may  yield  up  to  three  other  secondary  echoes, 
for  a  total  of  five  echoes  resulting  from  three  rf  pulses. 
The  origins  of  these  secondary  echoes  are  as  fallows. 

The  same  echo  formation  mechanism  which  results  in  the 
primary  echo  after  the  application  of  the  first  two  rf 
pulses  may  also  cause  the  formation  of  two  of  the  three 
secondary  echoes.  If  we  consider  the  three  rf  pulses  taken 
two  at  a  time,  then  there  are  three  unique  combinations,  the 
first  and  second  pulses,  the  first  and  third  pulses,  and  the 
second  and  third  pulses.  Each  combination  will  result  in  an 
echo,  with  the  first  case  simply  being  the  primary  echo.  The 
second  case  results  in  an  echo  at  a  time  t1+tS  after  the 
third  rf  pulse,  and  the  third  case  results  in  an  echo  at  a 
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time  t2  after  the  third  rf  pulse.  Each  of  these  echoes  will 
have  different  T2  weighting,  dependent  upon  the  time  the 
magnetization  spends  in  the  transverse  direction. 

The  final  secondary  echo  is  derived  from  the  PE.  At 
time  t1  after  the  second  rf  pulse  the  PE  has  an  amplitude 
maximum.  That  is,  at  this  time,  the  coherence  imparted  by 
the  first  rf  pulse  has  been  regained.  Indeed,  there  is  as 
much  coherence  amongst  the  isochromats  which  comprise  the  PE 
at  time  t1  after  the  second  rf  pulse  as  there  was 
immediately  following  the  first  rf  pulse,  ignoring  the 
coherence  lost  due  to  relaxation  and  the  diffusion  process. 
Therefare,  for  all  intents  and  purposes,  the  PE  at  time  t1 
after  the  second  rf  pulse  behaves  as  if  it  was  transverse 
magnetization  Just  after  the  application  of  a  rf  pulse. 
Thus,  after  the  time  t1  past  the  second  rf  pulse,  the 
isochromats  which  comprise  the  PE  will  begin  to  lose 
coherence  until  the  third  rf  pulse  acts  to  refocus  them, 
resulting  in  an  echo  at  a  time  t2-t1  after  the  third  rf 
pulse . 

Thus,  the  application  of  three  rf  pulses  may  result  in 
as  many  as  five  echoes.  Each  of  these  echoes  could  be  used 
for  image  formation,  although  the  simplicity  of  the  PE  and 
the  unique  Tl  dependence  of  the  STE  make  these  two  echoes 
ideal  for  image  formation.  The  addition  of  a  fourth  rf  pulse 
results  in  the  formation  of  eighteen  echoes,  and  an  example 
based  upon  this  extended  sequence  is  given  late  in  the 
chapter . 
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flaterials  and  Methods 
The  NHR  Imaner/Spectromster  Sustems 

All  proton  NMR  Bxperiments  uiere  conducted  on  one  of  two 
NtlR  imaBBr/spactrometer  systems.  The  majority  of  experiments 
were  performed  on  a  Nicolet  CGE-NflR)  NT-BO  spectrometer 
located  in  the  Department  of  Radiology  at  the  University  of 
Florida,  ft   single  experiment,  the  TART  imaging  sequence,  was 
conceived  and  developed  at  the  University  of  Florida,  but 
uias  implemented  on  a  prototype  General  Electric  CSI 
imager/spectrometer  located  in  the  applications  lab  of 
General  Electric  NflR  Instruments,  Fremont,  California. 

The  Nicolet  CGE-NMR)  NT-80  spectrometer  mas  coupled 
with  an  Oxford  80/310  MR  superconducting  magnet  operated  at 
1.B9T,  with  a  clear  bore  diameter  of  31  centimeters. 
Transient  gradient  pulses  were  generated  with  the  Oxford 
2320  room  temperature  shim  power  supply,  which  was  modified 
by  J.R.  Fitzsimmons  and  R.G.  Thomas  so  as  to  improve  its 
temporal  response  time.  The  power  supply  was  under  the 
direct  control  of  the  Nicolet  293C  pulse  programmer,  so  as 
to  allow  pulse  control  of  the  x,  y,  and  z  gradients. 
Controlling  software  was  written  in  Nicolet  1280  assembly 
language  by  T.H.  flareci ,  H.D.  Cockman,  and  R.G.  Thomas.  The 
DAC  outputs  to  X,  y,  and  z  were  dependent  upon  the  state  of 
six  TTL  inputs,  two  for  each  direction  so  as  to  allow  pulse 
sequences  to  be  implemented  which  required  the  sign  of  the 
magnetic  gradient  pulse  to  be  reversed  in  the  course  of  the 
experiment.  These  six  TTL  inputs,  U,  A,  U,  C,  W,  and  H 
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Ccorresponding  to  x,  y,  and  z  respect ively)  could  be  gated 
on  during  any  pulse  or  time  delay  in  the  imaging  pulse 
sequence. 

Imaging  pulse  sequences  ujere  generated  directly  with 
the  available  pulse  sequence  generating  software  resident  on 
the  Nicolet  system,  once  the  pulsed  magnetic  gradient 
controlling  hardiuare/sof tiuare  uias  implemented.  Additional 
assembly  language  soFtojare  hias  developed  by  n.D.  Cockman  to 
permit  the  sequential  incrementation  of  the  pulsed  magnetic 
field  gradients,  required  for  the  Fourier  imaging  technique 
used  in  this  investigation. 

The  Nicolet  spectrometer  uses  a  tuienty  bit  word  length, 
hence  single  precision  data  acquisition  resulted  in  twenty 
bit  deep  data.  In  conjunction  with  fl.D.  Cockman,  an  assembly 
language  routine  was  developed  which  scaled  the  twenty  bit 
deep  data  to  five  bits  (permitting  image  display  with 
thirty-two  gray  levels),  and  transfered  the  scaled  data  over 
a  RS23E  transmission  line  to  a  Cromemco  Z-ED  computer  for 
image  processing  and  display .  The  image  processing  and 
display  software  was  developed  by  L.T.  Fitzgerald.  All 
images  presented,  except  those  associated  with  the  TART 
experiment,  were  acquired  and  processed  in  a  like  manner, 
unless  otherwise  noted.  The  images  generated  with  the  TART 
experiment  were  acquired  and  processed  with  the  resident 
hardware  and  software  which  was  standard  on  the  General 
Electric  CSI  imager/spectrometer . 
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Residual  Gradients 

As  Dutlined  in  the  previous  section,  the  stimulated 
echo  imaging  scheme  of  figure  4-1  has,  as  its  backbone,  the 
rf  pulse  sequence  Tf/2-Tl-«/2-T2-iT/2 .  A  PE  image  is  obtained 
at  a  time  t1  after  the  second  Tr/2  rf  pulse,  with  relaxation 
weighting  expC-2Tl/T2) .  The  STE  image  is  obtained  at  a  time 
t1  after  the  third  rr/2  rf  pulse,  with  its  relaxation 
weighting  given  by  expC-2Tl/T2)expC-T2/T15 .  Since  the 
spatial  encoding  is  identical  for  the  PE  image  and  the  STE 
image,  the  two  are  in  spatial  registration,  though  residual 
gradients  may  act  in  a  manner  which  alters  this 
registration,  as  to  be  explained. 

Consider  the  effect  of  residual  gradients  in  the 
following  gedanken  experiment.  Phase  encoding  pulsed 
magnetic  field  gradients  are  applied  within  the  interval  t1 
of  the  basic  STE  imaging  sequence  as  depicted  in  figure  4-1 . 
Ideally  there  would  be  no  residual  gradients  present,  that 
is,  there  would  be  no  field  gradient  once  the  current  pulse 
which  produced  the  phase  encoding  gradient  had  died  down. 
Now  consider  the  effect  of  a  residual  gradient  which  remains 
even  after  the  current  pulse  ceases.  All  isochromats  in  the 
interval  t1  will  experience  a  phase  shift  B  as  a  result  of 
the  phase  encoding  pulsed  gradient.  At  time  t1,  half  the 
magnetization  is  rotated  to  the  longitudinal  direction,  and 
thus  is  no  longer  influenced  by  field  gradients.  The 
remaining  transverse  magnetization  experiences  a  phase 
reversal  of  sorts  at  time  t1 ,  which  accounts  for  the 
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formation  of  the  PE.  If  there  is  a  residual  gradient,  the 
transverse  isochromats  will  rephase  in  its  presence, 
resulting  in  an  additional  phase  shift  of  bICPED  .  The  stored 
longitudinal  magnetization  is  not  influenced  by  the  residual 
gradient  until  it  is  brought  back  to  the  transverse  plane, 
where  it  rephases  to  form  the  STE.  If  a  residual  gradient  is 
present,  there  will  be  an  associated  additional  phase  shift 
of  0CSTE) . 

Both  BlCPED  and  plCSTE)  are  negative  relative  to  B,  for 
they  were  incurred  after  a  phase  reversing  rf  pulse,  and 
|tfCPE)|  is  greater  than  |b1CSTED|,  since  the  magnitude  of  a 
residual  gradient  is  typically  a  monotonically  decreasing 
function  of  time.  Consequently,  the  PE  experiences  a  total 
phase  shift  which  is  less  than  that  of  the  STE  by 
|bjCPE)-0CSTE)  I  .  Since  ideally  there  is  a  one-to-one  mapping 
of  total  phase  shift  to  spatial  location,  the  different 
total  phases  of  the  PE  image  and  the  STE  image  result  in 
different  spatial  mappings  along  the  phase  encoding 
direction.  That  is,  since  the  total  phase  of  the  PE  is  less 
than  that  of  the  STE,  the  apparent  extent  of  the  phase 
encoded  direction  will  be  less  in  the  PE  image  than  in  the 
STE  image.  For  example,  points  separated  by  10  centimeters 
in  the  STE  image,  along  the  phase  encoded  direction,  might 
be  depicted  as  having  only  an  apparent  8  centimeter 
separation  in  the  phase  encoding  direction  of  the  PE  image. 
Since  the  readout  direction  is  essentially  unaffected  by  the 
residual  gradient,  the  aspect  ratio  for  all  STE  images  will 


95 
be.  idealltj,  one,  while  for  the  PE  images  the  aspect  ratio 
mill  vary  as  a  function  of  the  magnitude  of  the  residual 
gradient,  and  the  time  it  is  influential  on  the  transverse 
magnetization.  In  all  images  presented,  except  those 
corresponding  to  the  TART  experiment,  the  phase  encoding 
direction  is  displayed  along  the  horizontal  direction,  or  y 
direction. 

To  evaluate  the  residual  gradient  problem  on  the  NT-BO 
spectrometer,  the  following  experiment  was  conducted.  A 
single  one  centimeter  diameter  vial  of  paramagnetically 
doped  water  was  centered  in  the  rf  coil  of  the  spectrometer. 
The  magnet  was  then  shimed  with  the  sample  in  place, 
resulting  in  an  approximate  10  ppm  line  width.  A  pulse 
sequence  was  developed,  and  implemented,  to  evaluate  the 
duration  and  extent  of  the  residual  gradient. 

The  pulse  sequence  consisted  of  six  steps.  The  first 
step  was  simply  a  time  delay,  set  equal  to  one  second  so  as 
to  insure  the  entire  sequence  was  not  so  rapidly  repeated  as 
to  cause  signal  saturation.  The  second  step  was  the 
application  of  a  pulsed  magnetic  field  gradient  of  0.18 
mT/m,  for  a  duration  of  two  seconds.  This  duration  was 
chosen  so  as  to  insure  that  the  field  gradient  had 
sufficient  time  to  reach  its  full  magnitude.  The  next  step 
in  the  sequence  was  a  variable  timing  delay,  which  was 
incremented  from  ES    milliseconds  to  7E5  milliseconds  in  S5 
milliseconds  increments,  with  a  single  inital  value  of  one 
microsecond.  At  the  end  of  this  variable  time  delay,  the 
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fourth  step  in  the  sequence  cammenced,  a  90  degree  rf  pulse, 
with  the  next  step  simply  the  opening  of  the  receiver  gate. 
The  sixth  and  final  step  in  the  sequence  was  a  time  delay 
equal  to  the  total  acquisition  time,  as  dictated  by  the 
sweep  width  and  the  data  acquisition  block  size.  This 
permitted  the  receiver  gate  to  remain  open  for  the 
acquisition  of  the  free  induction  decay  CFID)  resulting  from 
the  90  degree  rf  pulse. 

The  significance  of  the  experiment  was  as  follows. 
Obviously,  if  there  was  no  pulsed  magnetic  field  gradient 
applied  in  step  two,  the  signal  intensities  of  the  Fourier 
transformed  FIDs  would  not  be  dependent  upon  the  variable 
time  delay  of  step  three.  Similarly,  if  there  was  a  gradient 
applied  in  step  two,  but  its  magnitude  was  reduced  to  zero 
prior  to  the  application  of  the  w/E  rf  pulse,  then  the 
resultant  signal  intensities  from  the  transformed  FIDs  would 
still  not  be  dependent  upon  the  variable  time  delay  of  step 
three.  Rather,  consider  the  case  where  a  gradient  is  applied 
in  step  two,  and  its  magnitude  had  not  reduced  to  zero  prior 
to  the  application  of  the  Tr/2  rf  pulse.  This  is  what  is 
meant  by  a  residual  gradient,  and  in  this  case  the 
corresponding  signal  intensities  of  the  transformed  FIDs 
will  be  reduced.  Thus,  by  incrementing  the  variable  time 
delay  it  was  possible  to  temporally  map  out  the  extent  of 
the  residual  gradient.  The  plot  of  signal  intensity  as  a 
function  of  the  variable  time  delay  is  given  in  figure  4-4. 
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Figure  4-4.  The  evaluation  of  a  residual  gradient, 
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It  was  apparent  From  the  results  of  this  experiment 
that  there  was  a  residual  gradient  present  after  application 
of  a  pulsed  field  gradient,  and  its  effect  uias  significant 
for  nearly  IBO  milliseconds.  Since  the  magnitude  of  the 
field  gradient  used  in  this  experiment  luas  typical  of  the 
values  used  in  the  STE  imaging  experiments,  and  since  the 
time  extent  of  the  residual  gradient  uias  comparable  to  the 
time  scale  used  in  the  experiments,  it  uias  concluded  that 
the  problem  of  residual  gradients  could  not  be  ignored  in 
this  investigation. 

Figure  1-5  illustrates  the  effect  which  the  residual 
gradient  has  upon  image  formation.  As  outlined  earlier,  the 
effect  of  residual  gradients  was  upon  the  spatial  mapping  in 
the  phase  encoding  direction.  The  images  are  of  two,  one 
centimeter  diameter  vials  filled  with  paramagnetically  doped 
water.  The  images  are  displayed  as  intensity  contour  plots, 
rather  that  gray  scale  images.  Horizontal  is  the  phase 
encoding  direction  and  vertical  is  the  readout  direction. 
The  image  in  figure  4-5a  was  formed  from  transverse 
magnetization  in  the  presence  of  a  residual  gradient  for  100 
milliseconds,  while  the  image  in  figure  4-5b  was  in  the 
presence  of  the  residual  gradient  for  BSD  milliseconds,  and 
figure  4-5c  for  350  milliseconds. 

As  predicted,  the  more  time  the  transverse 
magnetization  spends  in  the  presence  of  the  residual 
gradient,  the  less  is  the  apparent  separation  of  the  vials 
in  the  phase  encoding  direction.  Indeed,  in  the  limit  where 


93 


Figure  4-5.  Effect  af  residual  gradients  on  image  formatian, 
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the  integrated  effect  of  the  residual  gradient  just  balances 
that  of  the  phase  encoding  gradient,  there  would  be  no 
effective  phase  encoding  gradient,  and  hence  no  spatial 
separation  in  the  image  along  the  phase  encoding  direction. 
Phase  Cucling 

f^   drauback  to  the  STE  imaging  sequence  as  described  is 
its  insensitivity  to  signed  phase  information.  Because  there 
is  only  pure  amplitude  modulation  of  the  signal  as  a 
function  of  the  gradients  applied  during  the  interval  t1 ,  it 
would  normally  be  necessary  to  set  the  Tx  frequency  outside 
the  spectral  frequency  range,  employ  single  phase  detection, 
and  USB  only  one  sided  phase  encoding  gradients  to  avoid 
aliasing.  Fortunately  a  simple  modification  afforded  by 
phase  cycling  of  the  rf  pulses  and  the  receiver  circumvented 
this  problem  by  converting  the  amplitude  modulation  into 
phase  modulation. 

The  two  step  phase  cycling,  used  in  conjunction  with 
quadrature  detection,  consisted  of  the  first  two  acquistions 
as  depicted  in  table  4-1.  Here,  A,  B,  C,  and  D  refer  to  the 
phase  of  the  rf  pulses  and  the  receiver,  as  illustrated  in 
figure  4-1 .  Nagayamo  et  al .  have  used  a  closely  related 
phase  cycling  procedure  in  spin  echo  correlated  spectroscopy 
:Na73D . 

The  first  two  phase  cycling  steps  select  the  echo 
component  of  the  transverse  magnetization  over  the  narcissus 
component.  If  either  the  main  field  inhomogeneity  or  the 
readout  gradient  was  of  sufficient  magnitude  to  suppress  the 
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Table  4-1 

Four-Step  Phase  Cycling  Used  in 
Stimulated  Echo  Imaging 
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narcissus  component,  then  the  two  step  cycle  would  not  be 
needed,  for  only  the  echo  component  would  be  available  for 
acquisition.  Recall,  the  narcissus  component  is  comprised  of 
magnetization  which  continues  to  defocus  even  after  the 
phase  reversing  pulse  has  been  applied.  The  defocus  of  the 
narcissus  would  be  advanced  if  the  inhomogeneity  of  the 
magnetic  field  is  enhanced,  and  a  magnetic  field  gradient 
may  be  considered  as  Just  a  large  field  inhomogeneity.  Thus, 
the  narcissus  component  would  have  defocused  to  essentially 
zero  when  the  time  came  to  acquire  the  echo.  For  this 
scenario,  the  two  step  phase  cycling  would  be  superfluous. 

Although  the  two  step  phase  cycling  sequence  was 
adequate,  often  it  was  extended  to  four  steps  in  order  to 
suppress  artifacts  resulting  from  longitudinal  magnetization 
recovery  during  the  interval  t1  CBaBlD .  The  four  step 
phase  cycling  is  illustrated  in  its  entirety  in  table  4-1. 
The  phases  of  B  and  C  are  arbitrary,  as  long  as  they  remain 
equal  and  constant  throughout  the  entire  experiment.  For  a 
properly  adjusted  instrument,  and  with  all  rf  pulses  being 
slice  selective,  the  last  two  phase  cycling  steps  in  table 
f-l  could  be  ignored.  Indeed,  under  suitable  conditions,  STE 
imaging  could  be  conducted  with  a  single  step,  without  the 
need  to  phase  cycle. 

In  practice,  an  additional  spoiling  pulsed  magnetic 
field  gradient  was  applied  immediately  following  the  primary 
echo's  readout  gradient,  and  prior  to  the  third  Tr/2  rf 
pulse,  along  Gx .  This  reduced  artifacts  in  the  stimulated 
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acho  imagB  by  eliminating  the  three  other  passible 
overlapping  secondary  echoes  which  may  occur  after  the  third 
rf  pulse.  Since  the  STE's  magnetization  is  stored  along  the 
longitudinal  direction  during  the  time  interval  tB,  the 
spoil  pulse  dispersed  only  the  transverse  magnetization 
uhich  mould  have  gone  on  to  create  additional  secondary 
echoes  subsequent  to  the  third  rf  pulse. 
Image  Acquisition 

Unless  otherwise  noted,  all  images  acquired  with  the 
NT-BO  spectrometer/imager  were  done  so  in  a  like  manner. 
Although  figure  4-1  indicates  that  Gz  phase  encoding  was 
possible  in  the  STE  imaging  sequence,  none  was  performed  in 
practice.  This  resulted  in  z  direction  averaged  images  Cinto 
the  page) .  The  reason  for  this  was  to  save  time  and  data 
storage  space,  and  does  not  reflect  any  inherent  limitation 
of  the  imaging  system.  The  Gx  readout  gradient  was  0.5  mT/m, 
and  was  applied  for  32ms,  There  were  64  equal  steps  of  the 
6y  phase  encoding  gradient,  ranging  from  0.5  to  -0.5  mT/m, 
each  applied  for  16ms,  excluding  any  contribution  from 
residual  gradients.  The  data  matrix  for  each  image  consisted 
of  128  X  64  complex  data  points .  The  f ield-of-view  was  10  X 
10  centimeters,  corresponding  to  a  resolution  of  1.5 
millimeters  in  the  y  direction  and  0.75  millimeters  in  the  x 
direction.  Four  averages  were  done  to  take  advantage  of  the 
four  step  phase  cycling  given  in  table  4-1.  The  w/E  rf  pulse 
width  was  260  microseconds.  All  pulse  sequence  repetition 
times  were  set  so  as  to  avoid  signal  saturation  effects. 
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All  images  displayed  corresponding  to  the  TART 
Bxperiment  ojera  acquired  in  the  folloujing  manner.  The  two 
dimensional  images  uere  acquired  uiith  a  field  of  vietu  of  15 
X  15  centimeters,  and  a  slice  thickness  of  E  millimeters. 
Each  image  forming  echo  was  defined  with  a  data  block  size 
of  512  data  points.  There  were  256  equal  phase  encoding 
steps.  This  resulted  in  a  final  image  matrix  of  512  X  512 
points,  after  zero  filling  in  both  directions. 
Experiments 

Seven  experiments  which  exploited  the  unique  properties 
of  the  stimulated  echo  in  NflR  imaging  were  conducted.  These 
experiments  covered  a  wide  spectrum  of  application. 
Presented  in  this  section  are  the  methods  of  using  STE 
imaging  to  perform  Tl  contrast  imaging,  quantitative  Tl 
imaging,  water-lipid  contrast  imaging,  chemical  shift 
imaging,  multiecho  STE  imaging,  and  diffusion  coefficient 
weighted  imaging. 
Experiment  one:  Tl  contrast  imaging 

To  demonstrate  the  Tl  contrast  imaging  capabilities  of 
the  STE  imaging  sequence,  the  basic  imaging  sequence  given 
in  figure  4-1  was  used,  in  conjunction  with  the  phase 
cycling  scheme  of  table  4-1.  The  image  phantom  consisted  of 
a  red  Delicious  apple  with  a  one  centimeter  vial  of  copper 
sulfate  doped  distilled  water  embedded  in  its  center.  A 
short  t1  value  was  used  to  emphasize  the  poor  Tl  contrast  of 
the  PE  image,  while  the  long  t2  value  emphasized  Tl 
differences  dramatically  in  the  STE  image. 
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Experiment  tujo:  quantitative  Tl  Imaging 

In  addition  to  Tl  contrast  imaging,  STE  imaging  permits 
the  quantitative  measurement  of  Tl  values.  It  can  be  shown 
that  the  diffusional  damping  of  the  STE.  f CG, D, t1 , t2) ,  is 
functionally  dependent  on  t1  to  second  and  third  order,  for 
both  constant  and  pulsed  magnetic  field  gradients  CHaBO, 
Ta70D ) .  Hence,  an  analysis  of  equation  4-7  yields  the 
result  that  for  t1  very  small,  a  plat  of  the  logarithm  of 
the  stimulated  echo's  maximum  amplitude  versus  tE  uiill  yield 
a  straight  line  luhose  slope  is  a  measure  of  Tl .  If  the 
effect  of  diffusion  is  not  completely  negligible,  such  that 
corrections  for  higher  order  in  t1  must  be  made,  then  a 
better  value  of  Tl  can  be  obtained  by  noting  that  for  t1 
much  less  than  t2,  the  stimulated  echo's  maximum  amplitude 
decreases  with  t2  as  expC-rE/TCef f )D ,  where  the  effective 
relaxation  time  TCeff),  is  given  by 

1/TCeff)  -  1/Tl  +  gCB,D,Tl)  C4-8: 

with  gCG.D.Tl)  being  the  diffusional  contribution  to  second 

order   in  t1 .  This  term,  for  a  constant  steady  gradient,  is 

E  e   E 
given  by  gCG,D,Tl)  -  r  G  DtI  ,  while  the  functional  formed 

of  gCG,D,TlD  for  pulsed  magnetic  field  gradients  is 

different,  though  the  t1  squared  dependence  is  maintained 

CTa70D.  Hence,  a  plot  of  1/TCeff)  against  t1  squared 

yields  a  straight  line  whose  intercept  on  the  ordinate 

yields  the  value  of  Tl . 
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Experiment  three:  uiatsr-lipid  contrast  imaninn 

A  novel  application  of  the  Tl  weighted  STE  image  is 
either  the  enhancement  or  the  suppression  of  elements  in  the 
object  with  different  Tl  values.  In  a  conventional  spin  echo 
image  of  an  object  with  various  Tl  elements,  the  longer  Tl 
elements  may  be  suppressed  by  saturation  of  the  signal  via 
rapid  pulse  sequence  repetition  times.  In  contrast,  with  STE 
imaging  it  is  possible  to  suppress  the  shorter  Tl  elements 
via  Tl  weighting.  This  capability  of  STE  imaging  was  clearly 
demonstrated  by  imaging  a  hen's  egg,  where  the  difference  in 
the  Tl  values  between  the  albumen  and  yolk  Ccorresponding  to 
a  water  dense  and  a  lipid  dense  element,  respectively) 
allowed  for  easy  visualization  of  this  STE  image  property. 
Experiment  four:  chemical  shift  imaging 

The  STE  is  easily  integrated  into  chemical  shift 
imaging  schemes.  Any  existing  chemical  shift  imaging 
technique,  which  currently  makes  use  of  echo  formation 
Cna83D,  could  benefit  from  the  unique  properties  of  STE 
imaging.  Chemical  shift  STE  imaging  was  conducted  on  a 
phantom  consisting  of  a  vial  of  water  and  a  vial  of  oil, 
spatially  separated  in  the  y  direction,  but  not  the  x 
direction.  No  readout  gradient  was  employed  so  as  to  retain 
chemical  shift  information,  and  the  x  direction  was  the 
phase  encoded  direction. 
Experiment  five:  extended  STE  multiecho  imaging 

The  three  it/E  rf  pulse  STE  imaging  sequence  given  in 
figure  H-l  was  extended  to  a  Tl  weighted  multiecho  imaging 
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sequence  as  depicted  in  figure  4-6.  The  addition  of  a  fourth 
Tr/2  rf  pulse  results  in  the  formation  of  up  to  thirteen 
additional  tertiary  echoes  CTE5,  six  of  which  have  a  Tl 
dependence  in  their  relaxation  damping  term  CWoBlD . 

The  practical  factors  uihich  entered  into  the  decision 
as  to  uhich  TE  to  use  for  imaging  were  as  follows.  With  the 
possible  production  of  so  many  echoes,  it  was  important  to 
locate  one  which  could  be  isolated  from  the  other  echoes.  By 
the  proper  choice  of  the  time  intervals  t1,  t2,  and  t3,  it 
was  possible  to  create  some  echoes  which  did  not  overlap 
others.  Additionally,  it  was  desired  to  be  able  to  acquire 
the  PE,  the  STE,  and  a  TE  all  within  a  single  pass  of  the 
imaging  sequence.  Hence,  since  the  formation  of  the  PE  and 
STE  require  w/E  rf  pulses,  the  TE  chosen  for  echo  formation 
could  not  have  a  cosine  tip  angle  dependent  term,  as  is  seen 
in  equation  4-9  which  follows. 

Of  primary  concern  was  the  relative  intensity  of  the 
echo.  For  any  given  echo,  there  were  three  factors  which 
acted  to  reduce  the  echo  signal  intensity,  the  tip  angle 
dependence,  the  relaxation  damping,  and  the  damping  due  to 
diffusion  in  an  inhomogeneous  magnetic  field.  Although  there 
were  certain  definite  applications  of  the  tip  angle 
dependence  term,  this  investigation  was  concerned  primarily 
with  exploiting  the  applications  inherent  in  the  relaxation 
damping  term.  Hence,  it  was  advantageous  to  image  with 
echoes  that  had  a  tip  angle  dependence  as  close  to  one  as 
possible.  Although  the  diffusion  damping  term  was  of 
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Figure  4-6.  The  extended  stimulated  echo  imaging  sequence 
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interest,  it  uas  the  form  of  the  relaxation  damping  term 
which  finally  dictated  which  echo  was  to  be  used  for  image 
formation.  It  was  desirable  to  use  a  TE  with  as  simple  a 
relaxation  damping  term  as  possible. 

All  these  conditions  were  fulfilled  by  the  TE  whose 
maximum  occurs  at  a  time  t2-t1  after  the  application  of  the 
fourth  w/H  rf  pulse,  whose  amplitude  at  that  time,  flCTEJ,  is 
given  by 

nCTE5  -  -Cl/EDMisineisin  Ce2/2)sine3sine4  * 

BxpC-CETE/TS+T3/Tl)Df CG,D,t1,t2,t3)        C4-93 

where  (li  is  the  equilibrium  magnetization,  0i  represents  the 

tip  angle  of  the  i th  pulse,  T2  is  the  spin-spin  relaxation 

time,  Tl  the  spin-lattice  relaxation  time,  and 

f CG, D, t1 , t2, t3)  corresponds  to  the  diffusional  damping 

CWoBlD,  where  G  represents  the  magnetic  field  gradients 

and  D  the  translational  self-diffusion  coefficient.  The  time 

intervals  t1,  t2,  and  t3  are  defined  as  in  figure  4-6. 

This  echo  is  a  tertiary  stimulated  echo,  resulting  from 
the  storage  of  the  primary  echo  magnetization  by  the  third 
rf  pulse.  The  fourth  rf  pulse  rotates  this  stored 
magnetization  into  the  transverse  plane  for  echo  formation 
and  detection.  Using  all  tt/2  rf  pulses,  the  PE,  STE,  and  TE 
were  all  produced  and  available  for  image  formation  via  the 
execution  of  a  single  imaging  pulse  sequence. 

In  practice,  the  software  which  was  resident  on  the 
imager/spectrometer  permitted  the  acquisition  of  only  two 
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BchoBS  during  a  single  SBquencB,  hencB  thB  SBquBncs  had  to 
bs  rBpsatad  twice  so  as  to  acquire  all  three  echoBS.  ThB 
entire  sequBnce  tuas  usBd  Bach  tims,  so  this  limitation  did 
not  havB  any  effect  upon  the  results.  This  problem  has  been 
resolved  by  the  recent  addition  of  more  flexible  data 
acquisition  softiuare. 
ExpsrimBnt  six:  TftRT  multiecho  imaging 

The  extended  STE  multiecho  imaging  sequence  resultBd  in 
a  sBriBs  of  imagBs,  each  progressivsly  weighted  with 
spin-lattice  relaxation.  Although  the  underlying  concept  was 
correct,  some  of  the  intrinsic  problems  with  the  technique 
acted  to  limit  its  usefulness.  For  examplB,  with  so  many 
Bchoes  being  gensratsd  in  onB  pass  of  thB  SBqusncB  it  oftBn 
was  difficult  to  insurs  that  there  was  no  overlapping  of 
undasirsd  echoes  with  desirsd  echoes.  This  problem  would  be 
compounded  if  the  sequence  was  extended  to  five,  or  more, 
Tr/S  rf  pulses. 

Fortuitously,  anothsr  mBthod  of  forming  a  SBriBS  of  Tl 
waighted  images  was  suggested  by  earlier  work  in  the 
observation  of  Tl  rslaxation  CFrBBD  and  zsro  quantum 
cohsrBncB  CBaBOD.  In  both  thBSB  casBS,  the  spin  system  was 
manipulated  into  a  given  state,  and  the  subssquent  time 
evolution  of  the  system  was  monitorsd  with  a  SBriBS  of 
sampling  rf  pulses,  sufficiently  weak  so  as  to  minimize 
their  perturbing  effect  on  the  evolving  spin  system. 

It  was  possiblB  to  apply  the  storage  and  retrisval 
propBrtiBS  of  thB  STE  Calong  the  longitudinal  direction)  to 


Ill 

an  advantage  in  muitiecho  Tl  imaging.  Described  earlier,  and 
illustrated  in  figure  4-3,  the  application  of  the  first  two 
Tr/2  degree  rf  pulses  in  the  STE  imaging  sequence  resulted  in 
storing  magnetization  along  the  longitudinal  direction 
Cfigure  4-3b) .  In  the  basic  STE  imaging  sequence,  all  of 
this  longitudinal  magnetization  uias  rotated  into  the 
transverse  plane  after  a  Tl  ueighted  interval  tE,  to  form 
the  STE  image.  As  indicated  in  equation  4-7,  the  echo 
amplitude  is  dependent  upon  the  sine  of  the  third  rf  tip 
angle.  Indeed,  if  an  ideal  tt/E  rf  pulse  was  used,  then  the 
resultant  tip  angle  dependence  was  simply  unity.  Fundamental 
to  the  understanding  of  the  TART  sequence  is  the  fact  that, 
although  the  first  two  rf  pulses  of  the  basic  STE  imaging 
sequence  resulted  in  a  net  longitudinal  magnetization,  there 
was  no  reason  that  this  entire  magnetization  had  to  have 
been  brought  into  the  transverse  plane  all  at  once  by  the 
third  rf  pulse.  Indeed,  it  was  possible  to  bring  only  a 
portion  of  the  magnetization  into  the  transverse  plane  by 
using  tip  angles  less  than  or  equal  to  90  degrees. 

Figure  4-7  is  the  tip  angle  reduced  Tl  CTART)  imaging 
sequence.  In  this  version  of  the  sequence,  slice  selection 
was  determined  by  the  first  90  degree  rf  pulse,  applied  in 
the  presence  of  a  magnetic  field  gradient,  and  was  selective 
perpendicular  to  the  z  axis.  The  phase  of  the  first  rf 
pulse.  A,  and  the  phase  of  the  receiver,  B,  were  cycled  in 
the  same  manner  as  the  rf  pulse  and  receiver  phases  in  the 
basic  STE  imaging  sequence.  Following  convention,  the  time 
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Figure  4-7.  The  tip  angle  reduced  Tl  CTART)  imaging 
sequence . 
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from  the  first  rf  pulse  to  the  maximum  amplitude  of  the  PE 
uas  denoted  as  TE,  thus  the  time  betuieen  the  first  and 
second  rf  pulses  was  TE/S.  All  phase  encoding  uas  performed 
during  this  initial  interval  TE/2,  therefore  the  PE  and  the 
resultant  series  of  STE  images  were  in  spatial  registration. 
In  contrast  with  the  basic  STE  imaging  sequence  where  a 
single  90  degree  rf  pulse  was  used  to  form  a  single  STE 
image,  the  TART  imaging  sequence  used  a  series  of  rf  pulses, 
whose  tip  angles  were  30  degrees  or  less,  to  form  a  series 
of  Tl  weighted  STE  images.  The  reduced  tip  angle  rf  pulses 
were  denoted  by  al ,  aB,  a3,  and  generalized  to  an  in  figure 
4-7.  The  Tl  recovery  time,  from  the  second  90  degree  rf 
pulse  to  the  rf  pulse  with  tip  angle  ai  ,  was  denoted  as  RTi  . 
Thus,  the  i th  STE  image  had  a  Tl  weighting  given  by 
expC-RTi/Tl) . 

The  utility  of  the  TART  imaging  sequence  lies  in  the 
fact  that,  with  the  correct  choice  of  RTi  values,  the  entire 
spin-lattice  relaxation  curve  was  sampled  during  a  single 
pass  of  the  imaging  sequence,  with  each  sampling  point 
corresponding  to  one  of  the  STE  images  in  the  series.  The 
fact  that  each  STE  shared  a  common  TB  dependence,  given  by 
BxpC-TE/TB),  aided  in  simplifying  the  resulting  echo 
amplitude  expressions,  although  the  expressions  were 
complicated  by  the  fact  that  each  STE  had  a  different 
trigonometric  dependence  upon  the  rf  tip  angles. 

This  complication  was  eliminated  in  the  following 
manner.  Figure  4-8  illustrates  the  use  of  four  nonequal. 
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Figure  4-8.  The  formation  of  a  series  of  TART  images 
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reduced  tip  angles  to  generate  a  series  of  four  STE  images 
which  were  identical  in  tip  angle  dependence  and  T2 
relaxation  weighting,  differing  solely  in  Tl  weighting. 
Figure  4-Ba  is  the  longitudinal  component  of  magnetization 
immediately  after  the  second  rf  pulse  of  figure  4-7.  Indeed, 
figure  4-Ba  is  identical  to  figure  4-3b.  If  ideal  rr/B  rf 
pulses  were  used  for  the  first  two  rf  pulses,  and  TS 
relaxation  during  the  interval  t1  is  ignored,  then  this 
component  is  simply  Cl/B)ni  .  Thus  the  application  of  the 
first  tip  angle  reduced  rf  pulse,  with  tip  angle  al,  will 
rotate  into  the  transverse  plane  a  quantity  of  magnetization 
equal  to  Cl/EJflisinal,  leaving  a  quantity  (l/2)nicosal  in 
the  longitudinal  plane.  This  is  illustrated  in  figure  4-Bb . 
The  procedure  is  repeated,  as  in  figure  4-Bc,  when  the 
second  tip  angle  reduced  rf  pulse  with  tip  angle  aE   rotates 
a  quantity  of  magnetization  Cl/S5ni cosalsinaB  into  the 
transverse  plane,  leaving  a  quantity  Cl/B)ni cosalcosaB  of 
magnetization  along  the  longitudinal  direction.  Figures  4-Bd 
and  4-Bb  illustrate  the  application  of  the  third  and  fourth 
reduced  tip  angle  rf  pulses.  This  procedure  may  be 
generalized  to  the  application  of  n  reduced  tip  angle  rf 
pulses,  where  the  quantity  of  magnetization  rotated  into  the 
transverse  plane,  now  including  the  relaxation  terms  but 
ignoring  diffusional  effects,  is  given  by 


n-1 
Sn  -(l/B)nisino(nC  tr  cosai  )expC-CTE/TB  +  RTn/TlDD    C4-10D 

i-1 
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luhara  Sn  is  the  maximum  signal  intensity  of  the  nth  STE .  All 
other  terms  have  been  previously  defined. 

Equation  1-10  is  comprised  of  four  seperate  factors. 
One  factor  is  the  term  dependent  upon  the  total  longitudinal 
magnetization  available  prior  to  the  application  of  any 
reduced  tip  angle  rf  pulses,  simply  Cl/2)ni  .  This  factor  is 
the  same  for  each  STE  image  in  the  series.  A  second  factor 
which  is  the  same  for  each  STE  in  the  series  is  the  T2 
relaxation  damping  term.  A  third  factor,  the  Tl  relaxation 
damping  term,  is  different  for  each  STE  image  in  the  series, 
by  design,  for  this  is  the  source  of  image  to  image 
contrast.  The  final  factor  in  equation  4-10  is  the  tip  angle 
dependence.  Although  fundamentally  different  for  each  STE 
image,  it  was  possible  to  adjust  the  values  of  the  reduced 
tip  angles  in  a  manner  which  resulted  in  identical  tip  angle 
dependence  in  each  STE  image.  This  is  graphically  displayed 
in  figure  l-B,  where  the  tip  angles  are  such  that  a  constant 
amount  of  transverse  magnetization  is  produced  with  each 
application  of  a  tip  angle  reduced  rf  pulse,  ignoring  Tl 
relaxation.  Thus,  when  each  STE  in  the  series  is  formed,  as 
depicted  in  figure  4-8f,  the  only  difference  in  each  image's 
maximum  amplitude  is  the  different  exponential  Tl  weighting. 
Althought  a  qualitative  argument  based  on  figure  4-8 
indicates  that  each  successive  reduced  tip  angle  must  be 
greater  that  the  previous  one  Cfor  the  quantity  of 
longitudinal  magnetization  diminishes  with  each  rf  pulse, 
yet  it  is  desired  to  keep  the  quantity  of  induced  transverse 
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magnetization  a  constant),  it  is  possible  to  quantitatively 
determine  the  specific  values  of  the  reduced  tip  angles 
required. 

There  are  tuio  key  points  to  consider  in  the  development 
of  a  quantitative  expression  for  the  values  of  the  reduced 
tip  angles  which  will  result  in  each  STE  image  of  the  series 
having  the  same  tip  angle  dependence.  The  first  point  is  the 
most  obvious,  it  is  necessary  to  set  the  trigonometric  form 
of  the  tip  angle  dependence  of  the  nth  STE  equal  to  the 
trigonometric  form  of  the  tip  angle  dependence  of  the 
Cn-l)th  STE.  For  example,  if  n  equals  four,  then  setting 
sino(4cosa3cosa2cosal  equal  to  sino<3cosa2coscxl  results  in  the 
relationship  sina4  "  tano(3 .  The  generalized  recursive 
relationship  is  given  by 

aCn-1)  -  arctanCsinocCnJD  C4-11D 

where  aCn)  and  aCn-1)  are  the  nth  and  Cn-lDth  reduced  tip 
angles,  respectively. 

Equation  4-11  is  useful  in  assessing  the  relationship 
amongst  the  reduced  tip  angles,  but  it  does  not  suggest  any 
specific  values.  The  second  key  point  to  consider  is  the 
fact  that,  since  the  nth  STE  image  is  the  last  of  the 
series,  a  ir/S  rf  pulse  should  be  used  so  as  to  rotate  all 
the  remaining  longitudinal  magnetization  into  the  transverse 
plane.  Thus,  this  inital  Cactually,  final)  condition,  in 
conjunction  with  the  recursive  relation  of  equation  4-11, 
will  yield  all  the  reduced  tip  angles  required.  Listed  in 
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decreasing  order,  the  First  ten  reduced  tip  angles  are;  90 
degrees,  45  degrees,  35.26  degrees,  30  degrees,  26.57 
degrees,  24.09  degrees,  22,21  degrees,  20.70  degrees,  19.47 
degrees,  and  18.43  degrees.  Obviously,  this  list  could  be 
extended  for  any  value  of  n,  although  the  practical 
limitation  is  the  available  signal  to  noise  ratio  of  the  STE 
image.  It  is  interesting  to  note  that  these  values  are 
invariant  to  the  value  of  n.  That  is,  the  five  reduced  tip 
angles  of  a  n  equal  five  series  are  the  same  as  the  last 
five  reduced  tip  angles  in  a,  for  example,  n  equal  fifty 
series.  The  same  recursive  relationship  has  been  developed 
by  Mansfield  CHaB4]  for  use  in  a  different  application. 

The  TART  experiment  was  performed  on  a  phantom  of  seven 
vials  filled  with  water,  and  doped  with  copper  sulfate  to 
various  concentrations.  The  vials  were  2.5  centimeters  in 
diameter,  and  were  stacked  so  that  only  the  wall  thickness 
seperated  them.  The  Tl  values  of  the  vials  are  given  in 
table  4-2,  in  the  following  section.  Uial  seven  contained  a 
gel  solution  to  reduce  its  T2  value.  Thus,  any  unexpected  T2 
dependent  weighting  in  the  TART  sequence  would  result  in 
signal  intensity  variations,  in  the  images  of  vial  seven, 
which  could  not  be  explained  by  Tl  weighting.  The  value  of 
TE  was  15  milliseconds,  and  each  RTi  was  incremented  in 
steps  of  50  milliseconds,  with  an  initial  value  of  50 
milliseconds.  The  entire  pulse  sequence  was  repeated  every 
715  milliseconds. 
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Experiment  seven:  difFusion  coeFFicient  mapping 

In  addition  to  characterizing  an  object  by  its 
relaxation  times,  Tl  and  TB,  and  its  chemical  shiFt,  NflR 
experiments  aFFord  the  opportunity  to  characterize  the 
molecular  translational  selF-diFFusion  coeFFicient  oF  the 
sample.  One  may  utilize  the  Tl  dependence  oF  the  STE  image 
indirectly  to  aid  in  spatially  mapping  the  diFFusion 
coeFFicient.  Recently  reported  methods  oF  in  vivo 
diFFusion  coeFFicient  measurement  CTaBS,  U)eB4D  rely  upon 
the  w/S-T— IT  spin  echo  imaging  technique.  Since  in  vivo.  Tl 
is  typically  greater  than  T2,  the  relative  eFFect  oF 
relaxation  on  signal  intensity  might  be  lessened  by 
observing  the  diFFusional  attenuation  oF  the  STE  rather  than 
that  oF  the  PE  in  the  Tf/S-Tl-Tr/2-T2-Tf/E  sequence,  or  the 
spin  echo  in  the  tt/B-t-tt  sequence.  A  proposed  imaging 
sequence  is  presented  in  Figure  4-9. 

The  modiFication  to  the  basic  STE  imaging  sequence  oF 
Figure  4-1  is  the  addition  oF  diFFusive  pulsed  magnetic 
Field  gradients,  shown  in  Figure  4-9  as  a  part  oF  Gz . 
Assuming  the  inFluence  oF  any  constant  background  Field 
gradient  is  negligible,  the  eFFect  oF  the  diFFusive 
gradients  is  given  by 


2   E      B 
FCG.D.S.A)  -  expC-r  G  D6  CA-S/3)D  C4-ieD 


with  S  and  A  depicted  as  in  Figure  4-9,  where  A  -Al  For  the 
PE  image  and  A  - /fi  For  the  STE  image.  Hence,  the  Tl 
dependence  oF  the  STE  Facilitates  diFFusion  coeFFicient 
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Figure  4-9.  The  stimulated  echo-diffusion  coefficient 
imaging  sequence. 


determination  by  permitting  a  greater  range  in  a  2  values.  In 
sequences  uhere  relaxation  damping  is  strictly  TE  dependent 
Ce.g.  spin  echo  imaging,  primary  echo  imaging),  relaxation 
effects  often  dominate  the  diffusional  effect.  This  is  of 
primary  concern  in  systems  uhich  undergo  restricted 
diffusion  CUIoBBD ,  for  example  living  cells  CTaBBD,  uihen 
it  is  desired  to  extend  the  diffusion  effect  to  the  maximum. 

Results 
Introduction 

The  time  domain  response  of  a  spin  system  to  three  tt/S 
rf  pulses  is  given  in  figure  4-10.  This  is  simply  the 
stimulated  echo  rf  pulse  sequence,  uithout  applied  pulsed 
magnetic  field  gradients,  generated  ujith  the  NT-BO 
spectrometer  and  presented  here  for  illustrative  purposes 
only.  Follouiing  the  notation  of  figure  4-10,  Ca)  is  the 
response  to  the  first  rr/B  rf  pulse  applied  at  time  zero, 
with  a  FID  follouiing  immediately.  A  second  tt/E  rf  pulse  uas 
applied  after  a  time  interval  t1  equal  to  400  milliseconds, 
with  a  FID  following  immediately.  This  response  is  noted  as 
Cb3  in  figure  4-10.  (=1  PE  was  formed  at  time  EtI,  BOO 
milliseconds  from  time  equal  zero,  and  is  indicated  in 
figure  4-10  as  Cc) .  At  a  time  tE  equal  to  B75  milliseconds 
after  the  second  rf  pulse,  a  third  rr/E  rf  pulse  is  applied, 
with  a  FID  following  immediately.  This  response  is  noted  as 
Cd5  in  figure  4-10.  A  STE  was  formed  at  a  time  t1  after  the 
third  rf  pulse,  1475  milliseconds  from  time  equal  zero,  and 
is  noted  in  figure  4-10  as  Ce) .  Three  other  secondary  echoes 
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Figure  H-10.  The  response  to  the  stimulated  echo  sequence. 
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are  formed,  denoted  in  figure  4-10  as  Cf),  Cg),  and  Ch) .  The 
corresponding  times  were,  ErE    C1350  milliseconds  from  time 
equal  zero),  rl+ErE    C1750  milliseconds  from  time  equal 
zero),  and  SCtI+tSD  C2150  milliseconds  from  time  equal 
zero)  . 

One  problem  associated  uith  STE  imaging  is  highlighted 
in  figure  4-10.  It  should  be  noted  that  the  STE  is  not 
temporally  isolated  from  the  secondary  echo  which  formed  at 
time  2t2.  Artifacts  uiould  be  introduced  into  a  STE  image 
based  upon  this  particular  STE.  It  was  possible  to  adjust 
the  values  of  t1  and  tB  so  as  to  isolate  the  STE,  but  this 
approach  would  result  in  restricting  the  timing  parameters 
to  certain  values.  As  mentioned  earlier,  in  practice  a 
spoiling  magnetic  field  gradient  pulse  was  applied  after  the 
formation  of  the  PE,  but  prior  to  the  application  of  the 
third  rf  pulse.  This  spoil  pulse  eliminated  all  transverse 
magnetization  prior  to  the  application  of  the  third  rf 
pulse,  thus  no  secondary  echoes  were  formed. 
Experiment  one:  Tl  contrast  imaging 

Figure  4-11  presents  typical  results  from  a  basic  STE 
imaging  sequence.  The  phantom  shown  in  figure  4-lla  is  a  red 
Delicious  apple  with  a  one  centimeter  diameter  vial  of 
copper  sulfate  doped  distilled  water  embedded  in  its  center, 
dubbed  the  "William  Tell  phantom."  Figure  4-llb  is  the  PE 
image  CtI  -  20  milliseconds),  and  figure  4-llc  is  the  STE 
image  CtI  -  20  milliseconds,  t2  -  586  milliseconds). 
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Figure  4-11.  The  William  Tell  phantom.  Ca)  au  natural .  Cb) 
the  PE  image,  Cc)  and  the  STE  image. 
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Oiuing  to  a  short  t1  value,  the  doped  water  and  apple 
flesh  did  not  vary  greatly  in  intensity  in  the  PE  image. 
Most  of  the  variation  in  intensity  within  the  image  was  due 
to  volume  averaging,  except  for  the  region  of  low  intensity 
directly  beneath  the  water  vial .  This  appeared  to  be  the 
region  of  low  flesh  density  where  the  seeds  reside.  The 
source  of  high  intensity  at  the  base  of  the  apple  was  not 
known  with  certainty,  but  could  be  the  result  of  a  bruise. 

The  Tl  of  the  doped  water  solution  was  determined  to  be 
513±17  milliseconds  by  an  inversion  recovery  sequence 
utilizing  a  composite  tt  pulse  without  applied  gradients.  In 
a  similar  fashion  the  Tl  of  apple  flesh  was  determined  to  be 
less  than  SOO  milliseconds.  This  spread  in  Tl  values  was 
dramatically  demonstrated  in  figure  4-llc  by  the  STE  image. 
Experiment  two:  quantitative  Tl  imaginn 

In  addition  to  Tl  contrast  imaging,  STE  imaging 
permitted  quantitative  measurements  of  Tl  values  as  outlined 
earlier.  Figure  4-12  depicts  the  results  of  quantitative  Tl 
imaging  with  stimulated  echoes.  The  phantom  consisted  of 
three  one  centimeter  diameter  vials  of  copper  sulfate  doped 
water.  The  top  right  vial  had  a  Tl  of  250±9  milliseconds, 
the  bottom  left  vial  had  a  Tl  of  399±3  milliseconds,  and  the 
bottom  right  vial  had  a  Tl  of  512117  milliseconds.  These  Tl 
values  were  determined  individually  by  an  inversion  recovery 
sequence  utilizing  a  composite  tt  pulse,  without  applied 
magnetic  field  gradients. 
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Figure  4-12.  The  quantitative  use  of  the  STE  image 
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Figure  4-12a  is  the  PE  imagB  CtI  -  20  tnillisBconds) , 
uihilB  figures  4-12b  thru  4-lEf  are  STE  images  of  the  same 
phantom,  each  acquired  with  a  different  t2  value  CtI  -  20 
millisBConds,  t2  -  100,  250,  500,  1000,  and  1500 
milliseconds  respBctivBly ) .  Ths  Tl  valus  for  the  bottom  left 
vial  was  determinsd  in  two  mays,  bassd  upon  thB  maximum 
signal  amplitudes  of  the  imags  of  that  vial  in  each  of  the 
STE  images.  The  first  method,  plotting  the  logarithm  of  the 
STE's  maximum  amplitude  versus  t2,  yielded  a  value  for  Tl  of 
3B2±2  milliseconds.  This  calculated  value  represented  less 
than  a  five  percent  relative  error  when  compared  to  the 
actual  Tl  value.  As  outlined  earlier,  a  more  precise 
calculated  value  of  Tl  could  be  determined  if  the  effects  of 
diffusion  were  considered,  in  the  manner  of  equation  4-B. 
When  the  data  uere  analyzed  in  the  manner  specified  in  the 
text  folloLjing  equation  4-B,  a  calculated  value  for  Tl  of 
3aB±2  milliseconds  luas  obtained,  uihich  represented  less  than 
a  three  percent  relative  error  when  compared  uith  the  actual 
Tl  value.  These  results  uere  uithin  acceptable  experimental 
limits. 
Experiment  three:  uater-lipid  contrast  imaging 

Figure  4-13  is  a  uiater-lipid  Tl  contrast  image  of  a 
hen's  egg,  by  STE  imaging.  Figure  4-13b  is  the  PE  image  CtI 
-  20  milliseconds),  and  figure  4-13c  is  the  STE  image  CtI 
-20  milliseconds,  t2  -  186  milliseconds).  Consider  figure 
4-13c,  where  the  high  intensity  signal  is  from  the  egg's 
albumen  Cthe  water  dense  element),  the  circular  region  of 
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Figure  4-13.  A  luater-lipid  image  of  a  hen's  egg.  Ca)  au 
naturel .  Cb)  the  PE  image,  Cc)  and  the  STE 
image . 
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low  intensity  is  from  the  egg's  yolk  (the  lipid  dense 
element),  and  the  "missing"  semicircular  portion  near  the 
bottom  is  the  egg's  air  sac.  The  magnitude  of  the  readout 
gradient  used  in  both  figures  H-lBb  and  1-13c  was  0.5  mT/m, 
which  was  insufficient  to  overcome  the  lipid  chemical  shift 
from  the  yolk  CDiB4D,  most  apparent  in  the  PE  image.  This 
shift  appears  in  the  readout  direction  C vertical)  of  figure 
4-13b,  but  is  suppressed  in  figure  4-13c  due  to  the 
relatively  short  Tl  of  the  lipid  component  of  the  yolk. 
(Additionally ,  figure  4-13b  suffers  from  a  distorted  aspect 
ratio,  owing  to  the  influence  of  residual  gradients.  Hence, 
the  phase  encoding  direction  Chorizontal)  is  reduced  in 
figure  'i-13b  relative  to  figure  4-13c. 
Experiment  four:  chemical  shift  imaging 

Figure  4-14  is  an  example  of  STE  chemical  shift 
imaging.  The  phantom  used  in  this  experiment,  depicted  in 
figure  4-14a,  consisted  of  two,  one-centimeter  diameter 
vials  separated  spatially  in  the  y  direction  but  not  the  x 
direction.  One  vial  was  filled  with  vegetable  oil  having  a 
Tl  of  S04±2  milliseconds,  as  determined  by  an  inversion 
recovery  experiment  with  a  composite  rr  pulse,  without 
applied  magnetic  field  gradients.  The  second  vial  contained 
a  mixture  of  distilled  water  and  deuterated  water,  having  a 
Tl  of  S.a±0.1  seconds,  determined  in  a  like  manner.  The 
deuterated  water  was  added  to  reduce  the  signal  amplitude, 
so  that  the  unsaturated  signal  intensities  of  the  two  vials 
were  similar. 
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Figure  4-14 


A  STE  chemical  shift  image.  Ca)  The  phantom, 
Cb)  the  PE  image,  Cc)  and  the  STE  image. 
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To  retain  chemical  shift  information  Calong  the  cs  axis 
as  depicted  in  figure  4-14D,  no  readout  gradient  was 
employed.  The  x  direction  was  phase  encoded  (depicted  as  the 
X  axis  in  figure  4-14).  The  images  are  displayed  as  stacked 
plots  of  spectra.  The  PE  image  CtI  -  20  milliseconds), 
figure  4-14b,  demonstrates  the  near  equal  intensities  of  the 
two  vials  for  short  t1,  and  the  3.5  ppm  relative  chemical 
shift.  In  the  STE  image  CtI  -  EO  milliseconds,  t2  -  250 
milliseconds)  of  figure  4-14c,  the  oil  image  is  suppressed 
due  to  the  large  difference  in  Tl  between  the  water  and  oil. 
Each  image  was  scaled  to  maximum  peak  intensity . 
Experiment  five:  extended  STE  multiecho  imaging 

Figure  4-15  illustrates  the  results  of  the  extended  STE 
multiecho  imaging  sequence,  as  given  in  figure  4-B.  The 
phantom  consisted  of  three,  one  centimeter  diameter  vials  of 
copper  sulfate  doped  water.  The  top  right  vial  had  a  Tl  of 
E50±9  milliseconds,  the  bottom  left  vial  had  a  Tl  value  of 
339±3  milliseconds,  and  the  bottom  right  vial  had  a  Tl  of 
51E±17  milliseconds.  All  Tl  values  were  determined  by  an 
inversion  recovery  experiment  utilizing  a  composite  tt  pulse, 
without  applied  magnetic  field  gradients. 

Figure  4-15a  is  the  PE  image  CtI  -  EO  milliseconds), 
figure  4-15b  is  the  STE  image  CtI  -  EO  milliseconds,  tE  - 
148  milliseconds),  and  figure  4-15c  is  the  TE  image  CtI  -  EO 
milliseconds,  tE  -  14B  milliseconds,  and  t3  -  548 
milliseconds).  Note  the  Tl  contrast  imposed  on  the  STE  image 
and  the  TE  image  as  a  function  of  tS  and  t3,  respectively. 
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Figure  4-15.  An  extended  STE  multiecho  series  of  images.  Ca) 
The  PE  image,  Cb)  the  STE  image,  Cc)  and  the  TE 
image . 
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The  aspect  ratio  of  the  STE  image  was  one,  For  residual 
gradient  influence  along  the  phase  encoding  direction 
ChorizontalD  was  negligible  for  this  echo.  Both  the  PE  image 
and  the  TE  image  ujere  influenced  along  the  phase  encoding 
direction  by  residual  gradients.  The  amount  of  influence  uas 
a  function  of  houi  much  time  the  magnetization,  which  was 
used  in  echo  formation,  spent  in  the  transverse  direction 
after  initial  phase  encoding  within  the  t1  interval,  and  the 
magnitude  of  the  residual  gradient  during  this  time.  The 
influence  was  greatest  on  the  TE  image,  for  although  the  PE 
spent  a  time  t1  Cequal  to  SO  milliseconds!)  in  the  transverse 
direction  while  in  the  presence  of  a  residual  gradient,  the 
TE  spent  a  time  t2  Cequal  to  148  milliseconds)  in  the 
transverse  direction  in  the  presence  of  the  residual 
gradient.  Indeed,  the  time  spent  by  the  TE  overlapped  that 
of  the  PE .  All  images  were  scaled  identically. 
Experiment  six:  TflRT  multiecho  imaging 

Figure  4-lB  is  a  comparison  of  the  images  obtained  with 
the  conventional  spin  echo  imaging  sequence,  and  that 
obtained  from  the  PE  of  the  TART  sequence.  The  images  are 
presented  as  signal  intensity  contour  plots.  As  predicted  by 
equation  4-3,  the  signal  to  noise  ratio  of  the  PE  image  is 
reduced  by  a  factor  of  two  compared  to  that  of  the  spin  echo 
image.  The  images  in  figure  4-lB,  and  all  other  images 
presented  which  were  generated  with  the  TART  sequence,  were 
phase  encoded  in  the  x  direction  Chorizontal ) ,  and  with  the 
read  out  gradient  applied  along  the  y  direction  (vertical). 
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Figure  4-15.  A  comparison  of  the  spin  echo  image  and  the 
primary  echo  image.  Ca)  The  spin  echo  image, 
CbD  and  the  primary  echo  image. 
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The  Tl  UBighted  series  of  STE  images  generated  by  the 
TART  sequence  are  depicted  in  figure  4-17.  The  images  are 
presented  as  signal  intensity  contour  plots.  The  recovery 
time,  RT,  progressed  from  50  milliseconds  in  figure  4-17a, 
to  BOO  milliseconds  in  figure  4-17d,  in  50  millisecond 
increments . 

Qualitatively,  the  Tl  ueighted  series  of  STE  images 
represented  the  available  contrast  difference  due  to 
variations  in  spin-lattice  relaxation  rates.  The 
quantitative  nature  of  the  intensity  variations  in  the  TART 
imaging  series  is  represented  in  table  4-E.  The  first  column 
of  the  table  gives  the  Tl  value  determined  for  each  vial  by 
an  inversion  recovery  experiment,  uhere  eight  data  points 
were  fitted  to  a  three  parameter  exponential  curve.  Each 
vial  uas  measured  separately  by  positioning  the  vial  at  the 
center  of  the  rf  coil,  and  operating  the  CSI  unit  as  a 
spBctromatBr,  uithout  applied  magnetic  field  gradients. 
Column  tuo  is  simply  the  standard  error  in  these  Tl  values 
estimated  with  the  inversion  recovBry  tBchniqus. 

Column  three  of  table  4-S  gives  the  Tl  values 
determined  for  each  vial  by  fitting  the  signal  intensities 
from  the  four  STE  images  to  a  two  parameter  exponential 
curve.  The  intensity  values  were  determined  by  averaging 
pixels  over  a  fixed  size  region  within  the  image  boundaries 
of  each  vial .  Column  four  is  simply  the  standard  error  in 
these  Tl  values  estimated  with  the  TART  imaging  sequence. 
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Figure    4-17.    A    series    of    TART    images, 
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Table  4-S 


Calculated  Tl  Ualues  From  Phantom  in  TART  Experiment 


Tl  Ualues  Cmsec) 


Uial 

IR 

Standard 
error 

TART 

Standard 
error 

1 

115 

±0.7 

103 

±2.9 

2 

63 

±0.09 

66 

±1.2 

3 

58 

±0.04 

61 

±0.9 

H 

63 

±0.02 

65 

±2.3 

5 

115 

±0.7 

117 

±7.0 

E 

64 

±0.2 

63 

±1.6 

7 

246 

±1.7 

227 

±9.8 

: 
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Experiment  seven:  diffusion  coeFFiclent  mapping 

As  indicated  in  figure  4-4,  there  was  an  appreciable 
residual  gradient  present  after  the  application  of  a  pulsed 
magnetic  field  gradient.  Its  influence  on  image  quality  luas 
predominately  qualitative,  resulting  in  nonunity  aspect 
ratios.  Unfortunately,  the  residual  gradient  was  sufficient 
to  render  quantitative  diffusion  coefficient  mapping 
impractical,  although  an  example  of  qualitative  diffusion 
coefficient  contrast  imaging  has  recently  been  presented  by 
others  CfleBSD,  which  utilized  an  imaging  sequence 
identical  to  that  in  figure  4-9. 

Discussion 

Exploiting  stimulated  echoes  in  NtlR  imaging  involved 
the  judicious  use  of  the  ability  to  store  and  retrieve 
information  along  the  longitudinal  direction,  afforded  by 
the  stimulated  echo.  The  results  of  this  investigation  took 
the  Form  of  specific  examples  which  prudently  used  this 
unique  property  of  the  STE  in  NflR  imaging.  First,  it  was 
shown  that  in  addition  to  generating  Tl  contrast  images,  it 
was  possible  to  calculate  quantitative  Tl  information  from  a 
series  of  STE  images  in  which  the  storage  time  had  been 
systematically  varied.  Second,  a  novel  application  of  the  Tl 
weighted  STE  image  was  accessed;  the  enhancement  or 
suppression  of  elements  in  the  object  with  different  Tl 
values.  Third,  it  was  demonstrated  that  the  STE  was  easily 
integrated  into  chemical  shift  imaging  schemes.  Fourth,  two 
STE  methods  were  developed  which  permitted  the  acquisition 
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of  a  series  of  STE  images  within  one  imaging  sequence,  where 
each  image  was  prcgressively  weighted  by  increasing  Tl 
relaxation  clamping.  Finally,  a  method  of  in  vivo 
determination  of  molecular  translational  self-diffusion 
coefficients,  which  utilized  the  STE  to  lessen  the  effect  of 
T2  relaxation,  was  proposed.  A  detailed  disscussion  of  each 
of  these  results  follows. 
Experiment  One:  Tl  Contrast  Imaging 

The  PE  and  STE  images  of  the  William  Tell  phantom 
illustrated  in  figure  4-11  indicate  that  stimulated  echo 
imaging  not  only  works,  but  also  works  well  to  produce  Tl 
contrast  images.  It  is  worthwhile  to  elaborate  on  the 
unexpected  high  intensity  regions  in  the  images.  In  both  the 
PE  and  STE  image,  there  was  a  high  intensity  region  at  the 
base  of  the  apple,  hypothesized  to  be  a  bruise.  Also,  there 
is  a  region  of  high  intensity  surrounding  the  position  where 
the  cylindrical  vial  was  located.  Since  the  resolution  of 
the  image  was  1.5  X  0.75  millimeters,  and  since  edges  were 
well  defined  in  spite  of  volume  averaging,  it  seems  likely 
that  the  region  of  high  intensity  surrounding  the  vial  was 
not  a  structual  component  of  the  apple,  and  it  is  suggested 
to  be  associated  with  apple  flesh  which  was  bruised  while 
forcing  the  vial  through  the  apple.  The  reason  bruised  apple 
flash  would  result  in  increased  signal  intensity  is 
explained  by  the  following. 

The  full  wrath  of  the  Bloembergen,  Purcell  and  Pound 
CB14BD  theory  of  relaxation  need  not  be  brought  to  bear 
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against  this  problom,  rather,  Just  one  simple  result  of  the 
theory .  The  result  is  that  the  spin-lattice  relaxation  rate 
is  proportional  to  the  correlation  time  tc,  in  the  extreme 
narrouing  situation.  This  result,  coupled  uith  the  fact  that 
a  more  mobile  solution  has  a  lower  tc,  suggests  a  reason  why 
bruised  apple  flesh  resulted  in  increased  signal  intensity. 

An  apple  is  predominately  made  up  of  water,  although 
the  flesh  is  quite  cellular  in  nature,  offering  a  natural 
restriction  to  molecular  mobility.  When  the  apple  flesh  is 
bruised,  that  is,  when  the  cellular  structure  is  broken 
down,  the  molecular  mobility  increases.  This  increased 
mobility  results  in  a  lower  tc,  and  in  turn  results  in  an 
increase  Tl  value.  Hence,  bruised  regions  have  longer  Tl 
values,  and  would  therefore  appear  increased  in  intensity 
over  the  surrounding  unbruised  regions. 

Although  the  STE  in  figure  4-llc  admirably  displayed  Tl 
contrast,  the  relaxation  damping  of  the  stimulated  echo  is 
not  solely  determined  by  Tl ,  for  there  is  a  expC-ETl/T2D 
dependence  also.  To  obtain  a  "pure"  Tl  mapped  image,  simple 
post  acquisition  processing  is  required.  It  may  be  noted 
that  the  ratio  of  the  STE's  maximum,  given  by  equation  4-7, 
to  the  PE's  maximum,  given  by  equation  4-3,  MCSTE/PE),  is 
given  by 


nCSTE/PE)  -  Cl/e)Csin92sine3/sin  Cee/E)DexpC-Te/Tl )  • 

rCG,D,Tl,T25  C4-13: 
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where  all  the  terms  have  been  previously  defined,  with 
rCG,D,Tl,T2)  being  the  ratio  of  the  diffusional  damping  of 
the  STE  to  that  of  the  PE .  For  ideal  tt/E  rf  pulses  and 
negligible  diffusional  damping,  MCSTE/PE)  equals 
expC-T2/Tl),  assuming  the  pulse  sequence  repetition  time  is 
sufficiently  long  so  as  to  avoid  additional  Tl  dependence 
via  saturation.  Hence,  images  with  strictly  Tl  contrast 
could  readily  be  generated  from  a  single  pass  of  the  STE 
imaging  sequence  and  appropriate  post  acquisition  data 
manipulation.  Furthermore,  it  is  noted  that  fICSTE/PE)  is 
functionally  independent  of  91,  making  it  insensitive  to 
changes  in  91 . 
Experiment  Two:  Quantitative  Tl  Imaging 

The  results  of  this  experiment  indicate  that  STE 
imaging  is  not  only  capable  of  Tl  contrast  imaging,  but  that 
the  STE  imaging  sequence  may  be  as  quantitative  as  the 
conventional  inversion  recovery  or  spin  echo  techniques. 
Indeed,  STE  imaging  has  some  advantages  over  the 
conventional  techniques  which  yield  calculated  Tl  values. 

The  conventional  inversion  recovery  imaging  sequence 
introduces  Tl  weighting  into  the  image  by  inverting  the 
equilibrium  magnetization,  taken  to  be  positive  initially, 
along  the  negative  z  axis.  The  magnetization  is  sample  by  a 
SO  degree  rf  pulse  as  it  relaxes  along  the  z  axis  from  the 
negative  to  the  positive.  Conventional  state  of  the  art  NflR 
imagers  handle  this  situation  in  one  of  two  ways.  Some 
instruments  offer  the  capability  of  retaining  this  signed 
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information  by  performing  a  phase  sensitive  type  of  image 
reconstruction.  This  procedure  is,  at  best,  difficult  to 
optimize  when  the  object  imaged  is  so  diverse  as  the  human 
body .  The  second  option  offered  is  a  simple  magnitude 
reconstruction  of  the  image,  where  all  signed  information  is 
lost.  This  approach,  although  simple  to  implement,  could 
result  in  confusion  when  regions  of  the  object  which  are 
equal  in  intensity  but  opposite  in  sign  are  displayed  as 
identical.  Since  in  STE  imaging  all  Tl  weighting  is  one 
sided  Ce.g.  all  positive),  the  simple  magnitude 
reconstruction  suffices  to  display  all  available  Tl 
information  without  ambiguity . 
Experiment  Three:  Uater-Lipid  Contrast  Imaging 

The  images  of  figure  1-13  illustrate  a  unique  feature 
of  STE  imaging,  the  ability  to  suppress  the  shorter  Tl 
elements  in  an  object  by  Tl  weighting. 

Conventional  SE  imaging  techniques  introduce  Tl 
weighting  by  saturating  the  signal.  That  is,  the  sequence  is 
repeated  so  rapidly  that  there  is  little  time  between 
repetitions  for  Tl  relaxation.  Thus,  those  elements  with 
short  Tl  values  are  able  to  return  to  their  equilibrium 
value  prior  to  the  next  repetition  of  the  sequence,  whereas 
those  elements  with  longer  Tl  values  are  destined.  Just  as 
Sisyphus  was  CHaSBD,  to  never  fully  regain  their 
equilibrium  value  before  the  pulse  sequence  starts  over. 
Thus,  elements  with  longer  Tl  values  are  suppresed ,  whereas 
elements  with  shorter  Tl  values  are  seemingly  enhanced. 
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In  contrast,  the  simple  expC-rE/Tl)  uieighting  of  the 
STE  image  results  in  the  suppression  of  shorter  Tl  elements, 
while  longer  Tl  elements  are  seemingly  enhanced.  For 
example,  consider  the  case  where  t2  is  150  milliseconds,  and 
two  different  elements  exist  in  the  object,  one  with  Tl 
equal  to  100  milliseconds,  and  the  other  with  a  Tl  of  1000 
milliseconds.  Based  solely  upon  this  difference  in  Tl 
values,  with  all  else  being  equal,  the  longer  Tl  element 
would  appear  nearly  four  times  as  intense  as  the  shorter  Tl 
element.  The  clinical  utility  of  this  characteristic  of  STE 
imaging  has  yet  to  be  explored.  However,  since  many 
abnormalities,  and  in  particular  malignancies,  have  longer 
Tl  values,  this  raises  tantalizing  possibilities. 
Experiment  Four:  Chemical  Shift  Imaging 

It  has  been  demonstrated,  in  figure  H-14,  that  the 
stimulated  echo  imaging  technique  is  sufficently  flexible  so 
as  to  permit  its  incorporation  into  any  existing  chemical 
shift  imaging  scheme  which  relies  upon  spin  echo  formation. 
Not  only  may  the  basic  stimulated  echo  sequence  be  used  in 
chemical  shift  imaging,  but  so  may  either  of  the  two 
multiecho  stimulated  echo  sequences  presented  in  this 
investigation.  Additionally,  quantitative  analysis  of  Tl 
could  easily  be  done  for  each  of  the  chemically  shifted 
components  as  outlined  previously.  Extension  of  this 
technique  to  nuclei  other  than  hydrogen  would  be  beneficial, 
for  example,  phosphorus,  where  Tl  is  often  much  longer  than 
TS  CHo75D . 
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Experiment  Five:  Extended  STE  flultlecho  Imaning 

In  contrast  with  the  basic  STE  sequence,  the  extended 
STE  multiecho  sequence  allouis  for  greater  Tl  contrast  for  a 
given  imaging  time,  and  offers  a  convenient  method  of  Tl 
estimation.  The  spin  echo  multipulse  sequence  affords 
similar  advantages  for  TE  contrast  imaging.  The  Tl  weighted 
multiecho  imaging  scheme  has  the  added  advantage  that,  in 
addition  to  the  train  of  Tl  weighted  images,  the  TE  weighted 
PE  image  is  also  produced.  Thus,  this  single  sequence  yields 
a  TS  dominated  image,  and  a  series  of  Tl  weighted  images, 
all  in  spatial  registration.  This  would  be  quite  useful  in 
medical  applications,  where  typically  a  spin  echo  image  with 
a  short  echo  time  is  useful  for  anatomical  differentiation, 
and  images  with  high  Tl  contrast  are  useful  in  pathological 
differentiation . 
Experiment  Six:  TftRT  flultiecho  Imaging 

The  TART  imaging  sequence  provides  both  a  series  of  STE 
images  progressively  weighted  with  Tl  relaxation  and 
constant  T2  weighting,  and  also  a  TE  dominated  PE  image. 
Both  the  spatial  resolution  and  the  total  imaging  times  are 
comparable  to  the  conventional  NflR  imaging  techniques. 
Although  clinical  utility  has  yet  to  be  realized,  the 
flexibility  of  the  technique  promises  wide  applications.  For 
all  the  virtues  of  TART  imaging,  there  is  one  major  vice, 
the  reduction  in  the  signal  to  noise  ratio  of  the  TART 
images  over  the  corresponding  spin  echo  image. 
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NriR  imaging  makes  available  a  myriad  of  information 
concerning  the  object  under  investigation.  Unfortunately, 
there  is  aluiays  a  trade  off  involved,  the  most  common  one 
being  the  trade  off  betuieen  derived  information  on  the  one 
hand,  and  imaging  time  and  signal  to  noise  ratio  on  the 
other.  Typically,  if  one  uiants  to  glean  additional 
information  from  a  spin  system,  either  additional  time  must 
be  spent  or  a  reduction  in  the  signal  to  noise  ratio  must  be 
tolerated.  This  is  precisely  the  situation  in  TART  imaging. 
Additional  Tl  information  is  obtained  through  the  series  of 
STE  images,  but  only  by  paying  the  price  of  a  reduced  signal 
to  noise  ratio. 

The  signal  to  noise  ratio  of  the  PE  and  the  series  of 
STE  images  may  be  quantified  and  compared  to  the  signal  to 
noise  ratio  of  the  conventional  spin  echo  image.  If  the 
equilbrium  magnetization  is  fli  ,  then,  ignoring  relaxation 
and  diffusional  damping,  the  maximum  amplitude  of  the  spin 
echo  image  is  simply  fli  .  This  is  the  ruler  by  which  we 
measure  the  PE  and  STE  images  of  the  TART  sequence. 

As  indicated  in  equation  4-3,  the  maximum  amplitude  of 
the  PE  image,  ignoring  relaxation  and  diffusional  damping 
and  assuming  ideal  tt/E  rf  pulses,  is  C1/2)I1I  .  Thus,  the  PE 
image,  uhich  has  the  same  relaxation  ueighting  as  the 
conventional  spin  echo  image,  has  a  signal  to  noise  ratio 
reduced  by  a  factor  of  two.  The  other  half  of  the  signal  to 
noise  is  available  for  STE  image  production. 
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If  only  a  single  STE  image  is  produced,  then  its  signal 
to  noise  ratio  is  simply  a  half  of  that  for  the  spin  echo 
image.  If  the  stored  longitudinal  magnetization  is  used  to 
create  a  series  of  STE  images,  as  in  the  TART  sequence,  then 
the  available  Cl/SDfli  is  vectorially  divided  amongst  each 
STE  image.  Since,  in  the  TART  imaging  sequence,  each  STE  has 
the  identical  tip  angle  dependence,  this  factor  is  properly 
given  by  sinal ,  simply  the  tip  angle  dependence  of  the  first 
STE  image  in  the  series.  Recall,  the  vector  sum  of  all  the 
STE  images  in  the  series,  neglecting  relaxation  and 
diffusion  damping,  must  be  equal  to  Cl/E)f1i  .  For  a  series  of 
n  images,  the  vector  sum  of  all  the  STE  image's  maximum 
amplitudes,  EM,  is  given  by 
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uihich  may  be  simplified,  as  noted,  by  realizing  each  reduced 
tip  angle  term  equals  the  first  reduced  tip  angle  term, 
sinal.  Thus,  Efl  equals  the  square  root  of  nC  Cl/E)ni  sinalD  . 
Since  Efl  must  be  equal  to  Cl/EDni  ,  sinal  must  be  inversely 
equal  to  the  square  root  of  n.  Thus,  since  sinal  is  the 
reduced  tip  angle  factor  for  all  the  STE  images  in  a  TART 
sequence  with  n  images,  it  may  be  said  that  the  signal  to 
noise  ratio  of  each  image  is  reduced  by  a  factor  of  the 
square  root  of  n .  A  similar  result  was  expressed  by 
Mansfield  CnaB4D  for  a  different  application. 


CHAPTER  U 
SUnriARY  AND  CONCLUSIONS 


The  emphasis  of  this  uork  uias  tuo-fold.  One,  the 
accurate  and  precise  determination  of  the  spin-lattice 
relaxation  time  uith  a  conventional  state-of-the-art  NFIR 
imager,  and  two,  the  application  of  the  stimulated  echo 
CSTED,  uith  its  unique  Tl  dependence,  to  NflR  imaging. 

It  luas  shoun  that  under  certain  conditions  it  is 
possible  to  obtain  an  accurate  and  precise  estimate  of  Tl 
ujith  a  commercial  NflR  imager.  The  most  reliable  method 
involved  the  acquisition  of  three  fast  inversion  recovery 
images,  with  suitably  chosen  values  of  the  inverting  time  as 
given  in  equation  3-3.  If  the  signal  intensities  from  these 
images  are  fitted  to  a  three-parameter  monoexponential 
function,  then  an  estimated  Tl  value  was  determined  which 
was  insensitive  to  whether  the  object  was  in  an 
inhomogeneous  rf  field  or  not.  The  source  of  the  field 
inhomogeneity  may  be  the  rf  coil  itself,  or  the  rf 
attenuating  properties  of  the  human  body  which  act  to 
perturb  the  transmitted  rf  field. 

Additionally,  it  was  also  shown  that  three  spin  echo 
images,  with  suitably  chosen  values  of  the  pulse  sequence 
repetition  time  as  given  in  equation  3-5,  could  be  used  for 
Tl  estimation.  In  this  case,  the  signal  intensities 
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luould  bs  fitted  to  a  tuio-parameter  monoBxponential  function. 
Obtaining  estimated  Tl  values  in  this  fashion  could  be  more 
precise  and  time  efficient  than  the  fast  inversion  recovery 
technique,  but  utas  aluiays  much  less  accurate. 

The  results  of  this  uork  had  additional  value  in  that 
they  could  be  applied  a  priori  to  aid  in  designing  an 
entire  imaging  scheme  which  would  result  in  reliable  Tl 
estimates,  or  the  results  could  be  used  in  an  a  posteriori 
fashion,  to  evaluate  the  reliability  of  a  calculated  Tl 
value . 

Although  wide  in  scope  and  exacting  in  execution,  this 
investigation  did  have  some  necessary  errors  of  omission.  By 
limiting  the  forms  of  the  fitting  functions  to  those 
suggested  by  the  NMR  spectroscopy  literature,  certain 
refinements  unique  to  imaging  were  not  considered  CPyB3D . 
Additionally,  two-point  estimates  of  Tl  values  based  upon 
mixtures  of  inverting  and  noninverting  techniques  were  not 
considered  CLiBBD  .  Never — the-less,  the  conclusions  of  this 
investigation  are  correct  and  useful,  as  long  as  they  are 
utilized  within  the  context  of  the  study. 

This  work  was  also  unique  for  it  introduced  the  use  of 
the  stimulated  echo  for  data  acquisition  in  NMR  imaging. 
First,  it  was  shown  that  in  addition  to  generating  Tl 
contrast  images,  it  was  passible  to  calculate  quantitative 
Tl  information  from  a  series  of  stimulated  echo  images,  in 
which  the  storage  time  had  been  systematically  varied.  Also, 
a  novel  application  of  the  Tl  weighted  STE  was  demonstrated: 
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the  enhancement  or  suppression  of  elements  in  the  object 
with  different  Tl  values.  Additionally,  it  was  demonstrated 
that  the  stimulated  echo  could  be  integrated  into  chemical 
shift  imaging  schemes.  Perhaps  of  most  significance  were  the 
two  STE  imaging  methods  presented  which  permitted  the 
acquisition  of  a  series  of  STE  images  within  one  imaging 
sequence,  where  each  image  is  progressively  weighted  with 
increasing  Tl  relaxation  damping.  Stimulated  echo  imaging 
was  also  shown  to  have  the  potential  to  be  a  method  of  in 
vivo  determination  of  molecular  self-diffusion 
coefficients . 

STE  imaging  has  certain  characteristics  which  makes  it 
a  welcomed  addition  to  NfIR  imaging.  Since  all  STE  sequences 
use  tip  angles  of  90  degree  or  less,  the  instantaneous  and 
integrated  rf  power  deposition  in  the  human  body  is  less 
than  in  comparable  SE  imaging  techniques,  with  no  loss  in 
the  net  signal-to-noise  ratio.  This  feature  indicates  that 
STE  imaging  techniques  might  be  ideal  for  imaging  at  higher 
frequencies,  where  the  body  is  a  more  efficient  antenna. 

Since  the  STE  may  be  used  wherever  the  spin  echo  is 
currently  used  in  NfIR  imaging,  the  spatial  resolution  and 
the  total  imaging  time  may  be  the  same  as  for  a  conventional 
image.  Additionally,  STE  imaging  lends  itself  to  multislice 
imaging.  Thus,  multislice  and  Tl  weighted  multiecho  STE 
images  could  be  acquired  much  as  multislice  and  multiecho 
spin  echo  images  are  currently  acquired. 
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Bath  of  these  investigations,  although  different  in 
emphasis,  were  formulated  based  upon  the  same  fundamental 
premise,  that  it  was  possible  to  use  the  microscopic 
properties  of  matter  as  investigative  tools  with  which  to 
examine  the  macroscopic  properties  of  matter.  Indeed,  it  is 
Just  this  fact,  that  the  phenomenon  of  NMR  is  based  on 
interactions  with  the  atomic  nucleus,  which  gives  NUR 
imaging  the  potential  to  be  one  of  the  most  important  probes 
yet  devised  with  which  to  examine  matter.  Other 
investigative  techniques  which  rely  upon  interactions  with 
the  electron  energy  levels  of  an  atom  have  the  potential  to 
alter  the  chemical,  or  biochemical  environment  of  the  object 
under  investigation.  NflR  imaging  is  unique  in  that  it  only 
slightly  perturbs  the  nuclear  energy  levels  of  an  atom, 
although  the  price  paid  is  that  the  information  gleaned  is 
done  so  with  relatively  poor  sensitivity.  It  is  thought  that 
the  conclusions  of  these  investgations  into  spin-lattice 
relaxation  in  NtlR  imaging  have  contributed  to  a  more 
efficient  use  of  the  information  which  NflR  imaging  is  able 
to  supply. 
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